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Abstract 

Boundaries occur naturally in kinetic equations and boundary effects 
are crucial for dynamics of dilute gases governed by the Boltzmann equa- 
tion. We develop a mathematical theory to study the time decay and 
continuity of Boltzmann solutions for four basic types of boundary con- 
ditions: inflow, bounce-back reflection, specular reflection, and diffuse 
reflection. We establish exponential decay in L°° norm for hard poten- 
tials for general classes of smooth domains near an absolute Maxwellian. 
Moreover, in convex domains, we also establish continuity for these Boltz- 
mann solutions away from the grazing set of the velocity at the boundary. 
Our contribution is based on a new decay theory and its interplay with 
delicate L°° decay analysis for the linearized Boltzmann equation, in the 
presence of many repeated interactions with the boundary. 

1 Introduction 

Boundary effect plays a crucial role in the dynamics of gases governed by the 
Boltzmann equation: 

dtF + vV^F^Q{F,F) (1) 

where F{t, x, v) is the distribution function for the gas particles at time t > 0, 
position X G fl, and v G R'^. Throughout this paper, the collision operator 
takes the form 

Q(Fi,F2) = / / \v-u\^Fiiu')F2{v')qo{9)dujdu 
-II \v - upFi{u)F2iv)qoi9)dujdu 

= QgainiPl, F2) ~ QlossiFl, F2), (2) 

where u' — u+{v — u)-uj, v' = v—{v~u)-uj, cos 9 = {u — v)-lu/\u — v\, < 7 < 1 
(hard potential) and < q(){9) < C\ cos9\ (angular cutoff). The mathematical 
study of the particle-boundary interaction in a bounded domain and its effect 



on the global dynamics is one of the fundamental problems in the Boltzmann 
theory. There are four basic types of boundary conditions for F{t, x, v) at the 
boundary dO, : (1) In-flow injection: in which the incoming particles are pre- 
scribed; (2) bounce-back reflection: in which the particles boimcc back at the 
reverse the velocity; (3) specular reflection: in which the particle bounce back 
specularly; (4) diffuse reflection (stochastic): in which the incoming particles 
arc a probability average of the outgoing particles. Due to its importance, there 
have been many contributions in the mathematical study of different aspects of 
the Boltzmann boundary value problems [A], [AC], [AEMN], [AEP], [AH], [C2], 
[C3], [CC], [Dc], [Gui], [H], [LY], [MS], [M], [US], [YZ], among others, see also 
references in the books [CI], [CIP] and [Ul]. 

According to Grad (p243, [Grl]), one of the basic problems in the Boltz- 
mann study is to prove existence and uniqueness of its solutions, as well as their 
time-decay toward an absolute Maxwellian, in the presence of compatible phys- 
ical boundary conditions in a general domain. In spite of those contributions to 
the study of Boltzmann boundary problems, there are few mathematical results 
of uniqueness, regularity, and time decay-rate for Boltzmann solutions toward 
an Maxwellian. In [SA], it was announced that Boltzmann solutions near a 
Maxwellian would decay exponentially to it in a smooth bounded convex do- 
main with specular reflection boundary conditions. Unfortunately, we are not 
aware of any complete proof for such a result over the last thirty years [U2] . Re- 
cently, important progress has been made in [DcV] and [V] to establish almost 
exponential decay rate for Boltzmann solutions with large amplitude for gen- 
eral collision kernels and general boundary conditions, provided certain a-priori 
strong Sobolev estimates can be verified. Even though these estimates had been 
established for spatially periodic domains [G3-4] near Maxwellians, their validity 
is completely open for the Boltzmann solutions, even local in time, in a bounded 
domain. As a matter of fact, such kind of strong Sobolev estimates may not be 
expected for a general non-convex domain [G4]. This is because even for sim- 
plest kinetic equations with the diflferential operator v-V^, the phase boundary 
d^l X is always characteristic but not uniformly characteristic at the grazing 
set 7q = {{x,v) : x G Sfi, and v ■ n(x) = 0} where n(x) being the outward 
normal at x. Hence it is very challenging and delicate to obtain regularity from 
general theory of hyperbolic PDE. Moreover, in comparison with the half-space 
problems studied, for instance in [LY] , [YZ] , the complication of the geometry 
makes it difficult to employ spatial Fourier transforms in x. There are many 
cycles (bouncing characteristics) interacting with the boundary repeatedly, and 
analysis of such cycles is one of the key mathematical difficulties. 

The purpose of this article is to develop an unified l? — L°° theory in the 

near Maxwellian regime, to establish exponential decay toward a normalized 

i„|2 

Maxwellian fj, = e ^ , for all foTir basic types of the boundary conditions and 
rather general domains. Consec|ucntly, uniqueness among these solutions can 
be obtained. For convex domains, these solutions are shown to be continuous 
away from the singular grazing set 7q. 
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1.1 Domain and Chciracteristics 



Throughout this paper, we define Q = {x : £,{x) < 0} is connected and bounded 
with ^{x) being a smooth function. We assume V^(a:) 7^ at the boundary 
^{x) = 0. The outward normal vector at dCl is given by 

and it can be extended smoothly near dfl = {x : ^{x) = 0}. We say is real 
analytic if ^ is real analytic in x. We define f2 is strictly convex if there exists 
> such that 

%CWCV>cdCI' (4) 

for all X such that ^{x) < 0, and all C € R^. Wc say that has a rotational 
symmetry, if there are vectors xo and w, such that for all x G dQ 

{{x - Xo) X w} ■ n{x) = 0. (5) 

We denote the phase boundary in the space fl x as 7 = 90 x R-*, and 
split it into outgoing boundary 7_,_ , the incoming boundary 7_ , and the singular 
boundary 70 for grazing velocities: 



= {{x,v) e dnxii^ 
7_ = {{x,v) e dnxR^ 
7o = {{x,v) € dQxB? 



n 



{x) • « > 0}, 
n{x) ■ V < 0}, 
n{x) ■ V = 0}. 



Given {t, x, v), let [X(s), F(s)] = [X{s; t, x, v), V{s; t, x, v)] = [x+{s — t)v, v] 
be the trajectory (or the characteristics) for the Boltzmann equation (1): 

with the initial condition: [X{U t, x, v), V{t; t, x, v)] = [x, v]. For any (x, v) such 
that X G tt,v 0, wc define its backward exit time t-h{x,v) > to be the 
the last moment at which the back-time straight line [X{s; 0, x, v), V{s; 0, x, v)] 
remains in Cl : 

tb{x, v) = sup{T >0:x-TveQ,}. (7) 
We therefore have x — t^v € dQ and ^{x — ti,v) = 0. We also define 

Xh{x, v) = x{tb) = X — tbV e dCl. (8) 

We always have v ■ n{xh) < 0. 
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1.2 Boundeiry Condition and Conservation Laws 

In terms of the standard perturbation / such that F = H+ s/Jlf, the Boltzmann 
equation can be rewritten as 

{dt + v-V + L}f = T{f, /), /(O, X, v) = /o(x, v), 

where the standard hnear Boltzmann operator see [G] is given by 

Lf = vf-Kf=-^{Q{,ji,^f) + Q{^f,li)} (9) 

with the colhsion frequency v{v) = j \v ~ ii{u)qo{9)dud9 ^ {1 + \v\}^ for 
< 7 < 1; and 

r(/l, /2) = = rgain(/l, /2) - Tioss (/l , /2) • (10) 

In terms of /, we formulate the boundary conditions as 

(1) The in-flow boundary condition: for {x,v) G 7_ 

f\^_=g{t,x,v) (11) 

(2) The bounce-back boundary condition: for x G dQ, 

f{t,x,v)\^_=f{t,x,-v) (12) 

(3) Specular reflection: for x e dO., let 

R{x)v = v- 2{n{x) ■ v)n{x), (13) 

and 

f{t,x,v)\^_ = f{x,v,v-2{n{x)-v)n{x))= f{x,v,R{x)v) (14) 

(4) Diffusive reflection: assume the natural normalization, 

cj M.)|nW..|* = l. (15) 

J v-n(x)>0 



then for {x,v) € 7_, 



f{t, X, v)\^_ = c^^/ ii{v) 1 f{t, X, v')^/|J,{v'){nx ■ v'}dv'. 

J ■n{x)>0 



(16) 



For both the bounce-back and specular reflection conditions (12) and (14), it 
is well-known that both mass and energy are conserved for (1). Without loss of 
generality, we may always assume that the mass-energy conservation laws hold 
for t > 0, in terms of the perturbation /: 



/ f{t,x,v)^/fldxdv = 0, (17) 
vff{t,x,v)^dxdv = 0. (18) 



JnxR^ 
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Moreover, if the domain Q has any axis of rotation symmetry (5), then we 
further assume the corresponding conservation of angular momentum is vahd 
for alH > : 

/ {{x — xq) X w} ■ vf(t, X, v)y^dxdv = 0. (19) 

For the diffuse reflection (16), the mass conservation (17) is assumed to be 
valid. 



1.3 Main Results 

We introduce the weight function 

u;(w) = (l+pl^H/^e^l"!'. (20) 
where < 61 < i, p > and /3 € 

Theorem 1 Assume ui"^{l + |t;|}^ € in (20). There exists 5 > such that 
if Fo ^ H + i/m/o > 0, and 

IK0II00+ sup e''°'\\wg{t)\\oo < S, 

0<t<oo 

with Ao > 0, then there there exists a unique solution F{t, x, v) = fi+y/Jlf >0 to 
the inflow boundary value problem (11) for the Boltzmann equation (1). There 
exists < A < Ao such that 

sup e^*||M;/(t)||oo < C{||«;/o||oo+ sup e^°'\\wg{s)\U}. 

0<t<oo 0<t<cx) 

Moreover, iffl is strictly convex (4), and if fo{x,v) is continuous except on 70, 
and g{t, x, v) is continuous in [0, 00) x {dfl xR^\ 70} satisfying 

fo{x,v) = g{x,v) on^_, 

then f{t, X, v) is continuous in [0, 00) x {O x \ 70}. 

Theorem 2 Assume iu~'^{l + \v\}^ G in (20). Assume the conservation of 

mass (17) and energy (18) are valid for /q. Then there exists 5 > Q such that if 
Fo{x,v) = n + y/]lfo{x,v) > and ||w/o||oo ^ <5, there exists a unique solution 
F{t,x,v) = iJ. + -s/Jlf{t,x,v) > to the bounce-back boundary value problem (12) 
for the Boltzmann equation (1) such that 

sup e^*|K(t)||oo <C|Ko||oo. 

0<t<oo 

Moreover, iffl is strictly convex (4), and if initially fo{x,v) is continuous except 
on 7o and 

fo{x, v) = fo{x, -v) on dnxR^\ j^, 
then f{t, X, v) is continuous in [0, 00) x {fi x R^ \ 70}. 
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Theorem 3 Assume w~'^ {1 + \v\}^ E in (20). Assume that is both strictly 
convex (4) and analytic, and the mass (17) and energy (18) are conserved for 
fo- In the case of fl has any rotational symmetry (5), we require that the 
corresponding angular mom,entum (19) is conserved for /q. Then there exists 
S > such that if Fo{x,v) = iJ, + y/Jlfo{x,v) > and ||i«/o||oo < <5, there exists a 
unique solution F{t, x,v) = + ,/Jif{t, x,v)>Q to the specular boundary value 
problem (14) for the Boltzmann equation (1) such that 

sup e^*|K(t)||oo <C|Ko||oc. 

0<t<oo 

Moreover, if fo{x,v) is continuous except on^Q and 

fo{x,v) = fo{x,R{x)v) on dQ. 
then f{t, X, v) is continuous in [0, oo) x {f2 x \ 70} . 
Theorem 4 Assume (15). There is Ooii^o) > such that 

6o{i^o) < ^ < ^ ' ^'^'^ P sufficiently small (21) 

for weight function w in (20). Assume the mass conservation (17) is valid for 
fo. If Fq{x,v) = /i + Y^/o(x, u) > and ||t«/o||oo < S .sufficiently small, then 
there exists a unique solution F[t,x,v) = jj + y/Jlf{t,x,v) > to the diffuse 
boundary value problem (16) for the Boltzmann equation (1) such that 

sup e^*|K(0||oo <C|k/o||oo- 

0<t<oo 

Moreover, if ^ is strictly convex, and if fo{x,v) is continuous except on Jq with 



fo{x,v)\.y_ =c^y/iJ / fQ{x,v')^/|^{v'){n{x) ■v'}dv' 

J{nx-v'>0} 

then f{t, X, v) is continuous in [0, 00) x {Q x \ 70}. 
1.4 Velocity Lemma and Analyticity 

In section 2, we first establish some important analytical tools. The first Velocity 
Lemma 5 plays the most important role in the study of continuity and the 
cycles (bouncing generalized trajectories) in the specular case. It implies that 
in a strictly convex domain (4), the singular set 7q can not be reached via the 
trajectory = v, ^ = from interior points inside Q, and hence 7q does 
not really participate or interfere with the interior dynamics. By Lemma 6, 
no singularity would be created from 7q and it is possible to perform calculus 
for the back-time exit time ti,{x,v). This is the foundation for future regularity 
study. Moreover, the Velocity Lemma 5 also provides the lower bound away 
from the singular set 73, which leads to the estimates for repeating bounces in 
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the specular reflection cases. Such a Velocity Lemma 5 was first discovered in 
[G3-4] , in the study of regularity of Vlasov-Poisson (Maxwell) system with flat 
geometry. It then was generalized in [HH] for Vlasov-Poisson system in a ball, 
and it is the starting point for the recent final resolution to the Vlasov-Poisson 
in a general convex domain [HV] with specular boundary condition. 

Lemma 7 gives refined estimates for the operator K-u, with w fjr^f^ . Similar 
estimates were established in [SG]. Lemma 8 states that the zero set of a analytic 
function is of measure zero unless such a analytic function is identically zero. 
This provides a very convenient tool to verify certain geometric conditions of 
general domains for particularly specular reflections. 



1.5 Decay Theory 

Since no spatial Fourier transform is available, we first establish linear ex- 
ponential decay estimates in Section 3 via a functional analytical approach. It 
turns out that it suffices to establish the following finite-time estimate (Propo- 
sition 11) 

\\'Pf{s)\\lds<M^^j^ ||{I-P}/(s)||^ -I- boundary contributions! (22) 

for any solution / to the linear Boltzmann equation 

dtf + vVJ + L,f ^0, f{0,x,v)^ foix,v) (23) 

with all four boundary conditions (11), (12), (14) and (16). Here for any fixed 
{t,x), the standard projection P onto the hydrodynamic part is given by 



P/ = {ait,x) + b{t,x)-v + c{t,x)\v\''}^M^, (24) 
Paf = a{t,x)y^, V^f = b{t,x)v^/JI{v), VJ = c{t,x)\v\^^flI{v), 

and II • III/ is the weighted L"^ norm with the collision frequency ^{v). 

Similar types of estimates like (22), but with strong Sobolev norms, have 
been established in recent years [Gl] via so-called the macroscopic equations for 
the coefficients a, b and c. The key of the analysis was based on the ellipticity 
for b which satisfies the Poisson's equation Ab = — P}/, where d"^ is some 
second order differential operators. In the presence of the boundary condition 
b ■ n{x) = (bounce-back and specular) or 6 = (inflow and diffuse) at 9f2, 
such an ellipticity is very difficult to employ for the weak L^, instead of H^, 
estimate for b in (22). This is due to lack of regularity of b in (22), even the trace 
of b is hard to define. Instead, we employ the hyperbolic (transport) feature 
rather than elliptic feature of the problem to prove (22). By a method of 
contradiction, we can find fk such that if (22) is not valid, then the normalized 
Zk{t,x,v) ^ satisfies \\PZk{s)\\lds ^ 1, and 

l\\{I-P)Z,{s)\\lds<l. (25) 
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Denote a weak limit of Z]. to be we expect that Z = I'Z = 0,hy each of the 
four boundary conditions. The key is to prove that Z^ ^ Z strongly to reach 
a contradiction. By the averaging Lemma [DL2], we know that Zk{s) — > Z 
strongly in the interior oi fl. As expected, the most delicate part is to exclude 
possible concentration near the boundary dfl. Since is a solution to the 
transport equation, it then follows (Proposition 17) that near dfl, on set of 
the non-grazing velocity v ■ n{x) ^ can be reached via a trajectory from the 
interior of f2, which implies that Z^ can be controlled on such non-grazing set 
with no concentration. On the other hand, over the remaining almost grazing 
set V ■ n{x) ~ 0, thanks to the fact (25), we know that 

Zk ~ PZk = {ak{t,x) + bk{t,x) ■ V + Ck{t,x)\v\'^}^yiJ,{v). 

We observe that such special form of velocity distribution PZ). can not have 
concentration on the almost grazing set v ■ n{x) ~ (Lemma 15), and we 
therefore conclude (22). Clearly, the hyperbolic or the transport property is 
crucial to control boundary behaviors via the interior compactness of Zk. 

1.6 L°° Decay Theory 

Section 4 is devoted to the study of linear decay for all four different types 
of boundary conditions: in-flow, bounceback, specular and diffuse (stochastic) 
reflection. In order to control the nonlinear term r(/, /), we need to estimate 
the weighted L°° of wf. We recall that L = v — K, and study the L°° (pointwise) 
decay of the linear Boltzmann equation (23). We denote a weight function 

h{t,x,v) =w{v)f{t,x,v), (26) 

and study the equivalent linear Boltzmann equation: 

{dt + v-Va^ + u- Kyj}h = 0, h{0,x,v) = ho{x,v)=wfo, (27) 

where 

K^h = wK{-), (28) 

w 

together with various boundary conditions (11), (12), (14), or (16). In bounce- 
back, specular, diffuse reflection as well as the inflow case with g = 0,we denote 
the semigroup U{t)ho to be the solution to (27), and the semigroup G{t)hQ to 
be the solution to the simpler transport equation without collision ■ 

{dt + v-V^ + iy}h = 0, h{0,x,v) =ho{x,v) =wfo. (29) 

Notice that neither G{t) nor U{t) is a strongly continuous semigroup in 
[Ul]. 

We first obtain explicit representation of G{t) in the presence of various 
boundary conditions. We then can obtain the explicit exponential decay esti- 
mate for G{t). Moreover, we also establish the continuity for G{t) with a forcing 
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term g if is strictly convex (4) based on the Velocity Lemma 5. To study the 
L°° decay for U {t) , we make use of the Duhamal Principle : 

U{t) = G{t) + f Git- si)K^U{si)dsi. (30) 

"'0 

Following the pioneering work of Vidav [Vi], we iterate (30) back to get: 

U{t) = G{t)+ [ G{t-si)K^G{si)dsi+ [ [ ' G{t-si)K^G{si-s)K^U{s)dsdsi. 
Jo Jo Jo 

(31) 

The idea is to estimate the last double integral in terms of the L"^ norm of 
/ = ^, which decays exponential by decay theory in Section 2. The key 
difficulty lies in the presence of different boundary conditions which could lead 
to complicated bouncing trajectories. Each of the boundary condition presents 
different difficulties. 

Section 4.1 is devoted to the study of inflow boundary condition (11), in 
which the back-time trajectory is cither from the initial plane or from the bound- 
ary. Even though when <? ^ 0, the solution operators for (29) and (27) are not 
semigroups, for any {t,x,v), a similar representation as G{t — si)KwG{si — 
s)KwU{s) is still possible. With the compact property of Kyj (Lemma 7), we 
are led to the main contribution in (31) roughly of the form 

Iff \h{s,X{s;si,X{si;t,x,v),v'),v")\dv'dv"dsdsi. (32) 

Jo Jo ,^)''bounded 

The v' integral is estimated by a change of variable introduced in [Vi], (see also 
in [LY]) 

y = X{s; Si, X{si;t, x, v), v) = x — {t — si)v — (si — s)v' . (33) 

Since det(^) ^ almost always true, the v' and t;"— integration in (32) can be 
bounded by [h = wf) 

I \h{s,y,v")\dydv" <c( [ \f{s,y,v")fdydv")\ 

JQ^v" bounded \JQ.v" bounded / 

For bounce-back, specular or diffuse reflections, the characteristic trajecto- 
ries repeatedly interact with the boundary. Instead of X{s;t,x,v), we should 
use the generalized characteristics, defined as cycles, Xc\{s;t,x,v) in (32) as 
in Definitions 21, 30 and 36. The key question is wether or not the change of 
variable 

y = Xci{s;si,Xci{si;t,x,v),v') (34) 
is valid, i.e., to determine if it is almost always true 

^^(dX.,,..,.X^,.,..,:..)y)^^_ (35) 
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Section 4.2 is devoted to the study of the bounce-back reflection. The 
bounce-back cycles Xci{s; t, x, v) from a given point {t, x, v) is relatively simple, 
just going back and forth between two fixed points Xb{x, v) and Xb{xb{x, v), —v). 
Now the change of variable (34) and (35) can be established by the study of set 
Sx{v) in Section 4.2.2. 

Section 4.3 is devoted to the study of specular reflection. The specular cycles 
Xc\{s; t, x, v) reflect repeatedly with the boundary in general, and 
is very complicated to compute and (35) is extremely difficult to verify, even in 
a convex domain. This is in part due to the fact that there is no apparent way 
to analyze ^^<='(^'^i^c'J^^i'*'^'^)''" ) inductively with finite bounces. To overcome 
such a difficulty, in Section 4.3.2, det can be computed asymptotically in a 
delicate iterative fashion for special cycles almost tangential to the boundary, 
which undergo many small bounces near the boundary. It then follows that 
^^^^ dx,,is■,s^,x^,js^u.x,v),v') ^ ^ q j^^, ^j^^^^ special cyclcs. This crucial observa- 
tion is then combined with analyticity of ^ and Lemma 8 to conclude that the set 
of jg^.| rfXei(s;X.^(«i,x,»,),,;') | = is arbitrarily small (Lemma 34), and the change 
of variable (34) is almost always valid. The analyticity plays an important role 
in our proof. 

Section 4.4 is devoted to the study of diffuse reflection. The diffuse cycles 
Xci(s;i,x,w) contain more and more independent variables and (32) involves 
their integrations. Similar change of variable (33) is expected with respect to 
one of such independent variables. However, the main diflaculty in this case is 
the L°° control of G(t) which satisfies (29). The most natural estimate for 
G{f) is for the weight w = lJi~^, in which the diffuse boundary condition takes 
the form 

h{t,x,v) = Cfj, / h(t,x,v')ij{v'){v' ■ n{x)}dv' 

with Cfj_ f^,^^^^^^ fj,(v'){v' ■ n(x)}dv' = 1. But such a natural (strong) weight 

/x~2 exactly makes the whole linear Boltzmann theory break down (Lemma 7), 
so we have to use a weaker weight, which is closer to /u" 5 . This new weight will 
introduce serious new difficulty since no natural maximum principle is avail- 
able now from the new boundary condition (184). Moreover, for any {t,x,v), 
since there are always particles moving almost tangential to the boundary in 
the bounce-back reflection, it is impossible to reach down the initial plane no 
matter how many cycles the particles take. In other words, there is no explicit 
expression for G{t) in terms of initial data completely. To establish the L°° 
estimate, we make crucial observation in Lemma 37 that the measure of those 
particles can not reach initial plane after fc— bounces is small when k is large. 
We therefore can obtain an approximate representation formula for G{t) by the 
initial datum, with only flnite number of boimccs. 

Section 5 is devoted to the proofs of the main nonlinear decay and continuity 
results of this paper. 

Our contribution opens new lines of research about such interesting ques- 
tions as specular reflections in non-convex domains, decay for the soft potentials. 
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higher regularity in a convex domain, characterization of propagation of singu- 
larity in a non-convex domain, as well as charged particles interacting with 
fields. Moreover, our new — theory will shed new lights in other inves- 
tigations of Boltzmann equation, in which regularity of the solutions is difficult 
to employ [EGM] [GS]. 



2 Preliminary 

Lemma 5 Let ^ he strictly convex defined in (4). Define the functional along 
the trajectories ^^^^^ = V(s), ^^j^ =0 in (6) as: 

a{s)^e{X{s)) + [V{s)-Vax{s))r-2{V{s)-V'ax{s))-V{smXis)). (36) 
Let X{s) e fl for ti < s <t2- Then there exists constant (7j > such that 

< e^«(l^(*^)l+i)*=a(t2); (37) 
g-c,(|y(ti)l+i)tic,(t,) > e-^«(l^(*^)l+i)*=a(t2). 

Proof. Under the convexity assumption (4), we notice that a{s) > for ti < 
s <t2- Since ^^j^ = by (6), we compute the derivative of a{s) in (36) as 

^ = 2C(X(s))[VC(X(s)) • Vis)] + 2[Vis)V'aXis))VismVis) ■ Ve(X(s))] 

-2{Vis) ■ V'^Xis)) ■ t/(.)}[Ve(X(,s)) . V{s)] 
-2[Vis) ■V'aXisWis) ■VisMXis)). 

Note that the second and the third terms on the right hand side cancel exactly. 
By the convexity (4), there exists > so that we can further bound the last 
term by a{s) as: 

\2[Vis)-V'ms))V{s)-V{sM{X{s))\ 

< X \{-2V{s) . v'ax{s)) ■ v{sm{x{s))\ 

< Ci\V{t,)\a{s). 



We therefore have from the definition (36): 

< eiXis)) + [V^{Xis)-V{s)]' + C^\V{t,)Hs) 



ds 

< {C^\V{t,)\ + l}a{s). 
Our lemma thus follows from the standard Gronwall inequality. ■ 

Lemma 6 Let (t, x, v) he connected with {t—th,Xb,v) backward in time through 

a trajectory of (6). 

(1) The backward exit time <b(a;,u) is lower semicontinuous. 
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(2) If 



then {tb{x,v),Xh{x,v)) are smooth functions of {x,v) so that 

njxb) tbn{xb) 

Vxtb = 7 T, Vvtb ~ 



v-n{xb)' " v-n{xb)' 
V^Xb = 1 + VJb O V, VyXb = tbi + Vytb O V. (39) 

Furthermore, if is real analytic, then {tb{x,v),Xb{x,v)) are also real analytic. 

(3) Let Xi G dQ, for i = 1,2, and let (ti,xi,v) a,nd {t2,X2,v) be connected 
with the trajectory ^^^j^ = y{s), ^^^j^ = which lies inside Cl. Then there 
exists a constant > such that 

1 1 \n(xi) • v\ , . 

(4) Define the boundary mapping 

$b : {t ,x,v) — * {t ~ tb, Xb{x, v), v) G R,x{7q U 7^}- (41) 

Then $b and maps zero measure sets to zero-measure sets between either 
{t} X flxK^ and Rx{7(, U 7+} or Rx{j„ U 7+} ^ Rx{7p U 7_}. 

Proof. (1): We need to show that the set {{x,v) : x G fl and tb{x,v) > c} 
is open for all c e R. Let tb{xQ, vq) > c + e, for some £ > small. From the 

definition of t-h{x, v) in (7), a;o ~ svq G f2 for all < s < c + e < t\^{xo, vq)- Since 
fl is open, we deduce that for (x, v) close to {xq, vq), by continuity, x — vs G 0, 
for all c < s < c + e. This implies that tb{x,v) > c. Hence tb{x,v) is lower 
semicontinuous . 

Proof of (2): Since Xb € dn, ^{xbix,v)) = ^{x - tbv) = 0. But from (38) 
and the fact |V^(a;b)| ^ 0, we have 

dtMx - tbv)\t^ -V ■ V^(xb) > 0. 

Therefore, by the Implicit Function Theorem, we can solve for smooth tb{x,v) 
and deduce (39). Furthermore, when ^ is analytic, so are tb and Xb by Theorem 
15.3 in [D]. 

Proof of (3): Notice that for xi G dfl, 

lim \{^^-y}-^i^^)\ = 0. 

Hence, we have \{xi — y} ■ n{xi)\ < C^\xi — for all y € dCl. Taking inner 
product of xi — X2 = v{ti — ^2) with n{xi), we get 

Cskn*! - > C^\xi - X2\^ > \{xi - X2} ■ n{xi)\ = \v{ti - t2) ■ n(xi)|. 
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Wc thus deduce (40) by dividing \ti — t2|. 

Proof of (4): We define a map from RxSilxR'' to {t} x R'^ x R'^ as 

^t:{t',x',v')^it,x' + v'it-t'), v'). (42) 

Recall the boundary map ^]j{t,x,v) in (41). From the definition of tb in 
(7), $b is one to one from either from {t} x fl x R"^ or from Rx{7g U 7_|_} to 
Rx{7o U 7_}. Denote its inverse by $b^- 

In the case that v ■ n{xh) = 0, i.e., ^i,{t,x,v) € 7q, the grazing set, then 
{t,x,v) G 4'f(7o). In (42), 7q is characterized by the five-dimensional space: 
x' e dfl, v' ■ n{x') = 0, t' £ R. Since '^t is a smooth map, ^tilo) i^ ^1^° ''^ 
dimensional space locally at (x', v' , t') G 7o. This implies that ^ ^t(O'o) 

is a zero- measure set in {{} x x R''. 

In the case that v ■ n{xb) ^ 0, we consider the map where £,{x') = 0. We 
may assume that dxi^{x') ^ 0, and x[ = r]{x'2,x'^) with some smooth function 
T]. Now *t = {t,T]{x'2,x'3) + v[{t - t'),X2 + V2{t - t' ) , x'3 + v'^{t - t'),v') and we 
compute the Jacobian matrix of ^'t for {t' ,x'2,x'^,v') to get 

/ -v[ d,>r] d^'v t-a \ 



—v. 

-V' 



1 t-t' 
1 t-t' 



1 
1 
V 1 J 

and its determinant is exactly ±{v'-n{x')}y^l + |V?7p 7^ if (x', f ') <^ 7q. Hence, 
locally, \E'f is a smooth homeomorphism preserving zero-measure sets away from 
7q. Notice that from the uniqueness in part (2) of Lemma 6, = '^t and 
<I>^, = '^^^ locally if v ■ n{xh) ^ 0. Hence, for any A: > 1, by a finite covering for 
a compact set, <E>j, ^preserves the zero-measure sets on 

Ak = {(t',x',v')\ x' G dn,\v' -nix'M > y,\v'\ <k,\t'\ < k}n^b({t} xflxR^). 

k 

To prove preserves the zero-measure sets between {t} x fl x R^ and 

Rx{7q U 7+}, we take any set S € Rx{7q U 7_|_} with \S\ = 0. Clearly, S = 
{Rx7on5}U^^i{Afcn5}. Therefore, %\S) C *f (7o) ur=i *b'(^fe ^15), and 
|*b^(^fc n 5)1 = and |*t(7o)| = 0. On the other hand, if S e {t} x n x R^, 
we have $b(S') = {R x 70 n $b(S')} U^^ {Ak n $b(5')}. If |5| = 0, then 
l^fe n <&b(5')| = 0, because ^^^{Ak fl 4'b('S')} C S has measure zero and $b = 
maps zero- measure sets to zero-measure sets. 
To prove preserves the zcro-measmx; sets from Rx{7qU7^} Rx{7qU 
7_}, we take S* e R x {70U7+} with |5| 0. Consider the set *bM'^'b('S')\7o}- 
For any point {t,x,v) G $^^{$b('S') \ 7o}, we know that v ■ n{xh) 7^ 0. This 
implies from (40) that th > and t — tb < t. We can choose a fixed s between 
t and t — tb. Locally around {t, x, v), 

$b = $b(s)o*«. (43) 
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where defined in (42) which maps Rx{7o U7_|_} to the plane {s} x x R^, 
and $b(s) maps {s} x Q x R"^ to Rx{7g U7_}. Since $b is one to one, we have 
*bH^b(5') \ 7o} C S, so that |$bH^b(S') \ 7o}| = 0. Therefore, from previous 
arguments, ^'^ preserves zero-measure sets from Rx{7q U7^} to {s} x ilx R^, 
while <I>b(s) preserves zero-measure sets from {s} xilx R"^ set to Rx{7g U7_}. 
Hence, by (43), |$b(S)\7ol - |*b(s) o ^',[$-i{$b (5) \ 7o}] I = 0. We therefore 
deduce |<i>b(S')| = 0, for S inside a neighborhood of {t,x,v). 

On the other hand, we take S € Rx{7o U 7_} with |S'| = and we need to 
show |$b^('S')l = ■ For any point {t,x,v) G 'I'b^S') \ 7o, ' n-ix) ^ 0, so that 
if {t',x',v) = ^b{t,x,v), then t' < t. Hence, we again can find t' < s < t such 
that, near {t,x,v), ^j^^ = ° *b^("'*) (43)- Since |$b{$b^('5') \ 7o}l < 

\S\ = 0, for $^1(5) near it,x,v), {<i>^\S) \ 7o} - ^s' ° *b' W[*b{$b'(^) \ 
7o}]. It follows (43) again that \^^\S) \ 7o| = 0, |$b^('S')l = for 5 inside 
a neighborhood of (t — th,Xh,v). The general case follows from a countable 
covering for Rx{7g U 7^} and {s} x O x R^. ■ 

Lemma 7 Recall (9) and the Grad estimate [Gr2] for hard potentials: 

\K{v, v')\ < C{\v - v'\ + \v- ^;'r^}e"^'''"'''''"^ ""'i.^l'i^ ' . (44) 

Recall w in (20). Then there exists < e{9) < 1 and Cg > such that for 
< £ < e(6»), 

[{\v-v'\ + \v- ,'|-i}e-^l— '1^-^^^^^ < (45) 

J w{v') 1 + \v\ 

Proof. By (20), we first notice that for some Cp^^ > 0, 



w{v) 



w{v') 



<C[l + \v-v'\'y''^e 



_,/i2ii/3L-e{i«r-iH'} 



Let V — v' = T] and v' = v — r] in the integral of (45). By (44), we now compute 

'w{v') 



the total exponent in K{v,v') ^^^1 as: 



= -\\v\' + lv-^-\^-^-e{M'-2vv} 

Since 9 < j, the discriminant of the above quadratic form of jr^l and is 

A = (i + 2ef + 2{-e - 1) = 4^2 - ^ < 0. 
^2 4 4 
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Hence, the quadratic form is negative definite. We thus have, for e > suffi- 
ciently small, there is Ce > such that the following perturbed quadratic form 
is still negative definite 

< -C.l^ + lt^-r?]}. (46) 

Therefore, for given \v\ > 1, we make another change of variable rj^^ = i*?' j^} 
and ??_L = 7/ — ??|| so that \v ■ ri\ = \v\\r]^^\ and \v — v'\ > \rij_\. We can absorb 
{1 + |r?|2}l/'l, |r?|{l + by e^l"!', and bound the integral in (45) by (46): 

C0 f (^ + l)e-^H^|re-^«Hxl^iild|,y|||jdr;^ 
Jn? W j_\ L^— CO J 



We thus deduce our lemma since both integrals are finite. ■ 

Lemma 8 Let K{y) be a real analytic function of y G R" in a region such that 
K{y) is not identically zero. Then the set {y : K,{y) = 0} has zero n— dimensional 
Lebesque measure. 

Proof. We use an induction on the dimension n. When n = 1, we assume 
K(y°) = 0. Since k is analytic, for y near we have 

«(?/) = «(y°)+ET^(y-yT- 

fe=i 

Since K{y) is not identically zero, we can always assume a smallest ki such that 
^ ^ 0. We therefore can rewrite 

Hence K{y) = for y — sufficiently small implies y = y^ (an isolated point), 
which has zero one dimensional measure. 

Assume that the lemma is valid for m. For m+1 dimensional case, we assume 
K{y^) = 0. We first notice that by finite open coverings for any compact subset, 
it suffices to show that for any y^ such that k(?;°) = 0, then there is a ball 
{y:\y- y°\ < 6} such that \{y : \y - y°\ < S, K{y) = 0}| = 0. 
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Now for any y i^xp and \y — < 5, since K{y) is real analytic, we have 

where the multi-index k = [/ci, A:2, A;^], A:! = k\\k2\---km^., while {y — y°)^ = 
lii^iiyi — y?)'^'- Since K{y) is not identically zero, there exists k such that 
K^''\y^) 7^ 0. Without loss of generality, we can further assume that ki =/= 0, 
so that {y — 2/°)*^^ contains the factor (yi — y?)'^^. Furthermore, we can assume 
fci > is the smallest among those non-zero terms, so that every term {y — 
contains the common factor (j/i — j/i)*^^. We therefore can rewrite: 



<y) = {yi-y°ip[j2 

For yi ^ y\, K{y) = Q implies 



«.(ri-E^^^(!'-!'°)"--0. (47) 

Clearly, for any given yi, Ki(yi, j/2, y?ra+i) is an analytical function for m 
variables y — [j/2, ym+i]- Therefore, for fixed yi, we can expand Ki(yi,y) 
around [yl, -^yl^+i] to get 

n^yuy)^ E ^LVi|p^(y_yO)^ 

fe=[0,fe2,fe3,...,fem+l] 

Since by our choices of k and k\, the multi-index k — k\ = [0, ^2, ^m+i], and 
we can consider the term 

(fc-fci)! yy V ) 

We compute k^*^"'^^^ from (47) as 

by our choices of k and ki. From the continuity of Ki ^^\yi,y2, ym+i) with 
respect to yi, for |yi — y^] < ^ small, we deduce that 

4'"'^'(2/i,y2",-,y™+i)^o. 

Therefore, Ki(yi,y2, ■■■,ym+\) is an analytical function which is not identically 
zero for all |y — y"| < 5 sufficiently small. By the induction hypothesis, for 
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any fixed j/i, Ki(yi,y2, ■•■,2/m+i) = is a set of m— dimensional Lcbesque zero 
measure set for y2, ys, ...ym+i- We now apply the Fubini theorem to compute 

\{Ki{y) = 0,\y-y°\<6}\= i{i.,i{y)=o}iyi,y2, ■■■,ym+i)i\y-yO\<s 

= y '^{Kiiy)=o}{yi,y2,--,ym+i)l\y-yO\<sdy2dy3...dym+i^ dyi 

< / \{{y2,y3,-,ym+i) ■ Ki{yi,y2,...,ym+i) = 0, \y -y'^l < S}\dyi 

Jn 
= 0. 

Therefore, inside \y — y°\ < ^ for 5 sufficiently small, we have 

{y : K{y) = 0} C {yi = y?} U = 0}, 

and both sets have zero (m + 1)— dimensional Lebesque measure. ■ 
Lemma 9 Recall (20) and (10). We have 

\wT[g^,g2]{v)\ < C{w{v){\v\ + lV\gi{v)\ + \\wgr\\oo}\\wg2\U- 
Proof. First consider the second term T\oss in (10). We have 

/ \u-vy\tJi^'\u)g2{x,u)\du < C{\v\ + ir|H2||oo, 

Hence «;rioss[5'i,5'2] is bounded by 

w\g^\ J \u-vp\n^/\u)g2{x,u)\du<Cw{v){\v\ + ir\gi{v)\x\\wg2\U 

For Fgain in (10), by \u'\'^ + \v'\'^ = \u\'^ + \v\'^, w{v) < Cw{v')w{u'), and 

j qo{e)\u - v\^e-\"-^\"/*w{v)\gi{u')g2{v')\diodu 

< j qo{0)\u - v\'^e-\'^-'"\^/^w{u')w{v')\gi{u')g2{v')\dijdu 

< \\wgi\\^ X ||«;252||oo / \u- vpe-^^-^^'^'^du. 



Since < 7 < 1, this completes the proof. ■ 

3 Decay Theory 

We define the boundary integration for g{x,v), x G dQ, 



/ 9d'y= g{x,v)\nx-v\dSxdv (48) 
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where dSx is the standard surface measure on dil. We also define \\h\\y = 
ll'^lk^. + ll^lk- be the ^^(7) with respect to the measure \nx ■ v\dSxdv. 
For fixed x G dfl, denote the boundary inner product over 7_|_ in v as 

(51,52)7^(^,2:)= / gi{t,x,v)g2{t,x,v)\n{x) -vldv. 

J ±v-n{x)>0 

By (15), we also define a different -projection for any boundary function 
g{x, v) toward the unit vector y/c^^jijv) with respect to (•, •) as: 



P-yg = { 9{t,x,v) \J iJL{v)n{x) ■ vdvjcij, \J n{v). (49) 

J n{x)-v>0 

Our main theorem of this section is 

Theorem 10 Let f(t,x,v) G be the (unique) solution to the linear Boltz- 
mann equation (23) with trace G L^^^ (R_|_; ^^(7)). 

(1) If f satisfies the in- flow boundary condition (11), then there exists A > 
and C > such that for all <t < 00, 

e^^*||/(0ll'<2{||/(0)||2+ f\''^\\g{s)Wi^ ds}. 

Jo 

(2) Let f satisfy the bounce-back boundary condition (12), then there exists 
X>0 andC>0 such that supo<t<oo{e^^*||/WlP} < 2||/(0)||2. 

(3) Let f satisfy the specular reflection condition (14), and the mass and 
energy conservation laws (17) and (18). In the case Q. has any axis of rota- 
tional symmetry (5), we further require that the corresponding conservation of 
the angular momentum (19). Then there exists A > and C > such that 
suPo<t<^{e2^*||/(OlP}<2||/(0)|p. 

(4) If f satisfies the diffusive boundary condition (16) and the mass conser- 
vation (11), then there exists A > andC > such i/iai supo<j<oo{e^^*||/(f)||^} < 
2||/(0)|p. 

We remark that the existence of such a solution / with its trace f^ S 
-^loc (I^+; -^^(7)) (which guarantees the uniqueness) is in general not known for 
the bounce-back and specular boundary conditions within L^ framework. This 
is due to the possible blow-up of L^^^ (R^:L'^{j)) at the grazing set 7q. See 
[BP], [CIP] and [Ul] for more details. On the other hand, /J \\f{s)\\?^ds < 00 
will be established by the study of (27) in Theorems 28, 35 and 41 with property 
h = wf e L°° and h = wf G L°°{'y). We mainly will establish the following: 

Proposition 11 (1) There exists M > such that for any solution f{t,x,v) 
to the linearized Boltzmann equation (23), 

f \\Pf{s)\\lds < M{ f 11(1 - V)f{s)\\lds + f \\f{s%ds}. (50) 
Jo Jo Jo 
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(2) There exists M > such that for any solution f{t,x,v) to the linearized 
Boltzmann equation (23) satisfying the bounce-back boundary condition (12) and 
the mass-energy conservation laws (17) and (18), we have 

[' \\Pf{s)\\lds < M f \\(1-P)f{s)\\lds. (51) 
Jo Jo 

(3) There exists M > such that for any solution f(t,x,v) to the linearized 
Boltzmann equation (23) satisfying the specular reflection condition (14) and 
the mass-energy conservation laws (17) and (18), (in the case SI has any axis 
of rotational symmetry (5), we further assume conservation of the angular mo- 
mentum (19)), estimate (51) is valid. 

(4) There exists M > such that for any solution f{t,x,v) solution to the 
linearized Boltzm,ann equation (23) satisfying the diffusive boundary condition 
(16) and the mass conservation (17), 

f \\Pf{s)\\lds < M{ f 11(1 - P)/(5)||> + \\{I-P,}f{s%^ds}. (52) 
Jo Jo Jo 

We first show that Proposition 11 impHes Theorem 10. 
Proof, of Theorem 10: For any solution / to the Hnear Boltzmann equation 

(23), e^*/(i) satisfies 

{dt + v-V., + L}{e^\f} - \e^\f - 0. (53) 

hetQ< N <t<N + 1, N being an integer. We split [0, t] = [0, N] U[iV, t\. 
For the time interval [N,t\, since fj G L^^^ ; L'^ {^)) , \\f{s)\\'^ds < oo. 
We then establish the energy estimate for [N, t] as 



\\f{t)W'+ l\LfJ)ds+ I I f\s)djds^\\fiN)\f+ I I ns)djds. 
Jn Jn Jj_^_ Jn Jj_ 

(54) 

For the time interval [0, N], (we may assume N > 1), since fj G Lf^^ (R+i -^^(7))) 
/o Il/('*)ll7'^* < We multiply e^*/ with (53) and take energy estimate 
over 0<s<N: 

pN pN 

e''^ll/Wll'+ / e'^%LfJ)ds-X e'^^\\f{s)fds 
Jo Jo 

pN p pN p 

= \\fmf+ / e'^^f{s)d^ds- / e'^^f\s)d^ds. 
Jo J^_ Jo Jj^ 

Dividing the time interval into ^^^Q[k,k + 1) and letting fk{s,x,v) = f{k + 
s,x,v) for A; = 0, 1, 2. ..A'' — 1, we deduce 

Af-l »i 

e'"^||/(iV)|p+^ / {e2Mfc+4(L/fe,/,)-Ae2^{*^+^>||/fc(s)||2}ds (55) 
k=o •'^ 

= \\fm\'+j:ln e^'^'^+^^fl{s)d^ds- f ( e^^{*=+^>/,^(«)d7rf4- 
fc=o l-^o h- Jo Jj+ J 
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Notice that fk{k + s,x,v) satisfies the same Unearized Boltzmann equation (23) 
for < s < 1. 

In-flow boundary condition (11): Multiplying Jqc^'^*' with (50) to each 
fk{s,x,v) and then summing up over k yields 

^ fc=0 -^0 ^0 Jo 

(56) 

Note that j]y = jly + jly ) is the total boundary integration. Since 

(LhJk) > So\\{I - F}fk\\l = - P}fk\\l + y ll{I - P}fk\\l 

and e'^^''' > 1, we apply (56) to the first copy of - P}fk\\l in (55), and 

move the boundary integrals in (56) to the right hand side of (55). Hence, 

e^^^ll/WII^ + ^E / e'''\\Ph\\lds+'-^Y. e^^'=||{I-P}M|> 
fc=o -^0 

< ll/(0)ir + (l + |)E77 -''^'^'^9l{s)d^ds-{l-^-^)Y. f [ e^'^'+^^ fl{s)d^ds. 

Here we have used the fact || • || < Ciy|| • ||^ for hard potentials, and the in-flow 
boundary condition fk = gu on 7_. Combining P/^ with {I — P}/fe, and note 
g2Afe _ g2A{fe+s}g-2As > g2A{fe+s}g-2A^ obtain a positive lower bound in the 
left hand side: 

(niin{|, A}e-^ - C.a) E / ^''^'^^^ 1 1/^^(^)11 ^'^^ > « 

for C^A < min{^, ^je'^^. Changing back to fk{t) = /(i + k) and letting 
1 - ^ > 0, we deduce 

e^^^||/(iV)||2 < ||/(0)||2 + (1 + y) e'^^g^{s)d^ds. (57) 

Notice that e^^* < e'^Ht-N}^2\s s > N, and since t < N + 1, we can 
choose for Sq and A small such that e^^(*~^^(l + ^) < 2. Hence, multiplying 
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e^-^* with (54) and combining with (57) yields 

/' / f{s)d^ds<e^^'\\f{N)\\^+e'^' f / g\s)d^ds 

< 2{| 1/(0)1 r+ /* / e'^^g'{s)djds}. 

Jo Jj_ 

Bounce-back and specular reflections (12) and (14). In both cases, 
the total boundary contribution in (55) vanishes: 

N-l 1 Af-1 1 

J2 / e'^i'^+^yf,{s)djds~J2 / e''^>'+^yf,{s)d^ds = 0. (58) 

fe=0 "'0 •'1- k=0 "'0 

For N < t < N + 1, we use the same procedure as in the in-flow case (55) and 
the positivity (51) to get 

^ ^ fc=o 

For Cj^A = min{^, > , changing back to the original f{k + s) leads 

to 

e2^^||/(iV)||^< 11/(0)11^. (59) 
By (54) and (58), ||/(t)|| < ||/(iV)||. We deduce if e^^^*"^) < 2 

e'^'WfitW < e'^'WfiNW < e'^(*-^) 1 1/(0)1 p < 2||/(0)|p. 

Diff"use boundary reflection (16). We note from (16) and (49) that 
Ij f^{^)d'y = [P^f{s)]'^d'y, so that the boundary contribution in (55) is 

f' f e^^^'^+'^fisfd^ds-f'f e2^f=+«>/2(s)d7ds = - / e^^'[{I-P^}fis)fd'yds. 
Jo Jo J^^ Jo J^^ 

By the same procedure, we obtain from (55) and the positivity (52): 

e^^^||/(^)|P+(min{|,A}e-A_c^A^' |^,2M.+4||/,(,)||2^ 

^ ^ fe=0 

c N-l „i . 

< 11/(0)11^ - (1 - |) E / / e''f^^^y[{I-P,}f{s)fd^ds. 
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For % < 1 and C^A < min{%, ^} e'^^, we have 

e'^^||/(iV)||'< 11/(0)11^. 

Since for [N,t], we have < from (54). We therefore conclude 

the proposition for g^-^^*"^) < 2. ■ 



3.1 Strategy for the Proof of Prop. 11 

The rest of this section is devoted entirely to the proof of the crucial Proposition 
11. The proof of Proposition 11 is based on a contradiction argument. If 
Proposition 11 were false, then there arc no M exists as in Proposition 11 for 
every linear Boltzmann solution. Hence, for any A; > 1, there exists a sequence 
of non-zero solutions fk{t,x,v) to the linearized Boltzmann equation (23) to 
satisfy one of the following: 

(1) In the in-flow case: fk satisfies (11) and 



/' \\PMs)\\lds > k{ f 11(1- P)/,(s)||>+ WMsW^ds). 

Jo Jo Jo 

Equivalently, in terms of normalization Zk(t,x,v) = fk(.t-x.v) ^ have 



:2 j„ 



/ \\PZk{s)\\tds=l, (60) 
Jo 



and 



^'||(I-P)Zfc(.)||2d. + £||Zfc(s)||><i. (61) 
We also have from [dt + v- 'Vx + L,]fk = 0, 

[dt+vVx + L]Zk = 0. (62) 

(2) In the bounce-back case: fk satisfies (12), the mass-energy conser- 
vation laws (17) and (18), and 

/' WPMsWds >kf\\{l- F)Us)\\lds. (63) 
Jo Jo 

Hence, the normalized satisfies (60), (62), and 

j\\{l-V)Zu{s)\\lds<\ (64) 

(3) In the specular reflection case: fk satisfies (14), and the mass and 
energy conservation laws (17), (18). (If the domain Q. has any axis of rotation 
symmetry (5) then fk also satisfies (19)). We note that (63), (60), (62) and 
(64) are all valid for the normalized Z^- 
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(4) In the diffusive reflection case: fk satisfies (16), and the mass 
conservation (17), (62), and 

/' \\PMs)\\lds > k{ f ||(I-P)/fe(s)||>+ f \\{I-P^}fu{s)\\^ ds}. (65) 

JQ Jo Jo 

The normaUzed Zk satisfies (60) and 

11(1 - P)Z,is)\\lds + \\{I - P,}Zfc(s)||2^da < i. (66) 

In all four cases, there exists Z{t, x, v) such that 

.Zfe ^ Z weakly in / || • ||^rfs, 
^0 

since sup;. \\Zk{s)\\1ds<+oo, and from (61), (64), (66) that 

''||(I-P)^fe(5)||>-0. (67) 



/ 

Jo 



'0 

Notice that it is straightforward to verify 
-pZk FZ weakly 



in 

Jo 



Therefore (I - P)Zfe ^ (I - P)Z weakly, and (I - P)Z = from (67) so that 
Z{t, X, v) = {a{t, x)+v- b{t, x) + \vfc{t, x)}y/Jl. (68) 

Note LZk = L{I - P)Zk and ||(I - P)Zk{s)\\lds ^ in four cases. Letting 
A; — > oo in (62), we have, in the sense of distributions, 

dtZ + v VccZ = 0. (69) 

The main strategy is to show, on the one hand, Z has to be zero from (67) 
and one of the inherited boimdary conditions (11), (12), (14) and (16). On 
the other hand, Zj^ will be shown to converge strongly to Z in || • W^ds, and 
/g ||Z|pds 7^ 0. This leads to a contradiction. 

3.2 The Limit Function Z{t,x,v) 

Lemma 12 There exists constants ao,Ci,C2, and constant vectors bo,bi and vj 
such that Z{t, x, v) takes the form: 

^{y kl^ - bo - x + ao} + {-cotx - cix + w X X + bot + bi} ■ V + + cii + C2}|u|^^ 

(70) 

Moreover, these constants are finite: 

\ao\ + \co\ + \ci\ + \c2\ + \bo\ + \bi\ + |tu| < +oo. (71) 
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Proof. Wc first derive deduce (70). Notice that by plugging (68) into (69) and 
comparing coefHcients in front of y/Ji,Vy/Ji^ \ v\'^y/Ji, we deduce the macroscopic 



equations with h = (6^, 6^, h^): 

d^^c = 0, for i = 1,2,3 (72) 

dtc + d^,¥ = 0, fori = 1,2,3 (73) 

d.,¥ + d,,V = 0, fori^^j, (74) 

d^,a + dtb' = 0, fori = 1,2,3 (75) 

dta = 0. (76) 



Since fl is simply connected, from (72), c{t,x) = c{t). Similarly, from (73), 

6^(t,a;) = -dtc{t)xi +b^{t,X2,X3), 
b'^{t,x) = -dtc{t)x2+b'^{t,xi,X3), 
b^{t,x) = -dtc{t)x3+b^{t,xi,X2). 

To determine b^, we first make use of (74) to get 

da;2b'^{t,X2,X3) + dx,b'^{t,Xi,X3) = 

so that d^^b^{t,X2,X3) = 0. Therefore b^ is linear with respect to X2, and 

b\t,X2,X3) = j\t,X3)X2+ g\t,X3). (77) 
Similarly, we also have d^JP{t, xi,X3) = so that 

P{t,xi,X3) = -j\t,X3)xi+g^{t,X3). (78) 

Next, we make use of another equation of (74) to get 

dx3b'^{t,xi,X3) + dx^b^it, xi,X2) = 

with da;JP{t, xi,X3) = -dxJ^{t,X3)xi+dx^g'^{t,X3). From dll}^{t, xi,X2) = 0, 
we have 

P =f{t,xi)x2+g^{t,xi) (79) 

so that —dx3j^(t,X3)xi+dxsg'^{x3)+j^{t,xi)x2 = 0. Taking one more X3 deriva- 
tive, we get 

-dlj'{t,X3)xi+dl^g^{t,X3)^0 

so that d1^j^{t, X3) = d^^g'^{t, X3) = 0. Furthermore, taking two more xi deriva- 
tives, we have d^J^{t,xi) = 0. Hence and g^ can be expressed as 

jHt,X3) = l\t)x3 + h\t), 

f(t,xi) ^ l\t)xi-h\t), 
g\t,X3) = h'^{t)x3+m'^{t). 
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Plugging back into (77), (78) and (79), we deduce 

= {l\t)X3+h\t))X2+g\t,X3), 

P = -(l\t)x3 + h\t))x, + h^{t)x3+m^{t), (80) 
= l^{t)xiX2- h'^{t)x2 + g^{t,xi). 

Finally, from the remaining equation in (74), 

dxsb^{t,xi,X3) + dx^b^{t,xi,X2) = 0. 

By (80), l^X2+dx3g^{t,X3) + l^X2 + dxig^{t,xi) = 0. Hence = 0, g^ is a linear 
function of X3 and g^ is a linear function of Xi : 



Therefore, letting ^(t) = -[h'^{t),h^{t),h\t)] and m(t) = [■m^{t),m'^{t),m^{t)], 
we deduce from a direct computation that 



We also have d^b{t, x) = from (75) and (76). Hence d^c{t) = from (73), and 
c'{t) = Cot + c\ so that 

C= +Clt + C2. 

Hence x a; + m"{t) = and ^"{t) = m"{t) = 0. We can denote 

b = -{cot + ci]x + {w'{Q)t + zu] y. x + hot + 61. 
where zu is a. constant vector. Moreover, from (75), V x dfb = so that 

V X {-cox + w'{0) X a;} = 0. 
This implies that n7'(0) = and from (75) again, 



Lastly, to prove (71), we note that for 1 < i,j < 3, functions 

\xf^/JI, Xi^, tx-v^, x-v^, xxv^, tv^i, v^, l«Pv^ 
arc linearly independent. Therefore, their coefficients co,ci,C2,ao,bo,bi,w are 



g^{t,xi) = h^{t)xi+m^{t) 
g^{t,X3) = -h^it)x3 + m\t). 



b = -c'{t)x + va{t) X x + m{t). 



a = 




- bo ■ X + ao- 




which is finite. ■ 
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3.3 Interior Compactness 

Lemma 13 For any smooth function x{t, x) such that supx CC (0, 1) x fl, then 
up to a subsequence, linife^oo Jq Wxl^k — Z}{s)\\'^ds = 0. 

Proof. We multiply the equation (62) by x to get 

[dt+v- V^]{x^fe} = {[dt + V ■ Vjxl^fe - xLZk. 

Since \ \Zk{s)\\'^ds is uniformly bounded for the hard potentials, by (60) and 

(67), we deduce from the Averaging Lemma [DL], J x{t,x)Zk{v)Xv{v)dv are 
compact in L^([0, 1] x fi) for any smooth cutoff function X-u(^) (see [Gl]). It 
then follows that 

j xZk{v)[l,v,\v\'^]yJJidv 

are compact in ^^([0, 1] x il). Therefore, up to a subsequence, the macroscopic 

parts of Zk satisfy xP-^fe — > X^Z — xZ strongly in L^([0, 1] xfixR^). Therefore, 
in light of 11(1— P)Zfe(s)||^(is — ^ in (67) for all four boundary conditions, the 

remaining microscopic parts x-^fe satisfy limfe^oo llx{I ~ ^}Zk{s)\\'^ds = 0, 
and our lemma follows. ■ 

3.4 No Time Concentration 

We first establish in time estimate for Zk to rule out possible concentration 
in time, near either t = Q or t = 1. 

Lemma 14 &xvpo<t<i,k>i \\Zk{t)\\ < oo. 

Proof. Since ||/fe(s)||^ < oo, \\Zk{s)\\'i^ < oo. Therefore, by the standard 
estimate for (62), we obtain for < f < 1: 

\\Zk{t)\f + j\\Zkis)\f^^ds + 2 j\LZk,Zk){s)ds 

= \\Zkmf+ I \\Zk{s)\\^^ ds. (81) 
Jo 

We first derive an upper bound for Zk{t). In the case of in flow case (11), 
because of (61), (60) and L > 0, we deduce 

\\Zk{t)\\'<\\Zkm\' + l. (82) 

Note ||.Zfe(s)||^^(fs — Jq \\Zk{s)\\'^ ds for either bounce-back or specu- 
lar reflection (12) and (14), hence (82) is clearly valid. In the case of diffuse 
reflection (16), we deduce (82) because 

£ \\Zk{s)\f^_ds = \\P^Zk{s)\f^^ds < j\\Zk{s)\f^^ds. 
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Next, we derive an upper bound for Zk{0). We note that 

C 



< 

k 



j\LZk{t),Zk{t))dt <C f \\{l-V}Zk\\ldt 
Jo Jo 

In the case of the in-flow case (11), by (61) and (81), 

\\Zk{t)\f > \\Zkmf - j\\Zk{s)\f^Js- j\LZk,Zk){s)ds 

> \\Z,mf-^. (83) 

Note that /J | |Zfc(s)| j^^^ds = /q | |^fe(s)| ds for cither bounce-back and 
specular reflection (12) or (14), so that (83) is clearly vahd. In the case of 
diffuse reflection (16), (83) is valid because of (66): 

j\\Z,{s)\\^^_d8 - £ \\Z,{sW^^ds = \\{I - P,}Z,{s)\\^^^ds > -1 

Since \\Zkit)\\^dt < C \\Zkit)\\ldt < C{1 + i} for hard potentials, 
integrating (83) over <t < 1 yields 



wzmw < f\\zum'dt+^ 

Jo 



< cfwzumidt 

Jo 



k 
C 



10 ^ 
by (60) and (67). Our lemma thus follows from (82). ■ 



3.5 No Boundary Concentration 

The most delicate step is to prove that there is no concentration at the boundary 
dCt so that Zk ^ Z strongly in [0, 1] x H x R^. Let 

To this end, we will establish a careful energy estimate in the thin shell-like 
region near the boundary [0, 1] x {il \ fl^i} x R"^. 

Recall n{x) = ^ 0, well-defined and smooth on O \ 0^4 for e small. 

For m > 1/2, for any {x,v), we define the outward moving (inward moving) 
indicator function %+ (X-) as 

X+(.T,'y) = ln\Q^A^)'^{\v\<s-^,n{x)-v>E}i'") 

X-{x,v) = la\Q^4(a;)l{|„|<£-m_„(^).„<_5} (■{;). 
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Our main strategy is to show that the moving (non- grazing) part X±^k 
are controlled by the inner boundary values of Z/- on d^^i = {^{x) = — £^}, 
which are further controlled by the (compact!) interior parts of Zk. Hence, no 
concentration is possible. On the remaining almost grazing part {1 — x±}^ki 
thanks to the fact ||{I — 'P}Zk{s)\\^ds — > 0, no concentration can occur for 
the small velocity set {\v\ > e~"^} U {|n(a;) • v\ < e}. 



Lemma 15 



sup / / , , ^ ,^ \Zk(s,x,v)\^dxdvds < Ce. (84) 



or \v\>e 



Proof. LetPZfe = {ak{t, x)+v-bk{t, x)+\v\'^ Ck{t, x)}^/Jl. Since supk \\Zk{s)\\'^ds is 
finite and [1, v, \v\'^]^ are linearly independent, there is C > (independent of 
k) such that 

r\\ak{sWds+ ['\\bk{sWds+ [' \\ck{s)\\^ds<C r \\Zk{s)\\^ds<a 
Jo Jo Jo Jo 

(85) 

By (67), we can split: 

Jn\n, /i"(-)-l<^ \Zk{s,x,v)fdxdvds 



or Ivl^e 

or \v\'>e ^ or|i;|>e'^ 

< f [ / . ^ ,^ \FZk(s,x,v)\'^dxdvds + 



^ or |i;|>£ 



c 



Even with an extra weight {1 + (/ > 0), the first term can be bounded by 

the Fubini Theorem as 

or li'l^e"'" 
1 f 

2l„ „M , 17,2/ x| , I 2, 



< / / {\ai{s,x)\^\hi{s,x)\^\ci{s,x)\}^ 
Jo Jn\Q^4 

or |u|>£ ^ 

We note that the inner u— integral above is bounded, uniformly in x. In fact, by 
a change of variable t;|| = {n{x) ■ v}n{x), and v± = v — for \n{x) ■ v\ < e, the 
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' \n{x)-v\< 

and 



inner integral is bounded by the sum of 

/ {l+\v\y+^ndv < C r dv\\ [ e-l^^l'/^rfwi < Ce(87) 

J\n{x)-v\<s J-s 

/ {l + \vfy+^ndv < Cs. 

J \v\>6-^ 

Our lemma thus foUows from (85). ■ 

To study the non-grazing parts X±^k, we fix {x,v) G {O \ il^i} x R"^ and 
any moment s such that e < s < 1 — e, and . We define for backward in time 
< t < s : 

X+{t,x,v) = ln\n^i{x - v{t - s})l{iyi<s-^^ri(x-v{t-s})-v>e}{v); (88) 
and for forward in time < s < t : 

X-{t,X,v) = ln\Q^iix - V{t - S})l[\^\<e-m^nix-v{t-s})-v<-e}iv)- (89) 

Both x± solve the transport equations: 

5tx± + ^ • VxX± = 0' x±{s,x,v) =x±{x,v). (90) 
We first prove that 

Lemma 16 (1) ForO < s-e'^ <t < s, ifx+{t,x,v) ^ then n{x) • w > | > 0. 
Moreover, x+i^ — e^, w) = 0, for x G fi \ 17^4 . 

(2) For s < t < s + e'^ < 1, if x~{t,x,v) ^ 0, then n{x) ■ v < < 0. 
Moreover, x_(.s + e'^,x,v) = 0, for x £ n \ 51^4. 

Proof, ft suffices to prove (1), the proof for (2) being exactly the same. First of 
all, by (88), ifx+{t,x,v) ^ 0, then x — v{t — s) G fl\flgi, n{x — v{t — s}) ■ v > s, 
and \v\ < e~™. Hence for |t — s| < e^, for any < ^ < 1, 

\x - 9v{t -s)-{x- v{t - s)}\ < e-'"e^ < s, (91) 

for 2m < 1. Therefore x — 6v{t — s) is also near dfl and \Vn{x — 9v{t — s))\ is 
uniformly bounded. Now, 

n{x) ■ V = n{x — v{t — s}) ■ v — {n{x — v{t — s}) — n{x)} ■ v 
> £- sup \Wn{x-Ov{t- s})\x\t- s\\v\'^. 

0<9<1 

= e[l-Ce-2™+i] > |, (92) 

for 2m < 1. We thus conclude the first assertion. 

To prove the second assertion, let a; G \ ^1^4 so that — ^ ^{x) < 0. If 
X+{s-e^,x,v) > 0, by (88), -s* < ^{x - e^v) < and \v\ < e"™. But 

^{x - e'^v) = ^{x) - £^VC(a;) ■ v + e'^v ■ V^C(^) • ^'^v, 
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for some X is between x and x — e'^v. Since n{x) = ; there exists a constant 

> such that 

-e'V^x) ■ V = -e-'\Wax)\n{x) ■ v < ~\W^{x)\ < 
for a; G il \ ^^^4 . Here we have used the first assertion that n{x)-v> | . Therefore 

^{x-e'^v) < Q-^^+e\-V''^{x)-e''v 



2 

< -9^+C\e\^\ = -^{l-Ce^-^"'}<-e^ 



for 2m < 1, and small e. This is a contradiction to — < ^{x — e^v) < 0. 
Lemma 17 We have the strong convergence 



lim C \\Zk{s)-Z{s)fds = Q 
Jo 



and Jq\\Z{s)\\^ > 0. Moreover, Z defined in (70) satisfies the corresponding 
boundary conditions (11) with g = 0, (12), (14) and (16) 

Proof. By (90), we multiply x± with (62) to get 

[dt+v- V^]{x±Zk} = -x±LZk. (93) 

Since \ \Zk{t)\\'^^dt < oo, applying the estimate backward in time over the 
shell-like region [s — e^, s] x {f2 \ i^^i} x R'^ for outgoing part x+i we obtain: 

II x+Zk{s)\?+f \\X+Zm?,^dt- f ||x+Zfc(t)||2^dt = ||Z,x+(a-£')||' 
+ r \\x+Zkm^^_dt- r \\x+Zk{tW^s_dt- ix+LZk,Zk){t)dt, {M) 

where at the inner boundary ^ = {x : ^(a;) = — e''} x R'^, its normal vector is 
n{x) = 1^1^. We notice that by Lemma 16, = at s — e^, while = on 

7_ and 7!, since n{x\v > for x+(s, x, v) ^ 0. From ||(I— P)Zfc(s)||^rfs < |, 
we get for k large, 

/ {x+LZk,Zk){t)dt < [ f \L{I-P}Zk\\Zk\{t)dxdvdt 

Js-s^ Jo JJlxR3 

< c| /'ii{i-p}z.(i)ie| U'liz.miy 
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Therefore we can simplify (94) as 

C 

71' 



||x+Zfc(s)||2+ / \\x+Zkm'jt< \\x+Zk{t)f^.Jt 



Similarly, we use L"^ estimate forward in time [s, s + x {f2 \ } x for 
incoming part X- to get 



X_Zu{s + e^)f+ \\x.Z,{t)\f^^dt- \\x.Zkm',dt=\\x.Zk{sr' 

J s J s 



17; 

+ \\x_Zk{t%_dt- \\x-Zk{t%^Jt- {x-L{l-'P}Zk,Zk){t)dt. 

J s J s J s 

We notice that X- = at t = s + e^, while x_ = on the incoming part 7_|_ 
and 7^ by part (2) of Lemma 16. Therefore we deduce from (95) 

||x-^fc(s)lP + £^' \\x-Zk{t)\\t,_dt< j^''^ \\x_Zkm'^._dt + -^. 

By Lemma 16, the supports of x±Zk{t) on 75. are contained in \n{x) ■ v\ > ^. 
Combining the ± cases, we are able to estimate x±Zk in terms of the inner 
boundary contributions as 

/ , , s , \Zk(s,x,v)\^dxdv 

rs+e" 



< 



\\X.Zk{tW^s_dt + Jx+Zu{t%._^dt+-^ 

< j ^ l|l{|t>|<e-'", and |n(a:)-!;|>f }-^/s(i)ll7^f^* + (96) 



Since the outward normal at x € 557^4 is n(x), the set {{x,v) : x G dfl^i, \v ■ 
n{x)\ > |} is away from the singular set 7q = {{x,v) : x G d^^i, \v ■ n(x)\ = 
0}. Hence, by (40) in Lemma 6, both the backward or forward trajectories 
emanating from d^l^i x {|^^| < s""*, • 'n{x)\ > |} spend a positive period of 
time inside Cl^i. Since 

{dt + v- v.} {l{|,|<e-n.}(^fc - Z)] = -1{|,|<,-^}L{I - P}Zk, (97) 

we can apply Ukai's trace theorem (Theorem 5.1.1, [Ul]) to l{|^|<£-m}(Z/j — Z) 
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over 0^4 to get 

s+e" 



/ \\'^{\v\<e — ,\vn(x)\>^}{Zk{t) - Z{t)}\\^^sds 

= / \\'^{\vn{x)\>i}{'^{\v\<e-^}{Zk{t) - Z{t))}\\^^eds 

J s—e^ 

< '{l|l{H<^-'"}(^feW-^(i))llQ,.xR3 + l|l{k|<e-'"}mi-P}^fcW}|ln,.xR3} 

< C,j' ' \\Zk{t)- Z{t)\\l^^^^,dt + C, j\\{l-V}Zkmldt 

< cj' '\\Zk{t)-Z{t)\\l^^^^sdt + 



k ■ 

Therefore, for fixed e, we have from the interior compactness in Lemma 13 
m / 

k- 

Hence, for k large, and for any £ < s < 1 — £, by (96) 



hm / ||l{|^|<e-„j^.„(j,)|>£}{Zfe(t) - Z(t)}||Ldi = 0. 



/ / \Zk(s,x,v)\'^dxdv 



"|«-n(x)|>£ 



< 2 



Js-e^ l,).n.r?rM>P Js-£^ \„.n.(T:'\\> F V k 



\v'n{x)\>£ ^ \v'n{x)\>£ 

< £ + 



■ / mm',sds. 



But from Lemma 12, Z{s, x, v) is smooth so its trace is given by (70) as well. 
By (71), since the time interval is small, 

/ \\Z{t)\\'iM < 2e^ X sup \\Z{t)\\'is < C£^ 

Js-e^ 0<t<l 

where C depends on ao, cq, Ci, C2, 6o, 6i and w. We thus deduce that for £ < s < 
1 — £, for k large, 

/ i^.^^^^\Zk{s,x,v)\''dxdv<Ce. (98) 

We are now ready to prove compactness of Zk- We split 

III \Zk{s,x,v) — Z{s,x,v)\^dsdxdv = / + / + / 

Jo J Jn Jo Je Jl-e 
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By Lemma 14, we conclude that the integrals /J + f^_^ are bounded by Ce. On 
the other hand, we further split the main part ^ as 

2 



+2/ / [\Z{s,x,v)\''+ [ [[ \Zk{s,x,v)-Z{s,x,v)\\ 
Js Jn\n^4 Js J Jn^i 

Clearly, the first term is bounded by (98), the second term is bounded by Ce 
thanks to Lemma 15; by (71), the third term is bounded by 

/ [ \Z{t,x,v)\^dxdv <C\Q.\Q.ei\<Ce, 

where C depends on ag, Cq, Ci, C2, 60) ^1 and w. The last term goes to zero as 
A; ^ 00 by Lemma 13. We hence deduce the strong convergence 

III \Zk{s,x,v) — Z{s,x,v)'\^dsdxdv ^ Q 
Jo J Jn 

by first letting s small, then letting fc ^ 00. Prom our normalization ||PZ/j(s)||^ds = 
1 with PZfe = {ah+v-bh + Ch\v\^}^, there exists C > independent of k such 
that ^ ^ 

/ \\PZk{s)W'ds >C [ \\PZk{s)\\lds > C > 0, 
Jo Jo 

because both norms are equivalent to 

/ \ak{s,x)\'^dxds + / \bkis,x)\'^dxds + / \ck{s,x)\'^dxds. 
Jo Jo Jo 

Hence \\Zi.s)\\^ds = limk^^ \\ZkisWds > C> 0. 

Finally, we study the boundary conditions which Z satisfies. In fact, recalling 
(97) and \\Zkit) - Z{t)W'dt ^ 0, we use Ukai's trace theorem to conclude, 
for any fixed s > 0, 

< Clim[J^' ||l|,|<_i^{Zfc(s) - Z{s)}\\'ds + J^' m+v- V.]l|,|<_i,{Z,(s) - Z{s)}\fds] 
^ C-Jnn (^^'||l|„l<^{L{I-P}Zfc(i)}|rdi) =0. (99) 
For the in-flow boundary case, by (61) and the continuity of Z, 

/ \\'^{\vn(x)\>§}i\v\<e-'^ Z {s)\\'^ds = \im / ||l{|„.„(a;)|>|}l|^|<e-mZft(s)||^ds = 0, 
Jo fe— >cx) ^0 



33 



so that Z = on 7. 

For the bounce-back and specular reflections, Zk{t,x,v) = Zk{t,x,—v), or 
Zk{t,x,v) = Zk{t,x,R{x)v). Letting fc — > 00, we deduce that Z satisfies the 
same relation respectively for > |}. Therefore. Z{t,x,v) = Z{t,x,—v) 

or Z{t, X, v) = Z{t, X, R{x)v) respectively by the continuity of Z. 

For the diffusive reflection, notice that on 7_ , 

Zh{t,x,v) = C/^i Zk{t,x,v')y/Jin■v'dv'}^/ fi{v) = ak{t,x)Jc^ii{v) (100) 

Jn-v'>0 * 

Fix e > small and for any x G dCl, on the set {v ■ n{x) > s},Zk — > .Z^ in 
L2([o, 1] X 7). This implies that from (99) 



and 



Notice that for s small, / ^.rt(a;)>e fJ^dv is a finite non-zero constant, indepen- 

and |v|>£~"^ 

dent of a;. It follows that 



/ / |afc(t,a;)|^| dxdt <Ce <oo. 
Jo Jon J 



This implies that P^{Zk}-y^ = ak{t, x)^J Cfj,ii{v) are uniformly bounded in i^([0, 1] x 
7_,_). But from (66), {/ - P'y){Zk}^^ ^ in £^([0,1] x 7+), we deduce that 
{Zk}'y^ are uniformly bounded in i^([0, 1] x 7_|_) with a weak limit. But 
{Zk^-y^ Z strongly in i^([0, 1] x {7^ \ 7q}) by the trace theorem, so that 
{Zk}j^ — » Z weakly in L^([0, 1] x 7_,_) since 70 has zero measure. Hence 

/ Zk{t,x,v')^/]In■v'dv'}^/Jl(v) ^ Z{t, x, v')^n-v' dv'}y/ iJ.{v) 

weakly L'^{[0, 1] x 7_,_). We then recover (16) by letting A; ^ 00 in (100). ■ 



3.6 Boundeiry Condition Leads to Z — 0. 

Since Z now satisfies one of the boundary conditions Z^ = 0, (12), (14), and 

(16), we will show that Z in (70) has to be zero and this leads to a contradiction. 

In the case of in-flow boundary (11), since Z = on 7, from (70), for 

any t and x G dfl, and w e R^, by comparing the coefficients in front of the 

.2 

polynomials of v, we deduce that -|- cif -|- C2} = and 

{—cotx — cix -\- w X X -\- hot -\- 61} = { Y'la;^ — ho-x + ao} = 0. 
Therefore cq = ci = C2 = 0, and 60 = 0- Then ao = and izj x a; -I- 61 = 0, or 
w'^xs — vj^x^ + 61 = —zu^xs + zu^xi +bl = zu^X2 — vj^x\ -h 6f = (101) 
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for all X € dfl. Notice that since £_{x) = is two dimensional, so we may 
assume that {xi,X2) are (locally) independent. Hence zu^ = w"^ = 6f = then 
■cc7^ = 6f = 0, and finally 6} = 0. Therefore we deduce Z = 0. 

In the case of the bounce-back case (12), for any fixed t, because of (97), we 
apply Ukai's trace theorem over [0, i] x x R'^ to get, for any < t < 1, 

Jim ||l{|,|<,-i}{Zfe(i) - Z(i)}|| = 0. (102) 
Therefore, by (17) and (18), 

[ Z{t)y/]1 = lim / Zfe(t)v7^+ lim / Zk{t)y/]2 = 0, (103) 

( Z{t)\v\''^ = lim / \v\''Zk{t)^+ lim / \vfZk{t)^m) 

because the integrations over > j arc bounded by C\\Zk{t)\\ J^_^^^i n^^"^ = 
Ce, from the L°° estimates in Lemma 14. We therefore obtain that for all t, 

J{^\x\'' -bo-x + ao}^+{^+cit + C2}\v\'',/ll = 0, (105) 

{^\x\^-bo-x + ao}\v\^^+{^+c,t + C2}\v\^^ = 0. (106) 



/ 



This implies that cq = Ci = 0. Moreover, since from the bounce-back boundary 
condition Z{t,x,v) = Z{t,x, —v), we must have b{t,x) = in (24), or 

b{t, a;) = X a; + 6ot + 6i = 

for all X e 90. Clearly 6o = as a function of t. Prom the argument after (101), 
ZU2 = = bi. We therefore deduce that from (105) and (106) that 

J ao^dv + C2 j \v\^^dv = j aQ\v\'^^+C2 j \v\^^dv = Q. 

We thus have ao = C2 = 0, then Z = for the bounce-back case. 

The specular reflection is more delicate. Using the same mass and con- 
servation laws (105) and (106), we again have ci = cq = and b{t,x) = 
W2 X X + hot + bi. Now from the specular reflection, we have for any x G dfl, 
b{t, x) • n = or 

{w X x + bQt + 61} • n{x) = 0. 
Hence bo = for all x G dQ and 

{zux x}- n{x) + 61 • n{x) = 0. (107) 

In the case = 0, we have bi ■ n{x) = on dQ. We can choose x' G such 
that 61 1 1 n{x') by taking the minimizer of min|(3.)=o bi ■ x. Hence, bi ■ n{x') = 
and bi = 0. 
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For Tu 7^ 0, let's decompose bi = Pi-^ + fi^f}, where |r/| = 1 and 77 _L ro. 
Then ,7 = {^x ,7} x^. Hence 

^1 = Pi |— r + ^^2^ = /5ir~i +^2{|— 12 X ^} X = P\^\ - a;o X K7. 

|ro| |CC7| 

where xq = — j^y^^ x 77. By plugging this back into (107), we get 

/^ii — r • nix) + w X (a; — xq) • n[x) — 0. 

Once again, wc can choose a point x' G 90 such that w j| n[x'~) (e.g., look for 
minimizer of min^(^)=o ^ ' 2;)- We then deduce w x (x! — xq) • n(a;') = and 
hence /3i = 0. So 

Z = 07 X (a; - a;o) • (108) 

and 07 X (.X — .xo) • n(.T) = for all x G 90. If is not rotational symmetric, 
there is no non-zero w and exist, then we deduce that Z = for the specular 
case. 

On the other hand, if is rotational symmetric, there arc such non-zero 
w and xo for O so that (108) is valid, we then have to use the additional 
conservation law for angular momentum: 

J w X (x — xq) ■ vZ{t)^/Jidv = 

as fc — > 00 of the same expression for Z). (see the proof for (103)). Therefore, we 
combine (108) to get 

J {vj X [x — Xo) ■ v}^iidxdv = 

Therefore m x {x — xq) ■ v = and Z = from (108). 

In the case of diffuse boundary condition (16), because of (17), we have (105) 
and ci = Co = 0. Moreover, we have 

Z{t, X, v) = c^l / Z{t, X, v')y/Jln ■ v'dv'}yfjl. 

Jn-v'>0 

on 7_. Since Vy^, ,JJi, ^fji are linearly independent, this implies for all t and 
X e 90, 

6(i, x) = tn X .T + hi^t + 61 = 0, and C2 = 0. 

Therefore, 60 = 0, and = 61 = as in (101). Therefore, we have from (105) 
ao / ^JJidv = 0. Hence Z = 0. 



36 



4 L°° Decay Theory 

4.1 L°° Decay For In-flow Boundciry Condition 
4.1.1 G(t,0) and Continuity 

As outlined in Section 1.6, we study the L°° (pointwise) decay for the weighted 
h = wf of the linear Boltzmann equation (27) with the in-flow boundary con- 
dition. We first derive explicit formula for solution operator 0) for the 
homogenous transport equation (29) with in-flow boundary condition. Note 
that for non-zero in-flow datum at the boundary, G(t,0) in general is not a 
semigroup. 

Lemma 18 Let ho{x,v) G L°° and wg e . Let {G(t, 0)/io} he the solution 
to the transport equation (29) 

{dt+v-V^ + v}G{t, 0)/io - 0, G(0, 0)/io = /lo, {G{t, Q)ha]^_ = wg. 

For any {x,v), with a; € fi, let ti,{x,v) be its back-time exit time defined in 
Definition 7. Then for a. e. {x,v), 

{G{t,0)ho}{t,x,v) = lt-t^<oe~'''^''^*ho{x-tv,v) 

+lf-tB>oe-"(''^*'' {«;<?} (i - tb, x - ti,v, v). (109) 

Moreover, 

supe^°*||G(f,0)/io|U < \\ho\\oo+supe''°'\\wg{s)\\^. (110) 

f>0 s>0 

Proof. For almost every x,v, along the characteristic line ^ = ^ = 0, 

^{e''i-)-G{T,0)ho} = 0. 

Hence e''^'')'^G(r, 0)/io is constant along the characteristic. Choose any point 
{t,x,v) in [0,oo) X Q X with its backward exit point {t — tb,Xb,v). lit — 
t\i{x, v) < 0, then the backward trajectory first hits on the initial plane {t = 0}. 
On the other hand, ift — tb{x, v) > 0, then the backward trajectory first hits the 
boundary. Since {G(r, 0)/io}7_ = wg a.e., from part (4) of Lemma 6, (109) is 
clearly valid for almost every x,v, x E Ct, and estimate (110) follows immediately 
from (109) with tb=t- {t-tb). ■ 

Lemma 19 Let ft be convex as in (4)- Let ht^ix, v) be continuous in r2xR'^\7Q , 
g be continuous in [0, oo) x {dO. x R'^\7q}, q{t, x, v) be continuous in the interior 
o/ [0, oo) X Qx R^ anc? sup[o oo)xOxR3 l^^^^y^l < oo. Let h{t,x,v) be the solution 
of 

{dt + V -Vx + i^}h = q, h{0) = ho, h^_=wg. 
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Assume the compatibility condition on 7_, 

ho{x, v) = {wg}{0, X, v) (111) 
Then h{t, x, v) is continuous on [0, oo) x {fi x \ 7g}. 

Proof. Let {x,v) ^ 7q and denote its backward exit time [t — th,Xb,v]. Since 
^{e^('')^G(r, s)h} = q along the characteristic ^ = w, § = 0, for i - < 0, 

h{t,x,v)=e-''^''^^ho{x-vt,v)+ [ e-''^'"^^^-'^q{s,x-v{t- s),v)ds. (112) 
\it — t\3 > 0, we then have 

h{t,x,v) = e-''^''^*^{wg]{t-t^„Xb,v) + f e-'''^''^^*-'U{s,x - v{t - s),v)ds. 

Jt-ti, 

(113) 

Since {x, v) ^ 7o, ii x ^ d^, then ^{x) < 0; and if x G dQ, then v\/£,{x) ^ 0. 
This imphes in (36), a{t) > 0. Since (, is convex and ^(a;b) = 0, we now apply 
Velocity Lemma 5 to get 

a{t - tb) = {v ■ VC(xb)}^ > ca{t) > 0. (114) 

We thus conclude v ■ n{xh) 7^ and also th{x,v) > by (40). Therefore, by 

Lemma 6, th{x,v), Xh{x,v) are both smooth functions of {x,v). 

Now take any point {i, x, v) close to {t, x, v) and we separate three cases. If 
t — tb{t,x,v) > 0, when {i,x,v) is close to {t,x,v), i—tb{x,v) > by continuity. 
Therefore 

h{t,x,v) = e-''^'^^*''{wg}ii - tb,xb,v) + f e-''^^^^^-''U{s,x - v{t - s),v)ds. 

(115) 

From the continuity of g away from 7q, the second term above tends to the 
second term in (113). We split the third term into 



1=1+1 + I ' 

Jt-ti, Jt-e Jt-t-b+e Jt-tt, 



where e > is small. The first and the third parts above arc small since 
^ is bounded, from our assumption. Notice that x — v{t — s) is inside the 
interior of for f — fb + £ < s < t — e, the middle term above tends to 
It-tb+e c~'^'"^^''^~''^\{s, x — v{t — s), v)ds in (113), from the interior continuity of 
q. Clearly \h{t,x,v) — hif,x,v)\ ^ as {f,x,v) {t,x,v) in this case. 

In the case t — tb{x, v) < 0, x — vt ^ dSl. Then for {i, x, v) close to {t, x, v), 
we have t — tb{x, v) < 0, x — vt ^ dfl, and 

h{t,x,v)=e-''^''^^ho{x-vt,v)+ [ e-'^''^'^^-''U{s,x-v{t- s),v)ds. (116) 

Jo 
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Since ho is continuous away from 7q, and q is continuous and is bounded in 
the interior, we again deduce that h{t, x, v) h{t, x, v) by the same argument 
as in the first case t — tb > 0. 

Lastly, if t — t\j{x, v) = Q and (112) is valid. By (114), x\i = x — t\yV — x — tv, 
and (x — tv,v) ^ 70- Then for any {i,x,v) near {t,x,v), i — tb could be either 
> or < 0. If t— th < 0, then h{t,x,v) still has the same expression (116) as 
h(t,x,v) and h{i,x,v) h{t,x,tj) as before. On the other hand, if t — > 0, 
h{t,x,v) is given by (115). By the Velocity Lemma 5 and Lemma 6, we have 
that |t — fb| + |Sb — Xbl — > 0, so that by the previous argument, 

lim h{i,x,v) = e-'''^'''>'{wg}{Q,Xb,v)+ I e-'^''^'q{s,x - v{t - s),v)ds. 

But {wg}{0,Xb,v) = ho{xb,v) by the compatibility condition (111), hence this 
limit equals to h{t,x,v) given by (112). ■ 

4.1.2 Decay of In-flow U{t, 0) 

Theorem 20 Let {U{t,0)ho} be the solution to the weighted linear Boltzmann 
equation (27) as 

{dt + v-Va: + i^-K^}U{t,0)ho = 0, U{0,0)hQ = ho, {U{t,0)ho}j_ = wg. 
There exists < A < Aq such that 

supe^*||[/(i,0)/io|U<C{||/io||cx,+ sup e^°'\\wg{s)\\oo}. 

t>Q 0<s<oo 

Proof. By (112) and (113), we have {U{t,0)ho}{t,x,v) = 

lt-tv.<oe-"("^*/io(a: - vt, v) + lt_(^>oe-''(")*''{w;g}(< - t^, x^, v) 
+ f e-''^''^^*~'''^{Ky,U{si,0)ho}{sux - v{t ~ si),v)dsi. 

*/max{0,t— tb} 

Let Xi = X — v{t — .si), t'y^ be the exit time for {xi,v') and x[^ = Xi — v't[^. We 
now further iterate this formula to evaluate {KyjU{si,0)ho} as 

/ K^iv,v'){U{si,0)ho}{si,xi,v')dv' (117) 
= / KUv,v')l,^_t><oe-'"^'''^''ho{x,-v's„v')dv' 

+ / K^{v,v')lo<s,-ti,e-''^'''^''Hw9}{si-t'b,x'b,v')}dv' 

JR3 

+ / ' e-''(''')(«i-") f Ky,{v,v')Ky,{v' ,v"){U{s,0)ho}{s,xi - v'{si - s),v")dv'dv"ds 

We note that ||ii'uj/i||oo < C'H^Hoo from (45) in Lemma 7. Clearly, since 
iy{v), viv') > Vq > for hard potentials, 

g-i/(t-si)g-i'(j)')(si-s) < g-i^oit-s) g-i'(f)(*-«i)g-i'(^j')«b < g-''otgi^o(si-tb)_ 
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Plugging (117) back into {U{t, 0)/io}(i, x, v) yields that all the terms except for 
the last term in (117) are bounded by (0 < A < j/q): 

e-^^'WhoWoo + sup e^-'\\wg{s)\\^ + 

0<s<(X) 

+Ck f {e-'^°*||/io||oo+e-"°* sup e'">'\\wg{s)\\^}dsi 

Jmin{0,t-*b} 0<s<oo 

< C^{i+l}e-''»*{l|/Jo||oo+ sup e'''"\\wgis)\\^}. (118) 

0<s<oo 

Wc now concentrate on the last term in (117) and split the velocity-time 
integration into several regions. We first consider the case \v\ > N. 
CASE 1: For > N. Since from (45) with £ = in Lemma 7, 

Kw{v,v')Kw{v' ,v")dv' dv" < < 

l + \v\ N 

By Lemma 7 again, the double-time integration /Ia,{o,t-tb} /max{o,si-ti.} ^^^^ 
|u| > A'' is controlled by 



N 



[ f \-''''^'-'^\\U{s,0)ho\\oodsdsi < (119) 
JO Jo 



sup{e^||(7(s,0)/io||cx,} / / e"^^^T^dsdsi < — ^— sup{e^||(7(s,0)/io||<x>}, 

s Jo Jo -'^ s 

where we have split the exponent as 

e-.o(*-.)=e-i^e-^^e^, (120) 

and used the fact J^^ e~~^ ~dsdsi < +oo by a direct computation. 

CASE 2: For \v\ < N, \v'\ > 2N, or < 2N, \v"\ > 3N. Notice that we 
have either — u| > N oi \v' — v"\ > N, and either one of the following are 
valid correspondingly: 

\K^iv,v')\ < e-i''"\K^,{v,v')ei\--^'\"\, \K^{v',v")\ < e't^ V»(^^', ^")e*'"'"""'% 

(121) 

From (45) in Lemma 7, both / |is:^(w, wOeil"-"'!"] and / | A:^(w', w")etl"'-''"l' | 
are still finite. By (120), we use (121) to combine the cases of \v' — v\ > N or 
\v' - v"\ > N as: 

rsi 



lniax{0,t-tb} luia.^{0,si-t{^} |^i)|<A'',|u'|>2JV, |<2iV,|u"|>3JV j 

< Ck I r \\U{s,Q)ho\\^\ I \K^{v,v')\dv' + svip f \K^{v',v")\dv"\ 

Jo Jo (^•'|i'|<JV,|i)'|>2JV, v' J\v'\<2N,\v"\>3N J 

< a,ife-t^' / f \-'''>^'-'^\\U{s,0)ho\\oodsdsi 

Jo Jo 

< Ce,Ke-i^\-^sup{e'^'\\Uis,0)ho\\oo}- (122) 

s>0 
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CASE 3: si — s < £, for e > small. We bound the last term in (117) by 

r CKe-''°^'-'^\\U{s,0)ho\\oodsdsi 

*/inin{0,t— tb} si—e 

< CKe=^ f r e=''^{e^\\U{s,Q)ho\U]dsdsx 

Jo Jsi-e 



< Cife^~ sup{e~||L''(s,0)/io||cx)} X / / e ^ dsdsi 

s>0 Jo Jsi — s 



< CifC^ sup{e~||C/(s,0)/io||oo} X e / e = dsi 

s>0 Jo 

< CK£e^sup{e^\\U{s,0)ho\\oo}. (123) 

s>0 

CASE 4. si- s > e, and |^;| < N, \v'\ < 2N, \v"\ < 3N. This is the last 
remaining case because i{\v'\ > 2N, it is included in Case 2; while if \v"\ > 3N, 
either \v'\ < 2N or \v'\ > 2N arc also included in Case 2. We now can bound 
the integral of the third term in (117) by 

C f f r ' e-'°^*-'^K^{v,v')KUv',v"){Uis,0)ho}{s,xi-isi-s)v',v")\ 

imax{0,t-tb} J B J ma,x{0,si-t{^} 

where B = {\v'\ < 2N, \v"\ < 37V}. By (44), K^{v,v') has possible intcgrable 
singularity of , we can choose Km{v, v' , v") smooth with compact support 
such that 

sup / \KNip,v')-K^{p,v')\dv' <^ (124) 

\p\<3N J\v'\<3N -'V 

Splitting 

Ku,{v,v')K^{v',v") = {K^{v,v')-Kn{v,v')}KUv\v") 

+{KUv', «") - Kn{v', v")}Kn{v, v') + Kn{v, v')Kn{v', v"), 

we can use such an approximation (124) to bound the above s\,s integration by 



rt 



sup{e-^||[/(s,0)/io||oo} X <^ sup \K^{v',v")\dv" + sup j \KNiv,v')\dv'}} (125) 

s \\v'\<2N J |j;|<2JV. 



+C [ f e-'"^^'-'^\KN{v,v')KN{v'y)\{U{s,())ho}{8,xr-{s^-s)v',v") 

./ max{0,t— tb} J B ./max{0,si — 



Note that xi — (si — s)v' G Cl for either si — tj, < 0, s > 0, or < si — fj^ < s. 
Split 

rsi — e j-si—e 

/ = / {lsi-t;<0 + lo<si-t;,<s}- 

Jmax{0,si-t;,} Jo 
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for the last main term in (125). Since Kj\f{v,v')Kj\f{v' ,v") is bounded, we first 
integrate over v' to get 

Cn [ {ls,-t'^<r{v') + lo<s^-t'^<siv')}\{Uis,0)ho}is,Xi - (Si - s)v' , v" )\dv' 

J\v'\<2N 

< CnII ln{xi-{si-s)v')\{U{s,0)ho}{s,xi-{si-s)v',v")\^dv'\ 

\^J\v'\<2N J 



< 



N 



/'M2, 



\{U{s,0)ho}{y,v")\ 



1/2 



Here we have made a change of variable y = xi — (si — s)v' € fl, and for 
Si — s > £, ^ > p-. Denote U{s, 0)ho = wf{s) so that / is a solution to the 
linear Boltzmann equation (23) with /(O) = ^ and = g. We then further 
control the last term in (125) by: 

% r \-Mt-s) f U \{U[s,Q)K}{yy)\^dy \ ' dv"dsds, 

^ JmaxfO.t-tb} Jo J\v"\<3N Un ) 



< 



< 



C 



N 



max{0,t— tb} • 

t /^si—e 



c 



N 



JO 

t psi — e 



JO 



-iyo{t—s) 



-i/o(*— s) 



'|u"|<3JV JQ 



I I 

J\v"\<3N Jn 



s. 1/2 

\{U{s,0)ho}{y,v')\^dydv"[ dsdsi 

\ 1/2 

\f{s,y,v")\''dydv'\ dsdsi 



< —^e s 



l\v"\ 



up{e^1l/(s)||} / re-'?(*-«)dsdsi = %e-^*sup{e^1|/(s)||} (126) 

>0 Jo JO ^ s>0 



1/2 



\\fm\ + y^e'''%ie)\f^_de\ 

where we have used crucially part (1) of Theorem 10 with some < A < ^ 
in the last line. Note that since {1 + g ^^(R^), ||/(0)|| = ||w"^/io|| < 

C||/io||cx>, and \\g\\-y_ = \\^g\\-y_ < C\\wg\\^, and e^^^-^o^'^dO < oo, where 
Ao is in Theorem 1. We can then further bound (126) by 



Cn,x ^-xt 



||/io||oo+ sup e^'^'Wwgis) 

0<s<oo 



In summary, replacing j/q, by A and combining (118), (123), (119), (122), 
(125) and (126), we have established, for any £ > and large > 0, 

sup{e^*||C/(t,0)/io||oo}< {£+%}sup{e^^||f/(s,0)/io||oo}+CKSupe'^''^||«;5(s)||„o+Ce,Jv||/io| 

t>0 -'^ s>0 0<s 

First choosing e small, then N sufhcicntly large so that {e + ^} < ^, 

snp{e^'\\U{t,T)h\\^} <2Ck sup e^«'\\wg{s)\\^ + 2Ce,N\\ho\\oo, 

t>0 0<s<(x 

and we conclude our proof. ■ 
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4.2 L°° Decay for the Bounce-Back Reflection 
4.2.1 Bounce-Back Cycles and Continuity of G{t) 

Definition 21 (Bounce-Back Cycles) Let {t,x,v) ^ 'Jq. Let {tQ,XQ,VQ) = 
{t, X, v) and inductively define for k > 1 : 

{tk+l,Xk+l,Vk+l) = {tk - tb{Xk,Vk),Xh{Xk,Vk), -Vk). 

We define the back-time cycles as: 

Xci{s;t,x,v) = Y,Mt,+utk){s){xk+{s-tk)vk}, Vci{s;t,x,v) = Y,Hk+i,t,){s)vk- 

k k 

(127) 

Remark 22 Clearly, we have Vk+i = {—l)'^^^v, for k > 1, 

Xk = \ ' x^ + '^y> X2, (128) 

and let d= ti~ t2, then for k> 1, 

tk - tk+i = fi > i - <b > 0. (129) 

We follow the outline in Section 1.6 and first establish an abstract lemma. 

Lemma 23 Let Ai he an operator on L°°{'y_^_) — > L°°{j_) such that \ \M\\c{l°°,l°°) ■ 
1. Then for any e > 0, there exists h{t) G L°° and h^ € L°° solving 

{dt + v ■V:o + v}h = Q, hj_ = {l-e)Mhj^, h{0,x,v) = ho G L°° . 

Proof. Fix e > 0, we construct a solution by the following iteration (with 
h°^ = 0) for k = 0,1,2.... 

{dt + v-V^ + u}h''+^ =0, /i^+i = (l-£)M/i^^, h''+\0,x,v) = ho. 

We now show h'' and h'^ is a Cauchy sequence. Taking differences, we get 

{dt + v-\7, + z/}{/i'=+i - h''} = 0, /i^+i - = (1 - e)M{h';^ - /i^-i}, 
with zero initial datum — h'^}t=o = 0. Note that from Lemma 18, 

sup||/i^+i(s) - h';^{s)\\^ < (1 - 6)snp\\h';^{s) - h';-'is)\u 

Repeatedly using such inequality for = 1, 2, we obtain 

snpm+\s) - /i^^(s)|U < (1 - ersnp\\h\ is) - h?^ {s)\\^. 

s s 

Hence {h'^^} is Cauchy in L°°(R x 7_), and then both and {h''} are 

Cauchy respectively by Lemma 18. We deduce our lemma by letting fc — > oo. ■ 



43 



Lemma 24 Let ho € i°°(f2 x R^). There exists a unique solution G{t)ho of 
{^t + 1; • Vx + u}{G{t)ho} = 0, {G(0)/io} = ho, (130) 

with the bounce-back reflection {G{t)ho}{t,x,v) = {G{t)ho}{t, x, —v) for x G 
dCl. For almost any {x,v) GflxR^^\ 7o, 

{G{t)ho}{t, x,v) = J2 h.+ut.me-'^''^'ho (Xe,(0), yei(O)) . (131) 

k 

Moreover, e'^''*||G(t)/io||oo < \ \ho\\oo- 

Proof. For any e > 0, by Lemma 23, there exists a solution h^ of 

{dt + v-Va: + i^}h^ = 0, h^{t,x,v) = {1- e)h^{t,x,-v), h^{0,x,v) = ho- 

with finite \\h^{t, ■)\\<x and supj ||/i^(f, OIloo- Such a solution is necessary unique. 
This is because wc can choose + |?;|} G i"'^ so that = — e is a 

solution to the same equation in (130) with the same boundary condition, with 
an additional property Jq ||/^(s)||^rfs < oo. Then uniqueness follows from the 
energy identity for f^. 

Given any point {t,x,v) ^ 7" and its back-time cycle [-^01(5)1 ^1(5)]- We 
notice |14i(s)| = \v\, for all s, and -j^G{s)ho = —uG{s)ho along the back-time 
cycle [Xci{s),Vci{s)] for tk+i < s < tk- Together with the boundary condition 
a.t s = tk, and part (4) of Lemma 6, we deduce that for almost every {x,v), 
e'^('')*G(s)/iQ is constant along its back-time cycle [^ci(s),K;i(s)] in (127). If 
(x, t;) e fi X \ 7o, then th{x, v) > 0, and 

h%t,x,v) = J2 l[t.+„t.)(0)[l - e]'=e-^M*fto (Xe,(0), ^1(0)) , 
fe 

where the summation over k is finite for finite t by (129). For all e, 

e"°*||/i^(t)||oo < ll/iolU, sup \h'{t,x,v)\< sup \h'{t,x,v)\<\\ho\\oo, 

uniformly bounded. We thus can construct the solution h to (130) with the 
original bounce-back boundary condition by taking w — * limit: h{t, x, v) = 
lime__>o h'^{t, X, w), and hj(t, x, v) = lim^^o ^^(ij 2;, v). We thus deduce our lemma 
by letting £ ^ 0. Once again, such a solution h{t, x, v) is necessarily unique in 
the class because = e if„^(L=^(7)). ■ 

Lemma 25 Let ^ be convex as in (4)- Let ho be continuous m O x R^ \ 7^ and 
q{t,x,v) be continuous in the interior of [0,00) xflxH^ (''i^'isup^o,oo)xnxR^ < 
00. Assume the compatibility condition on 7_ : ho{x,v) = ho{x,—v). Then the 
solution h{t, x, v) of 

{dt+v-V^ + iy}h = q, h{0,x,v)=ho, (132) 

with h{t, X, v) = h{t, X, —v), x e dO, is continuous on [0, 00) x {J^ x R^ \ 70}. 
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Proof. Take any point (t, x, v) ^ [0, oo) x 7q and denote its backward exit point 
[ibjajbji'] along the trajectory. Recall its back-time cycle and (24). Assume 
tm+i < < tm- Since ^{e^^^^h} = q along the characteristics, h{t,x,v) takes 
the form 

e-'^^^'/lo {Xm - tmVm, Vm) + (133) 

/ e^''^''''^''^''''q{s,Xk + {s - tk)vk,Vk)ds + / e^"^"^^''^''^ q {s,Xm + {s - tm)Vm,Vm) ds. 
fe=0 "'0 

Since O is convex and {x, v) ^ 7o, then from the Velocity Lemma 5, n(a;i)-t;i ^ 0. 
Notice that Xk € 90 and ^{xk) = for > 1 so that 

a{tk) = {vk ■ V^(a;fe))2. 

We now apply the Velocity Lemma 5 to conclude a{tk) > Ca{t) > and 

Vk • n{xk) 7^ for all fc > 1. By Lemma 6, tk, and Xk for I < k < n are smooth 
functions of {t,x,v). For any other point {i,x,v) which is close to {t,x,v), we 
deduce that tm > 0. 

In the case that tm+i < 0, or equivalcntly. x„i — t^Xm G from continuity, 
im+iit,XiV) < 0. Therefore, h{t,x,v) has the same expression as h{t,x,v) in 
(133). Therefore, h{t,x,v) — > h{t,x,v) because Xk — > Xk, and tk — > tk,Vk — > Ufc 
for 1 < A: < m + 1, as in the proof of Lemma 19. 

On the other hand, if tm+i — 0, or equivalently, Xm+i = Xm — tmXm G 90,, 
and {xm+i,Vm) ^ 7o- Then by continuity, we know that tm+iit,x,v) is close 
to zero. In the case that im+i{t,x,v) < 0, then (133) is again valid and the 
continuity follows as before. However, if im+iit, x, v) > 0, then im+2{i, x,v) <0 
(because f„_|_2 < 0)) ^^id we have a different expression for h{i,x,v) as: 

e~''^'"^^ho {Xm+l - tm+lVm+l,Vm+l) + X] /- e"''^'')^*"*)^ (s, Xk + [s - tk)Vk, Vk) ds 
rim+i 

+ / e-'''^'''>^*-''\{s,Xm+l + {s-im+l)Vm+l,Vm+l)ds. (134) 

Jo 

The last term in (134) goes to zero because tm+i — > 0, and the second term 

in (134) tends to the q integration in (133) since /g" = J*"^^^- We now show 
that the first term (134) tends to the first term in (133) as well. Since Xm+i — 
tm+iVm+i Xm+l Xm+l = Xm — tmVm G 90, the first term in (134) tends 
to 

^0 i^Xm tmXm, Vm) — ^o{Xm tmXm, Vm), 

which is exactly the first term in (133), from the compatibility condition ho{x, v) = 
ho{x, —v) on 7. We therefore conclude the continuity. ■ 

4.2.2 Non-Grazing Condition \Soc\ = 0. 

The following lemma is due to Hongjie Dong: 
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Lemma 26 For any a; e O, define the set 

S^{v) = {v&'R? w- n{xb{x, v)) =0} = {v eB? w- V^{x - tb{x, v)v) = 0}. 

(135) 

If dfl is C^, then \Sx{v)\ = 0, where \ ■ \ is the Lebesaue measure. 

Proof. Wc first note that if w G Sx{v), tlien kv G Sx{v) for all > 0. It 
therefore suffices to show that the surface measure \Sx{v) fl S^| = 0. 

We fix .X G fj and recall Xh = x — thV G dil. li x G 90, we require v-n{x) ^ 
(a zero measure set). Hence = — jf^Efy and x\, ^ x. Our goal is to show that, 
if e Sx{v), then |^ is a critical value of the mapping from dO. — > 



(138) 



\y-x\ 

at y = Xh- Since 0(y) is smooth for y ^ x, hy Sard's theorem, has zero 

measure in and our lemma is valid. 

Indeed, we assume locally around Xb, 50 = (j/i, 2/2, '7(yi, 2/2)) and if u G 
Sx{v), then at (2/1,2/2) = (a;bi,a;b2) : 

V 

= — • n{xb) = {xbi - xijdiri + {xb2 - X2}dir] - ri{xbi, Xb2) + X3 = 0. (137) 
Clearly, since Xb ^ x, from (137), [a;bi — xi, Xb2 — X2] 7^ 0. But [dy^(p, dy^cf)] = 

1 _|_ (v-L-x-Lf+(vi-xi)(r)-X3)dir) . (y-L-xi){y-2~X2 + (ri-X3)d2Vi} 

\v-x\ ~'~ \v-x\^ \v-x\^ 

(V2-X2){yi-Xi + {rt-X3)diri} 1 . (y2-X2)^ + (y2-X2){r]-X3)d2r) 

\v-x\'-^ \v-x\ ly-^l" 

diT) I (rj-X3){vi-xi + {r]-X3)di_r]} 8211 , {rj-X3){y2-X2 + (ri-X3)d2v} 

\y-x\ \y-x\ ly-x^ 

by (137), a direct comupation yields 

{xhi - xi}dy^(l) + {xb2 - X2)dy^(j) = (139) 

at (2/1,2/2) = (a;bi,a;b2)- This implies that is a critical value of 0. ■ 

Lemma 27 Assum,e \v\ < 2N. Then for any e > 0, there exist Se^N > 0, and 
le,N,^ balls B(xi;ri), B(x2, r2)..-, B{xi;ri) C O, as well as open sets Oxi , Ox2,---Ox, 
of the velocity v with \Oxi\ < s for 1 < i < I, such that for any x G Q, there 
exists Xi so that x G B{xi;ri) and for v ^ Ox^, 

\v ■ n{x — th{x,v)v)\ > 5s,N > 0, \v ■ n{x + tb{x,—v)v)\ > 5s^N > 0- (140) 

Proof. Fix e > 0. For any x G O, since \Sx\ = by Lemma 26, there exists an 
open set 0+ such that |0+| < e/2, and \v ■ n{x — ti,{x,v)v)\ ^ 0, for v ^ 0+. 
But from part (2) of Lemma 6, this implies that v ■ n{x — th{x, v)v) is a smooth 
function on the compact set {|w| < 2A^}n{0+}'^. Hence, there exists Sx,e,N > 0, 
such that on the set {|v| < 2A^} n {0+Y, 

\v ■ n{x — tb{x, v)v)\ > Sx,s,N > 0. (141) 
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In particular, for —v ^ 0+ and \v\ < 27V, or equivalently, v ^ —O^ = {—v' : 
v' G 0+}, \v ■ n{x + tb{x, -v)v)\ > d^^e,N > 0. We define = 0+ U {-0+}, 
clearly \0x\ < £• But by part (2) of Lemma 6, for such v ^ O^, both tb{x,v) 
and th{x, —v) are smooth functions of both variables x and v. In other words, 
there exists B{x;rx) such that if y G B{x]rx) and v ^ Ox 

\v ■ n{y T ib(y, ±v)v)\ > Sx,e,N/2 > 0. 

Now for any .x G V., all B{x, Vx) form an open covering for the compact set Cl, 
hence there is a finite /— subcovering B{xi;ri), B{x2,r2)---, B{xi;ri). From our 
construction, for any x G 0,, there exists i, so that x G B{xi,ri) and moreover, 
\v ■ Vn{x =F tb{±v)v)\ > (5xj,£,Ar/2 > 0. We conclude our lemma by choosing 

ds,N = mmj ■ ■ 

4.2.3 L°° Decay of Bounce-back U{t) 

Theorem 28 Assume w~'^{l + \v\} G L^. Let ho = wfo G L°°. There ex- 
its a unique solution f{t,x,v) to the linear Boltzmann equation (23) satisfying 
f{0,x,v) = /o, and h{t,x,v) = U{t)ho to the weighted linear Boltzmann equa- 
tion (27) satisfying h{0,x,v) = ho, both with the bounce-back boundary condi- 
tion. Then there exist A > and C > such that 

e^*||?7(i)/io|U < C-ll/iolloo. (142) 

For the well-posedness for both problems, we know from the Duhamel prin- 
ciple (30), there exists a L°° solution h{t) = U{t)ho to the weighted linear 
Boltzamnn equation (27). By Ukai's trace theorem, it follows that /;,, is also in 
L°^. Therefore, since w-^{l + \v\} G LK f = ^ e and f^ = !^e ^Ll^^W) 
is a solution to the original linear Boltzmann equation (23), which is unique by 
the standard energy estimate. Based on the decay estimate for /, to prove 
the decay estimate, it suffices to establish a finite-time estimate (143). 

Lemma 29 Assume that there exists A > so that the solution f{t, x, v) of (23) 
satisfies e^*||/(t)|| < C||/o||. Let ho = wfo G L°° and h{t) = U{t)ho = wf(t) is 
the solution of (27) where w~'^ G L^. Assume there exist Tq > and Ctq > 
such that the satisfies 

||C/(To)/io||cx> < e-^^«||/io|U +Cto r \\f{s)\\ds. (143) 

Jo 

Then (142) is valid. 

Proof. It suffices to only prove (142) for t > 1. For any m > 1, we apply the 
finite-time estimate (143) repeatedly to functions h{lTo + s) for / = m — 1, m — 
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2,...0: 

\\h{mTo)\\^ < e-^^'>\\h{{m-l}To)\\^+CTo [ " \\fi{m-l}To + s)\\ds 

Jo 

= e-^^°||/i({m-l}To)|U+CTo / \\f{s)\\ds 

J{m-l}To 

/■{m-l}To 

< e-2^^''||M{m-2}To)|U+e-^^°CTo / \\m\\ds 

J{m-2}To 

l-mTo 

+Cto / \\f{s)\\ds 

J{m-l}To 

m-1 /■{m-fc}To 

< e-™^^o||MO)||oo+CTo^e-'=^^'' / 

fe=0 J{m-fe-l}To 

where = U{t)ho. Now by the decay assumption, in the interval {m - 
k-l}To<s< {m- fc}To, we have ||/(s)|| < e-^«||/o|| < e-M'"-'=-i>^''||/o|| 
Hence, ||/i(toTo)||oo is further bounded by 

™-l /-{m-kyTo 

-"^^"ll/iolloo +Cto Y1 / e-^{'"-'=-i>^''||/o||ds 

fc=0 J{m-fc-l}To 



< e-^^^ll/iolloo +CToe^^°mroe-'"^^°||/o|| 

< Gro.Ae ^\\ho\\oo, 

where by u;-^ e L^, ||/o|| = ||w-i/io|| < C||/io||oo and mToe"'"^'^'' < e"^. 
For any t, we can find m such that mTo <t< {m + Ijlo, and 

IIMOIloo < C||MmTo)|U < CTo.Ae-'^ll/iollcx, < {CTo,Ae^}e-^*||/io|U, 
since e 2 <e2'e2.B 

Proof, of Theorem 28: By Lemma 29, wc only need to prove the finite-time 
estimate (143). We use the double Duhamal Principle (31) for semigroup U{t) 
and G{t). We first estimate the first term in (31) by Lemma 24, 

e"°*||G(t)/io||oo < ||/io||oo. (144) 

For the second term in (31), we note that | li^uj/il loo < C'H/iHoo, then by 
Lemma 24, 

/ G{t - si)K^G{si)hodsi <[ e-'">^*-'^^\\K^G{si)hQ\\^dsi = Cte-'">*\\ho\ 
Jo 00 Jo 

(145a) 

Wo now concentrate on the third term in (31) with the double time integral. 
We now fix any point (i, a;, v) so that {x, v) ^ 7q with its bounce-back cycle. 
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Using (131) twice, we obtain 

G{t - si)K^G{si - s)K^h{s) 



where X'^^{s) = Xci{s;sx,Xc\{si),v'), and yj,(.s) = Vc\{s;si,Xc\{si),v'). In the 
case that \v\ > N, we use the same argument in Case 1, (119) in the proof of 
Theorem 20 to conclude: 

C ^ * '^O ^ 

I |G(i- 51)1^^^(51 -5)1^^/1(3)1100^5^51 < — e"— sup{e— ||/i(s)||oo}- 

(146) 

Moreover, since |T4i(si)| = \v\, |V^'i(si)| = \v'\, hence as in Case 2, (122) in 
the proof of Theorem 20, for < N, \v'\ > 2N or \v\ < N, \v'\ < 2N, \v"\ > 3N, 
we deduce for e small, 

I,|<^,|,,|>2^ e-''°^'-^^KUV.i{s,),v')KMis),v")his,X:,iis,), 

or |i;|<Ar,|u'|<2JV>"|>3iV, 

< a,;ve-t^'e-^sup{e^«||Ms)||oo}. 

s 

We need to only consider the case \v\ < N,\v'\ < 2N, \ v"\ < 3N, for which 
we can use the same approximation in Case 4, (125) to obtain an upper bound 

^e-^sup{e^^||Ms)||oo} (148) 



+Cn f r I e-''o(*-'')|/i(s,X^,(5),t;")|. 

Jo Jo J \v\<N,\v'\<2N,\v"\<3N, 



Recall 0^4 = {x : ^(.x) < — £^}. We focus on the second main term in (148) 
and separate two cases: 

CASE 1: Xci(si) and {v' : \v' ■ \Zf(xl\t[])\ I ^ ^'^^^ case, 

since | V^(Xci(si))| ^ for £ small, the second term in (148) is bounded by 



Cn / / e-^°(*-^) / _ dv'\\h{s) 



lO Jo J{v':\v'- \<s,\v'\ + \v"\<5N} 

< Cnsg 2 sup{e 2 ||/i(s)||oo} 11^ ^ dsdsi 

s Jo Jo 

< Cjv£e-^sup{e'^||/i(s)||oo}. (149) 
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CASE 2: \v\ < N,\v'\ < 2N,\v"\ < 3N and either Xei(.si) G a.i or 
^ci(.si) en\n^4 but {v' : \v' ■ > e}. We denote such a set of v, v' 

and v" by A. By using the formula for cycles (128), wc get 

t t' f h—\ 

f e-MC*-) /■ ^ r ip_^^j(,)/,(,,Xe,(si) + X:(i^+i-iD(-l)V + (s-4)(-l)V,^;" 

(150) 

We first claim that the number of bounces are bounded on A : 

k < CTo,N,e. (151) 

Proof of the claim (151): In the first case Xci{s\) e VL^i, we have from the 
mean-value theorem, 

= e(Xel(si) - t'b(^cl(si), V')V') = ^(Xcl(si)) - t'y^v' ■ Ve(S). 

Since \v'\< 2N, and |V$(x)| < C, 

> ie(^ci(gi))i > ^ 
b'-ve(x)i - 

Because i'^. - t'^^^ >t'^,k< ^2|^ and (151) is valid. 

In the case Xci(si) e n \ ^1^4, and {v' : \v' ■ ivf(xl\(ll))\ I > denote 
*bK) = 4(^ci(si), v') and t'b(-t;') = 4(Xci(si), -t;')- We expand 

= ^(Xe,(Si) - 4(t;')^^') = aXolisl)) - 4(^;')VC(Xel(Si)) • v' + i^b(|0}!^,^2^(^^)^.^ 
= ^(Xel(si) + = C(^cl(Sl)) + 4(-«')V^(Xel(si)) • + ^^b(-^0}\ /^2^(^_ 

Since V^$(x±) are bounded, < 2iV, and — e"' < ^(Xci(si)) < 0, for some 
constant Cjv, we have 

-4(^;')VC(Xe,(si))-z;' + Cw{4(^;')}' > 0, 

We thus have f^(z;') > ^^(^g^;^))'^' and4(-v') > -YM^^i)}:^, Since |Ve(Xei(si))- 
> C^e, either t[^{v') > C^e or t'y^{—v') > C^e. But for bounce-back cycles, 

t'k - 4+1 =ti-h = 4("') + 4(-«') > (152) 

for all A; > 1. We therefore have verified the claim (151). 

We are now ready to estimate (150). By Lemma 27, for the given £ > 0, 

there is 5^ ^ > 0, and [_B(.Xi, r^), O^,.] for i = 1, ...,l, \Oxi \ < e. For Xci{si) G O, 
there exists i such that Xci(si) & B{xi, r^) and for v' ^ 

\v' ■ n(Xci(si) T ib(±«')«')l > Ss,N > 0. 
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Hence t\ = si — t\,{x,v'), t'2 = ti — tt,{x,v') — ti,{x, —v') are both also smooth 
with bounded derivatives over {s,v') G [0,ro] x {O^. n \v'\ < 2N}. It thus 
follows from Lemma 6 that t'l and x'l are all smooth functions of v' ^ ■ We 
then split {v' :\v'\< 2N} into 



J A JAn{v'eOa,.} JAn{v'eo^.} 



(153) 



Since J2k It''' ~ lo^ ' ^^^^ P^^^ bounded by 



JO Jv'eO:,,.\v"\<3N I. Jt'^, 



< CNsnp{e'^^\\h{s)\\^} f e-^°' [ ^ 

s Jo Jv'eO^^,\v"\<3N J. Jt 

s Jo Jo 

= Cjvee"T"^ X sup{eT"*||/i(s)||<x>}. 



e " ^ds 



(154) 



By (151), we therefore only need to consider the second part in (153): 



c- 



; e-"'^*-'h[o^,^-^{s)\h{s,X'^i{s),v")\dsdsidv'dv". 

Jv'<^0^.,\v'\<2N,\v"\<3N f. Jt'^_ 



k + l 



We now wish to change variables as 

fc-i 



1=0 



Since for 1 < Z < A:, tj is a smooth function of v' on the set of integration: 
{v' ^ Oxi, \ < 2N}, we can expand the determinant as a cubic function of s : 



det 



dy 



dv' 



[-Ifs'' + q^{v')s'' + q2{v')s + q3{v'), 



where qi{v) are smooth functions of v' . Therefore, by the analytical formula of 
the algebraic cubic equation, there exists up to three (real) continuous functions 
'n^{v') for 1 < j < 3 so that 



(s, v') : det 
For e > 0, we then split 



dy 



dv' 



q \ = \Jj{s:s = i^Av')}. 



10 Jv'iO^,,\v'\<2N,\v"\<3N Jt'^^^ 



X] /, lu,{|s-7,^.|<e} + ln.{|s-r,^-|>£}- 
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The first part with a small s interval {\s — rij\ < s} is bounded by 



< 



^0 Jv'tO^„\v'\<2N,\v"\<3N Jt'^^^ 

Csup{e^^||/i(s)||oo} /* / e--^"* J V l{\s-r,A<efi<s<s,}{s)e'^'ds 

Jo J\v'\<2N,\v"\<3N \^ k,j -^^'k+i 



< CNee-'^*snp{e'^'\\h{s)\\^}. 

s 

On the other hand, for the second main term, we notice that on the compact 
set in < s < Tq and v' G r\j{\s-r]j\ >e}n {v' ^ O^, , \v'\ < 2N}, the function 
J{-§p-} is uniformly continuous, with uniformly bounded derivatives for s,v'. 
There exists a at To > such that 

J{^}>Ce,iV,To>0. (156) 

And for any point {s,v') e [0,To] x {nj{\s - r^j] > e] r\ {v' e O^,, \v'\ < 2N}, 
there exists open set Os.v' such that v' y is one-to-one and invcrtiblc. We 
therefore have a finite covering (depending on e, To,N) Os^.v'^ such that 

nj{|s - r^l > e} n {v' e \v'\ < 2N} C U„0,„,„;^ 

and v' ^ y is invertible on each Os^,v'^. We therefore can change variable 
v' ^ y locally as 



f\-Mt-s)y- f I ^o,s,]{s)his,X'^,{s),v")dv'dv"dsdsi 

< Cn f e-'^'T f"^^'^ f I l^o,s^]{s)e''°'h{s,y,v")—^dydv"dsds^ 

Jo t:;^Jt'^+.(v)Jo,^y^J\v"\<iN J{g^} 

< CN,e,To I e-"""' I ' / e-'^'l [ h^{s,y,v")dy} ' dv" dsdsi 

Jo Jo J\v"\<3N Uq J 



< c. 



N,e,To 



[ e"^"* / ' e''"" I / h^s,y,v")dydvA dsdsi 

Jo Jo yJnx\v"\<3N J 



< CN,e,To f e-'"' r e-'"^\\f{s)\\dsds, (/ = -) 

Jo Jo w 

< CN,e,To^ I " \\f{8)\\d8, 

Jo 

where k < Cto,n,e and m < CTo,N,e- We thus conclude from (155), (154), (149), 
(148), (147), (146), for t < Tq, e^*\\h{t)\\oo is bounded by 

CK(l+t)e-'^*||/io||oo + (^ + C^e)supe'?^||/i(s)||oo+C^,e,To / ' \\f{s)\\ds. 

a<t Jo 
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We first choose Tq so that 2Ck{1 + To)e^—'^" = e"^^", next choose A'' large, 
then e sufSciently small to get + Cn£) < \- We therefore conclude 

sup e^*||/i(t)|U <2CK(l+To)||/io||oo+2C;v,e,To r \\f{s)\\ds. 

0<t<To Jo 

We thus conclude our theorem by letting t = Tq on the left hand side. ■ 

4.3 L°° Decay for SpeculEir Reflection 
4.3.1 Specular Cycles and Continuity of G{t) 

Definition 30 Fix any point {t,x,v) ^ 7o, and define {to,xo,VQ) = {t,x,v), 
and for k> 1 

{tk+i,Xk+i,Vk+i) = {tk - tb{tk,Xk,Vk),Xb{xk,Vk),R{xk+i)vk), (157) 

where R{xk+i)vk = Vk — 2{vk ■ n{xk+i))n{xk+i)- And we define the specular 
back-time cycle 

^ci(s) = Mtk,tk+i)is) {xk + Vk{s - tk)} , V;,(s) = 

fe=i fe=i 

(158) 

Lemma 31 Let be convex (4). Let ho € L°° and G{t)ho solves (130) with 
specular boundary condition h{t, x, v) = h{t, x, R{x)v) for x € dfl. Then for 
almost all {x,v) ^ 7q, 

{G{t)ho}{t, X, v) = e-'^^HoiX^iiO), y,,(0)) 

= Yht>.+ut,){0)e-''^''^'ho{xk-tkVk,Vk). (159) 

k 

And e''«*||G(i)/io||oo < \\ho\\oo- 

Proof. The existence and uniqueness of the solution follows exactly the argu- 
ment in the proof in Lemma 24, with the bounce-back condition replaced by 

the specular reflection. 

If {x,v) G O X R''\7q, then ti,{x,v) > 0. We consider the back-time specular 
cycles of {t,x,v) as [X^^{s),V^i{s)] as in (157). Clearly, |V'ci(.s)| = v. Since 
j^{e~'^'^^^G{s)ho} = for tk+i < s < tk, any k, by part 4 of Lemma 6 and the 
specular boundary condition at tk+i and tk, it follows that e~'^^'"^G{s)ho is a 
constant along the cycle [^^.[(s), V^,(s)]. 

Wc now show for fixed t, the number of bounces k is finite. Since (x, v) € 
^ X \ 7o, by (36), a{t) > 0. By repeatedly applying Velocity Lemma 5 
along the back-time cycle [-''^ci(*)' Ki('5)]> have for all fc > 1 : e~'^*'=a(tfc) > 
e-<^*'=-ia(ife_i) > ... > e-'^*a(i) > 0. But a{tk) = {vk ■ V^{xk)y, we then have 

{vk ■ n{xk)y > Ga{t) > 0, (160) 
for all k > 1, where C depends on t and v. Therefore by (40) in Lemma 6, that 
tk — tk+i > c(tv)\v\^ ^ ^- ^° summation over k is finite. ■ 
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Lemma 32 Let ^ be convex as in (4). Let ho be continuous in Clx R^\7q and 
q(t, X, v) be continuous in the interior o/[0, oo)xrixR'^ and sup[Q oo)xnxR3 I ^^^^y^l < 
oo. Assume that on j_, ho{x,v) = ho{x, R{x)v). Then the specular solution 
h{t,x,v) to (132) is continuous on [0, oo) x {fi x R^\7q}. 

Proof. We only sketch the proof, which is similar to that for Lemma 25. 
Take any point (i, x, v) ^ [0, oo) x 7q and consider its specular back-time cy- 

[^ci('*)' i'^ (159). By repeatedly applying the Velocity Lemma 5 
and Lemma 6, it follows that tk{t, x, v), Xk{t, x, v) and Vkit, x, v) are all smooth 
functions of (t,x,w). We assume that im+i < < t„i, then h{t,x,v) is given 
by (133) with specular cycles \tk,Xk,Vk] € [X^y{s),V^^{s)]. For any other point 
(t, X, v) which is close to [t, x, v). We now show that h{i, x, v) is close to h{t, x, v) 
by separating two different cases. 

In the case that t^+i < 0, or equivalently, Xm — {tm — s)vm G il, away from 
the boundary. By continuity, tm+i < 0. Therefore, we have the same expression 
for h{i, X, v) as h{t, x, v) in (133) with tk, Xk, Vk replaced by tk, Xk, Vk- Therefore, 
since \vi\ = \v\, h{t,x,v) — > h{t,x,v) following from the continuity of ti ti, 

Xi Xi, vi Vi. 

On the other hand, in the case tm-i-i = 0, Xm+i = Xm — tmVm & d^. From 
(160), {xk+i,Vk) ^ 70- Then by continuity, we know that > 0, and im+i is 
close to zero. In the case that tm+i{t, x, v) < 0, then (133) is still valid and the 
continuity follows. However, if im+i > 0, then im+2 < 0, due to tm+2 < tm+i = 
0. Therefore hit, x, v) is given by a different expression (134) with specular cycles 
[tk,Xk,Vk] S [-'^cil*)' notice that since tm+i — > 0, the g— integrals in 
(134) tend to g— integrals in (133) because of /q" = /j • On the other hand, 

since Xm+i - im+iVm+i —>■ Xm+i, Vm+1 —>■ Vm+1 = R{xm)vm, the first term in 
(134) tends to the first term in (133) as 

from ho{x, v) = ho{x, R{x)v) on 7. We thus complete the proof. ■ 
4.3.2 det (1^) Near Sfi 

Assume O is convex as in (4). We now compute det(|^) for a carefully chosen 
specular back-time cycle near the boundary dCl. We assume Xi € dCl. Given £0 
small, we choose vi such that 

|t;i|=£o, vi-n{xi)= i^^^^^^i =^0- (161) 

We shall analyze the specular back-time cycle of (0, a;i,t;i) : {tk,Xk,Vk)- 
Letting Sk = ti,{xk,Vk), we have ^{xi — siVi) = 0, X2 = xi — siVi G dfl and for 
k>2: 

k 

£,{xi-^SjVj) = Q, Vk = R{xk)vk-i, Xk = Xk-i - SkVk G dfl. 
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Proposition 33 For any finite k> 1, 

^ = <5ii + mn\x^)n\xi) + O(£o), (162) 
where O depends on k, and Q is defined as = 0, 

fc-2 fe-2 

C(A;) = 4^(-l)'=-f+i+4^(-l)'=-i-^'C(p)+2 + 3C(A;-l), fork>2. (163) 

In particular, C,{k) is an even integer so that 

det = {Q{k) + 1} + 0(£o) 7^ 0. 

Proof. From n{xj) = |^||^^, since vj = vj-i — 2{n{xj) ■ Vj-i}n{xj), we define 

dj = Vj ■ V^Xj) = -Vj-i ■ V^Xj) > 0. (164) 

By the Velocity Lemma 5 and our ciioice of vi in (161), if X)^=i ■^i < 

liave Cia(O) < Q!(Ep=i Sp) < C'2a(0), for all j = 1,2,..., k- 1. But a(0) = 

{vi ■ V^(xi)}^ ^ Eg (^iJ2p=i = {^j ■ ^^(^j)}^> W6 then have 

Wj • n{xj) = —Vj-i ■ nixj) ^ Ce^. (165) 

We therefore deduce that, by denoting Uj = n{xj), 

\Vj-Vl\ < \Vj - Vj-i\ + \Vj-i - Vj-2\ + - + \V2 - Vi\ 

< 2\vj-i ■ nj\ + 2\vj-2 ■ nj-i \ + ... + 2\vi ■ n2\ 

< 2jCel. (166) 

With the assumption J2j~^ < C, by = eo, we deduce that 

fe-i 

\xk-xi\<C^\vj\<Ckeo. (167) 
i=i 

We first estimate the next Sk- Note that for fc > 2, 

^{xk + Sk-iVk-i) = 0, ^{xk - SkVk) = 0. 
We then use Taylor expansion at Xk to get 

^{xk + Sk-iVk-i) = ^{xk) + Sk-iVk-1 ■ V^(.Xfe) + ^sl_j^Vk-i\7'^S,{xk)vk^i + 0(s|_iWfe_i); 
^{xk - SkVk) = ^{xk) - SkVk ■ V^{xk) + ^slvkV'^^{xk)vk + 0{slvl). 
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But ^|g^ = nk, ^{xk) = 0, \vk\ = \vk+i\ = 0{eo), by (164), (165) and (166) , 
we have 

1-2*^-1 + C>(£o)sfc_i = 0, 

1 VkV'^^{Xk)Vk , ^, . 2 n 

2^'' Jk = °' 



where dk ^ Eq- Therefore, by (166) and (167), 



sfe-i = — ^^^7r~\ — + Q(go)= 2^t^ +Q(£o), 

= , +0{e,)= +0(so), (168) 

Vk^ i{Xk)Vk V-iV t,{Xx)Vx 



so that Sfe+i — Sfc = O(eo), for finite k. 

We now compute from Vk = Vk-i — 2{n{xk) ■ Vk-i}n{xk), 

It = ^^-2(^fe-i-«feRi4-2(t^fc-i-9„;nfe)4-2(^gj^-nfe)4. (169) 

To compute d^i^ri^ in (169), we note n^iy) = ^■^^(y)\ and 

fe-i 

= d,.^{n^{x,-Y,SjVi)} 
i=i 

fc-i fc-i 
= d,n"^{xk)x{-Y,d.{^3^1-T.'0Q:j} 

To compute d^i^sj, we recall ^(xj+i) = $(xj) = so that 
Taking their v[ derivatives, we split i'^to + 

^ao^(a;j+i){-^^Sp-^w°a^!Sp} = ^5oe(a;i+i){^Sj +i;°a„,sj, 

o p— 1 1 p— 1 o 1 

o p=l 1 p=l 
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Subtracting these two identities, we deduce 



dv° 



=1 1 p=i 



o p 

But by the Taylor expansion and (167), 

do^{Xj+i) - doCiXj) = doe^{Xj){x^j+-L - X^j) + 0{\Xj+i - Xjf ) 



Rewriting ^o^o^{xj-^-l)v° = Vj ■ V^(xj-)-i), we therefore have 



Vj ■ V^(a;j+i Sj = - 2^ do^{xj+i) —j sj 



o,e p—1 1 p=l 

p=l ^ p=l 

Since vj ■ Vaxj+i) = -dj+u from (168), Eo,e ^"'^-rf^If'^"^ = "2 + O(£o) and 

o ^ p=l o,e p=l ^ 

p=l ^ p=l 

We first claim that for 1 < j < 



\d.isj\<^, \g^J<C,. (172) 



We shall prove this via an induction of j. In fact, when j = 1, = 6oi, and 
from ^(a;i — sivi) = 0, we deduce 

„ dl^{x2)si _2x /1vo^ 

dv[Si = = O(eo ). (173) 

And a simple induction leads to the desired result (172). 
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Prom (172) and (171), we have 

o ^ p=l 

By letting Aj = X]p=i ^v' moving one copy of ^v' the right 

hand side, we deduce 

so that we can obtain exphcit formula for Aj as 

= i:E(-r'^g».+o,i,. ,174, 

we have used the fact by (173), Ai = 9^; si = ^^^^f^*i- finally we 
recall (170), 3^m^ = O(^), so that the second term on the right hand side in 
(169) is of the order 0{eo)- Hence 

-2(^-nfc)n; + 0(eo). 
Since ^^fe!li'^™(2;fe) = 0(eo)) the second term on the right hand side is 

= ^ '"'fvcfiof' {E^.i«iK - • m)ni + O(.o) by (167) 



^ , g(-l)--^M:^^,^„«^ + 0(.o) by (174). 



^ot® d^^"^'»9^(2;i)'^i = 4 + O(eo) from (164) and (168), we deduce 
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Clearly, ^(1) = 0, and assume (162) is valid up to A; — 1. Then 

dv[ 



_fc = 5;.+^(fc_i)n>',+4^(-l)'=-P-i<(,5„i+C(p)<n'i)ni 



p 

-2{5io + C(fc - iXni) • n>i + 0(eo). 

Notice that njnf = 1, and splitting = + X]j=fe-i ^■iid collecting 

terms, we conclude our proposition. ■ 

4.3.3 Decay for U{t) 

We now fix any point (i, x, v) so that (x, w) ^ 7o- Let the back-time specular cycle 
of [t, X, v) be [Xp,(si), V^i(si)]. By (31), we use twice (159) to derive h{t, x, v) = 

e--M*/io(X,,(0),K.(0)) 

+ J K^{V,,{s^)y)e-'^^^'>-h^{X',MXm)dv' (175) 

+ f r [ e--('')(*-''^)--(''')(«^-«)i^^(yei(si),t;')^«.(K'i(s),^^")/^(^;i(s).^^")- 
where the back-time specular cycle from (si, X^j(si), w') is denoted by 

= Xei(s;si,Xe,(si),z;'), V^{s) = 1/ i(s; si, X,,(si), z;'). (176) 

More explicitly, let tk and tj^, be the corresponding times for both specular 
cycles as in (157). For tk+i < si < tk, t'^z+i < s < t'j., 

X',,{s) = X,,{s;s^,X,,{s^)y)^x'^,+{s-tM' (177) 

where a;'^, = Xc\{tk';s-i_,Xk + {si-tk)vk,v'),v'^, = Vci{tk';si,Xk + {si-tk)vk,v'). 
Recall a in (36) and define naturally 

a{x, v) EE a{t) = f{x) + [v ■ V^{x)f - [vV^^ix)v]^{x). (178) 

We define the main set 

Aa = {{x,v) -.xen, ^ <\v\<N, and a{x,v) > (179) 

We remark that for x is near dfl, det{|^} ^ 3 in Lemma 33 for vi is 
almost tangential to n(x). On the other hand, by (39), it is easy to compute 
for vi = n{x), det{|^} ^ — 1 since tb ^ 0. This implies from continuity that 
there is vi such that det{|^} = even after one specular reflection. However, 
as shown next, such a zero set is small if ft is both analytic and convex. 

Lemma 34 Fix k and k'. Define for tk+i < si < ifc, s S R 

J = Jk,k'{t,x,v,s,,s,v') ^ det (^^^i^^k±^_^kKJ^ . 
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For any e > sujficiently small, there is S{N, e,TQ, k, k') > and an open 
covering U'^iB{ti, Xi, vi; rj) of [0, Tq] x Aa and corresponding open sets Ot^^xi,vi 
for [tk+i + s,tk — s] X R X with \Oti,xi,vi \ < such that 

\Jk,k'{t,x,v,si,s,v')\ >S>0, 

for < t < To, {x, v) e Aa and (si, s, v') in 

Olx„v, n [tk+i+e,tk-e] X [0,To] x {\v'\ < 2N}. 

Proof. Fix {t, X, v) such that x,v ^ Aa in (179), and fix k, k' . Since x,v ^ 7q, by 
Velocity Lemma 5, we deduce that a{tk) = {n{xk)-Vk]'^ 7^ so that tk—tk+i > 
from (40). We note since > from (179), tk — t^+i < A^diamfi. Since for 
tk+i + ^ < Si < tk — ^, Xk — {si — tk)vk S ^, the interior of the domain with 
s sufficiently small. Prom (178) a{xk + (si — tk)vk,v') > for all v'. This 
implies that along its back-time specular cycle [^ci(^)' 'K:i(*)]' '^i'^'i') > ^^'^ 
v'l, ■ n{x'i,) ^ from the Velocity Lemma 5. Clearly, by the Velocity Lemma 5 
and part (2) of Lemma 6, t;, xj, i>; are analytical functions of s\, s, v'. Therefore 
the function Jk.k'{'t,x,v, si, s,v') is well-defined, and analytic for all v' G R^, 
s G R, and tk+\ + § < Si < ife — §• Moreover, expanding as a polynomial of s, 
we obtain 

f dv' ,\ 

J{t,x,v,si,s,v') = det ( j +Pis^ +P2S+P3 

where pi = Pi(t, x, v, si , v') is an analytical function of (si, v') G {tk+i + ^,tk — 
|) X R'^. But at si = tk+i, -'^ci(si) = Xk+i G From Proposition 33, there 

exists v'q with a{tk+i) = {v'q ■ n{xk+i)Y = > such that det {^^^ W^vo 
0. Since vq ■ n{xk+i) ^ 0, by the Velocity lemma 5 and Lemma 6 for such 
v'q, det is continuous with respect to y near Xk+\- In particular, along 

Ki(s)>Ki(s)] in (176), det |„.=„„ ^ for some si at tk+i + f , for 

£ sufficiently small so that Xk + (si — tk)vk ^ Xk + (ife+i — tk)vk ^ Xk+i- 
Therefore, det is an analytical function which is not identically zero, so is 

Jk,k'{t, X, V, si,s, v') as an analytical function of (si, s, v') G {tk+i -I- |, tfe — §) x 
R X R^. By Lemma 8, for each {t, x, v), there exists an open set Ot^x.v of si, s, v' 
in {tk+i + f ,ife - §) X R X R3 such that \Ot,x,v\ < £, and for (si, s,v') ^ Ot^x.v, 
J{t, x, V, Si,s, v') ^ 0. Therefore, by continuity of J{t, x, v, Si,s, v') with respect 
to si,s,v' , there exists 5t,x,vM,Ta,e,ei,k,k' > 0, such that 

\J{t,X,V,Si,S,v')\ > 5t,x,v,N,Ta,e,k,k' > 

for the compact set: 

{s^,s,v') e 0^%,, n [tk+i + f - f ] X [0,To] X {|^;'| < 27V}. 

Since a{x, v)> j^,hy the Velocity Lemma 5 and par (2) of Lemma 6, tk,Xk, 
and Vk are analytic functions respect to {t, x, v), and x^, and t'f,, are analytic with 
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respect to {t, X, v) as well. Therefore, there exists an open ball B{t, x, v; r(^t^x^y^g)) 
such that if {T,y,w) G B{t,x,v;ri^t,x,v,s)), 

tk+i{T,y,w) < tk+iit,x,v) + ^<si<tk{t,x,v)-^<tk{T,y,w), 

E S 

tk+i{t,x,v) + - < tk+i{T,y,'w) +e, tk{T,y,w) - e <tkit,x,v) - -{180) 

Hence by (180), Jk,k'{T,y,z,s\,s,v') is well-defined, continuous, and we may 
then assume 

IT/ /M ^ ^t,x,v,To,N,e,k,k' ^ „ 

\Jk,k'{r,y,z,si,s,v )\ > ^ > 0, 

in B{t, x, v; r(j,,^,,^_,)) x Ol^ ,^ Ci [tfe+i(t, x,v) + ^, tk{t, x, t^) - §] X [0, Tq] x < 
2A''}, and clearly also on the smaller set (by (180)): 

B{t, X, v; r^t,x,v,e) ) x Ot%,„ n [tk+i (r , y, z) +e, tk (r , y, z)-e]x [0, Tq] x { < 2N}. 

Now by a finite covering for the compact set [0,To] x Aa by such B{t,x,v;r), 
there are [ti,xi,vi], ...[tm,Xm,Vm] such that [0,To] x C U^LiBiU, Xi,Vi;ri). 
For any point {t, x, v), there is i so {t, x, v) G B{ti, Xi, Vi; ri{e)) and 

I T/i /\| ^ • ^i,To,N,e.k,k' ^ „ 

\J[t, ,x,v,si^s,v )\ > mm — — — ^>0 

l<i<m 2 

for {si,s,v') e Ol^^^^^^ n [tk+1 + e,tk - e] x [0,To] x < 2N}. m 

Theorem 35 Assume w~'^{l + \v\} G L^. Assume that (, is both strictly convex 
(4) and analytic, and the mass (17) and energy (18) are conserved. In the case 
offl has rotational symmetry (5), we also assume conservation of corresponding 
angular momentum (19). Let ho G L°°. There exits a unique solution to both 
the (23) and (27) with boundary specular condition, and the exponential decay 
(142) is valid. 

Proof. The well-posedness follows from the exact argument in the proof of 
Theorem 28. Thanks to Lemma 29, we only need to establish the finite time 

estimate (143). Recall A^ in (179). 

STEP 1: Estimate of h(t, x, f )1ao • We first express and estimate the main 
part h(t,x,v)lA through (175). As in the case of bounce-back reflection, the 
first and the second terms in (175) are bounded by (144) and (145a) respectively. 

For the third main contribution in (175), notice that along the back-time 
specular cycles Ki(s),Ki(s)] and K,(s), K'i(s)] in (176), \VM\ = 1^1 and 
= l^'l- Hence, the integration over \v'\ > 2N or \v'\ < 2N but \v"\ > 3N 
are bounded by (147). By using the same approximation, we only need to 
concentrate on the bounded set {\v'\ < 2N and \v"\ < 3N} as in (148) of 

/ / '/ e-''°^'-''hAjh{s,X'^i{s),v")\dv'dv"dsids 

Jo Jo J\v'\<2N,\v"\<3N 



I a(X^i(si),v')<e / a(X^i{si),v')>e 
\v'\<2N,\v"\<3N, |i;'|<2JV,|i;"|<3JV, 
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where wc have further spit into a{Xci{si),v') < e and a(Xci{si),v') > e. 

In the case a(Xci(si), u') < e, f{Xci{si)) + [v'- VC(Xci(si))]2 < e. Hence 
for s small, Xci(si) ^ dQ, and |V^(Xci(si))| > 5 . The first part integral is 
bounded by 



Cn [ r e-''°^'-'^\\h{s)\\^dsds^ / ,^ , , 

Jo Jo J «(^ci(si),i; )<£ 

|t;'|<2JV>"|<3JV, 

< C;vsupe-'^^(*-^)||/i(s)|U / , <(:7jv£supe-'?(*-'^)||/i(s)||oo. 

*>« J\v'- \ Zl{xfXi))\ l<2e>'|<2iV>"|<3JV t>s 

Finally, from (177), we bound the first main term a{Xci{si),v') > £ as 



\v'\<2N,\v"\<3N 



fe,fe' |t,'|<2Ar,|'u"|<3Af 

where [fj^, , x^/, v'j^,] is the back-time cycle of {si,Xk + (si — tk)vk,Vk), for i^+i < 

Wc now study x'^,, + (,s — . By the repeatedly using Velocity Lemma 5, 

we deduce for {t,x,v) e Aa and a{X{si),v') > e, \v'\ < 2N : 

a{ti) = {vi ■ > e-{^«^-i>^«a(si) > Cto,c,jv > 0; 

a{t'i) = {v[ ■ n,j}2 > e-{^^^-i>^°a(Xci(si),t;') > Cto,jv,££ > 0. 

Therefore, applying (40) in Lemma 6 yields ti — t^+i > ^^^1'^ and t'l — t'lj^-^ > 
'-^^^ so that 

< = CTo.C,iV, fc' < = CTo,«,iV,e. (181) 

We therefore further split the si— integral as 

Ck,n r f Yl /*" lA„e--°(*-^)|/i(s,4, + (s-4)4,,t;")| 

^t.+i ^|v'|<2JV,|,;"|<3JV fe<CTo,iV,fe'<CTo.JV.c-^*'fe' + i 



r^k—s r^k—s r^fc+i 
Since ^j., J^*,''' = JJ^ , the last two terms make small contribution as 

fc' + 1 

eCK,N sup e-^«(*-«)||/i(s)|U / = eCK,N sup e-^°(*-«)||/i(i 

0<s<t Jo i|i;'l<2JV.Ii;"l<3JV 0<s<t 
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For the main contribution By Lemma 34, on the set O'l. .^. ,,. fl [tfe+i + 

e, ffc — e] X [0, To] x {\v'\ < N}, we can define a change of variable 

y = x'k, + {s - t'k,)v'k, , 

so that det(^) > S on the same set. By the Implicit Function Theorem, there 
are an finite open covering Uj^^Vj of O^.,^.,^. D [tk+i +e,tk—e]x [0, Tq] x {\v'\ < 
N} , and smooth function Fj such that v' = Fj{t,x,v,y,si,s) inVj. We therefore 
have 



k,k' Jtk+i+e J\v'\<2N, Jt'^,^, k,k' J^k+i+e J\v'\<2N, Jt'^,^, j^k,k' •^k'l<2iV, 



Since J2k' It'''' ~ lo^ — Io° ^^'^ \Oti,xi,vi\ < £, the first part is bounded by 

fc' + 1 

CM£e~~^* sup^{e^^||/i(s)||oo} from Lemma 34. 

For the second part, we can make a change of variable v' ^ y = x'^ + {s — 
t'i.i)v'y on each Vj to get 

Ce,n,NY. I I e-''('')(*-«)|ft(s,4,+(s-4K',OI 

= C,,To,N^ [ I e-(")(*--) |fe (5, y, v") I ] dydv"dsds, 
yJv,J\v"\<m |det{g^V}| 

'e,T„N /■* r\-.ot r g.o.j f h''{s,y,v")dyy^\v"dsdsi 

Jo Jo J\v"\<3N [Ja J 



5 

^£,To,N 



< a 



[ \\f{s)\\ds, 
Jo 



where / = ^- We therefore conclude, summing over k and k' and collecting 
terms 

\\h{t,x,v)lAj\oo < {l + CKt}e-"°*||/io||oo + {^+Cjv,To£}supe-'?{*-«>||/i(s)||oo 

+C,,N,To f \\f{s)\\ds. (182) 
Jo 

STEP 2: Estimate of h{t, x, v). We now further plug (182) back in: h{t, x, v) = 
G{t,s)ho + J* G{t, si){Kyjh{si)}dsi to get 

mmoo < e-'^'WhoW^ + f e-^o^'-'^^WK^hWMds. (183) 

Jo 

But {Kyjh}{si,x,v) = J Kyj{v,v')h{si,x,v')dv' and we split it as 

j Kyj{v,v')h{si,x,v'){l- lA„(x,v')}dv' + j Kyj{v,v')h{si,x,v')lA„(x,v')dv'. 
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By the definition of Aa in (179), the first term is bounded by 

(/ \K^{v,v')\dv' + ( \K^{v,v')\\\\h{si)\\^. 

\J\v'\>N,ot\v'\<^ Ja{x,v')<^ J 

By (45) and (124) if necessary, ij^,|>^_ \v'\<^ \K^{v,v')\dv' = o(l) as N ^ 
oo. From a{x,v') < j^, ^^(a;) + [v'- V^(a;)]2 < i. For AT large, x ^ dfl and 
|V^(a;)| > 5 so that 

/ \K^{v,v')\dv' < f \K^{v,v')\dv' = o{l) 

as N ^ oo. We apply (182) to the second term to bound ||{-K'u,/i}(si)||(x> as 
{l+CxSi}e-^«^^ \\ho\\oo+{o{l)+^+CN,Toe} sup e-"^^''-'^ \\h{s) \ \^+C,,n,To H 1 1/(^)1 M^- 

s Jo 

Hence, by (183), ||/i(t)||oo is bounded by 

e-"°*||/io||oo + t e-'°'{l + CKSi}\\ho\\oo + 
Jo 

+ f e-o{*-i}{o(l) + £+C;v,To£}supe-^{«-«>||Ms)||cx> + Ce,JV,To H \\f{s)\\ds}dsi 
Jo -'^ s<t Jo 

< {1 + CKt^je-^^'WhoWoo + C{o{l) + ^+ Cjv,To£}sup{e-T{*-«}||/i(s)||^} + C,,n,To f \\f{ 

s<t Jo 

We choose Tq large such that 2{1 + CKT^}e-^'^° = e"^"^", for some A > 
0. We then further choose N large, and then s sufficiently small such that 
C{o(l) + jf+ Cn,To£} < 5- We there have 

sup {e^^||/i(s)|U} < 2{1 + Cxi'}||/io||oo + Cto f \\f{s)\\ds. 

0<s<t Jo 

Choosing s = t = Tq, we deduce the finite-time estimate (143), and our theorem 
follows from Lemma 29. ■ 

4.4 L°° Decay of Diffuse Reflection 

4.4.1 Infinite Cycles and Bound for Diffuse G{t) 

In this section, we study the L°° decay of the diffuse reflection. Define h = fw 

to satisfy 



{dt + v-Va; + i^}h = 0, h(t,x,v)L_ = [ h(t,x,v')w(v')da[184:) 

W{V) Jv(x) 

1 



(^) JVix) 

where V{x) = {v' eR^ :v' ■ n{x) > 0}, w{v) = ^ ,(185) 



64 



and by (15), the probability measure da = da{x) is given by 

da{x) = c^n{v'){n{x) ■ v'}dv'. (186) 
For ^ — ^ > and 6 > 6o, and for p sufficiently small, 

„ f 2 11 

* = /! * i2Nfl > 1' / '^'^da - / e"^^l''l'{n(a;) • v}dv = 2 < 



(l+p|z;|2)/3 '- ^' J^"" ^'^ ' Jy'^ L-n-/ -J- 16^2 -^g^2- 

(187) 

We have used the normalization (15) and a change of variable v = ^J^v' to 
evaluate the integral. 

Definition 36 Fix any point {t, x, v) ^ 7o, and let {to, xq, vq) = {t, x, v). Define 

the back-time cycle as 

{tk+i,Xk+i,Vk+i) = {tk-tb{xk,Vk),Xb{xk,Vk),Vk+i), forvk+i e V^+i = {vk+i-n{xk+i) > 

(188) 

And 

Xci{s;t,x,v) = ^lit^:+,,tk){s){xk+{s-tk)vk}, Vci{s;t,x,v) = XI ^ [*/=+!.*<= 

k k 

We define the iterated integral for k>2 

I nfr/da, = I ... J / duk-i \ dai (189) 

We note that each vi {I = 1, 2, ...) are all independent variables, however, the 
phase space Vi implicitly depends on {t,x,v,vi,V2, ■■■Vi-i). We first show that 
the set in the phase space lifZiVi not reaching t = after k bounces is small 
when k is largo. 

Lemma 37 For any £ > 0, there exists ko{e,To) such that for k > ko, for all 
{t,x,v),0 <t <To, X gD, and v e R^, 

/ ^ , '^{tk{t,x,v,vi,V2...,Vk-i)>0}^'l=l'^'^l - ^■ 

Proof. Choosing < (5 < 1 sufficiently small, we further define non-grazing 
sets for 1 < Z < A; — 1 as 

Vf = {vi GVi: VI- n{xi) > 6} D {vi G V; : |^;;| < 

Clearly, by the same argument in (87), 

/ dai< I da I + [ dai < Cd, (190) 

Jv,\Vf Jvrn{x,)<S J\'>'i\>7 
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where C is independent of I. On the other hand, if vi € Vf , then from diffusive 
back-time cycle (188), we have xi — xi+i = {U — ti+i)vi. By (40) in Lemma 6, 
since < y, and vi ■ n{xi) > S, 



{ti - ti+i) > 



Therefore, if ife(t, x, v, vi, v^..., Wfe-i) > 0, then there can be at most 
number of Vi G V/ for 1 < / < A; — 1. We therefore have 

/ ■■■{ 'i-{tk>o}d'^k-i} dak-2---d(7i 

JVi [JVk-i J 

J {There are exactly j of vi^ SVf, , and k—l—j of vi^ ^Vf. } 

..... 



+ 1 



. "'^ )|sup / dai\^ J sup / d(Ti 
3 J I Jvf \ I JVi\vf 



Since da is a probabihty measure, J^s dai < 1, and 



N k-j-1 

dai } < 



fe-2- 



vAvf 



dai 



< 



But 



C]^) <{k- Ip' <{k- we deduce that 



For £ > 0, our lemma follows for C6 < 1, and >> 



C;To 



1. 



Lemma 38 Assume that h, - E L°° satisfy {dt + v ■ Vx + v}h = q(t, x, v), with 
the diffuse boundary condition (184)- Recall the diffusive cycles in (188). Then 
for any < s < f, for almost every x, v, if ti{t, x, v) < s, 

h{t,x,v) = e''^'-*^h{s,x-v{t-s),v)+ / e''^^-^^q{T,x-v{t-T),v)dT; (191) 

J s 

If ti{t,x,v) > s, then for k > 2, 

h(t,x,v)= / e'''^''~*>q(T,x-v(t-T),v)dT-\ -— 



H 
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where H is given by 

k-l 

X] ^{ti>s,ti+i<^}f^i^^^i + - ti)vi,vi)dT,i{s) 



1=1 



, J s 



1=1 

k-l ^ti 



+ X] / '^{ti+^>s}qij,xi-\- {t -ti)vi,vi)dY.i{T)dT 
1=1 

+'^{tk>s}h{tk, Xk,Vk-i)dY,k-i{tk), (192) 

and dYik-i{tk) is evaluated at s — tk of 

dUiis) = {n^tiVi'^^j}{e'^"'^^'"*''*(^Orf'^;}nj=\{e'^("^)(*^+^-*^^da,}. (193) 

Proof. When k = 2, ifti{t,x,v) < s, then (191) is clearly valid. If ti{t,x,v) > 
s, 

h{t,x,v)l{t,>s} = e''^^'>^''-'^h{ti,xi,v) + I e''^^'^^^-'\{T,x + {T-t)v,v)dr. 

Jti 

(194) 

Since '^^'^^^ = e^^g along a trajectory ^ = v,^ = almost everywhere, 
the first term can be expressed (almost everywhere) by the diffuse boundary 
condition (184) and part 4 of Lemma 6 as 

h{ti,xi,vi)w{vi)dai 



w{v) 

^v{v)(h-t) 
W 



Vi 



v)(ti-t) f 

/ Mt,>sM<s}<i"^'''^^"~'''^Ks,xi+vi{s-ti),v^)w{vi)dai 



+ ^-r^ / -^{t,>sM<s}e-''^'"^^^~*'\{T,xi+vi{T-ti),vi)w{vi)daidT 

+ — / la,>4e"(''iH*^-*i>/i(i2,a;i+vi(t2-ii),vi)w(wi)rfai 



^u(v)(ti-t) fli 



+ , / / lit,>s}e''^'"'>^^-''\{T,xi+v,{T-ti),v,)w{vi)da,dT. 
w(v) ' '■■ ^ ^ 



t2 JVl 



Therefore, the formula (192) is valid for k — 2. Assume that (192) is valid for 
k>2, then for A; + 1, we further split the last term in (192) with tk > into 

h{tk,Xk,Vk-i)w{vk-i) = / h{tk,Xk,Vk)w{vk)dak = / l{tk>s,tk+i<s}+i{tk+i>s}- 

For the first term, we further integrate along the characteristics ^ — v, ^ = 
to reach the plane t = s as 

I lt.+i<.<tJe'^('"=)(«-*'=)/i(s,Xfe+(s-tfc)i;fe,i;ft)+ / ' e''^'"''^^^-'''\{T,Xk+{T-tk)vk,vu)dT}w{vu)dak; 

JVk Js 
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For the second term, we integrate along the characteristics ^ = ^ = to 

t = tk+i > s to get 

/ lt,+,>s{e''^'"'^^*''+'-*''^h{tk+i,Xk+i,Vk)+ [ " e''^'"'^^^-*''U{r,Xk+{T-tk)vk,Vk)dT}w{vk)dak. 

•'Vfc •'tfc+l 

We then deduce our lemma by adding Jy^ dak = 1 inside the rest of the terms 
so that all the integrations are over n^^j^V; instead of nJ^J^j^V;. ■ 

Lemma 39 Let ho € L°° and assume (21). There exits a unique solution h{t) = 
G{t)ho e L°° to (29) with the diffuse boundary condition (184) o.^d 

sup {e''''*||{G(i)/io}lt,>o||oo} < e^^ll^^lloc, ||{G(t)/io}lt,<o||oc < ||e-*''('')/io||cx> 

0<t<l w 

(195) 

In particular, 

supe^*||G(0/io||cx> <e'^°max{||^|U,||e-'^(^)+'^°/io||cx>}. (196) 

t>i w 

Proof. Given any m > 1, we first construct a solution to {dt+v-W x+f}h"^ = 0, 
with the following approximate boundary and initial conditions as 

h'^{t,x,v) = [l- —\ [ h'^{t,x,v')w{v')da{x), (197) 
t m J w[v) Jy 

h'^{0,X,v) = hol{\v\<m}- 

Then /i™ = h'^w satisfies {9* + 1; • + iy}h"' = but with 

h"'{t,x,v) = i^l- J hr{t,x,v')da{x). (198) 

Clearly, since J da = 1, this boundary operator maps L°° to with norm 
bounded by 1 — and initially 

||/l™(0)||oo = SUp|/l™(0,X,f;)w| = ||/lol{|t,|<TO}W||oo < CmMl^oWoo < oo. 

Therefore, by Lemma 23, there exists a solution h"^{t,x,v) G L°° to (29) with 
(198), so that we have constructed h™ = ^ with (197), which obviously is 
bounded. Such a solution is unique by the transformation = £ with 

j^\\.r{sw^ds< oo. 

In order to take m — > cx3 in (197), we need to obtain an uniform bound 
(195) and (196) for /i™, which is more delicate. We first claim that it suffices 
to show (195) to derive (196) for /i™. Letting <; = 1 in two parts of (195), since 
w > 1 from (187) and e^^e'^" > 1, we have 

||/i'"(l)||oo <e-^max{||MO)||c,o,||e-^(^)+^°MO)||oo}<e-^||/i(0)||oo. (199) 
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For any I < t < I + 1, we can repeatedly apply (199) to get: 

||/i"(Olloo<e-^||/i'"(Z-l)||oo. 
Therefore, by (195), wo deduce (196) as 

||/i™(t)||oo < e^||/l'"(0||oo<||/l'"(/-l)||oo 

< e-^^||/z'"(l)|U 

< maxlll^ll^, ||e-^W+^°/i(0)||oo} 

< e''«e-'?*max{||^||oo, ||e-''('')+'^°/i(0)||oo}. 

The rest of the proof is devoted to the vahdity of (195). If ti{t, x, v) < 0, we 
have {G(t)hl^}it,x,v) = e-'^('')*/iS^(x - tv,v) and (195) is clearly valid. 

We now consider ti{t, x, v) > 0, then the back-time trajectory first hits the 
boundary. As in the proof of (192) with q = 0, we can ignore powers of the 
factor 1 — -to bound \h™'{t,x,v)\ by 



■X 



k-1 „ 

y\ / l{t,>o,ti+i<o}|/i'"(0,a;,-tm,t;0|£^5]i(0) 

+ [ l{t,>o}\h'^{tk,Xk,Vk-i)\di:k-i{tk)- (200) 

Over the second small set l^^^^o}) choose any e{i'o) > such that 

{l-2^/i)e^ >1, (201) 

then choose ko{s) by Lemma 37 with Tq = 1. By Lemma 37 with Tq = 1, for 

k = fco(e), /n'=-iv, ^{tkX^j^iZidf^i < ^- We further split the second integral in 
(200) into {tk > 0,tfe+i < 0} and {tk+i > 0} in Vk with dcfe. Integrating along 
the characteristic for the first part {tk > 0, tk+i < 0} yields: - — x 

/ l{tfc>o,tfc+i<o}l^™(0,a;fe -Vktk,Vk)\dT,k{0) 

,Xk,Vk-i)\dT^k-i{tk)d(Jk- (202) 

Since ti{tk,Xk,Vk) > over the set {tk+i > 0}, we deduce that 

'i-{tk+i>o}\h"'{tk,Xk,Vk-i)\ < sup|/i"'(tfe,a;,t;)l{ti>o}|. 
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From (193), the exponential in dT,i{s) is bounded by e''"'^" Since fy^ dak = Ij 
by Lemma 37, the last part in (202) is then bounded by 



|/i'"(tfc)lti>o||oo / ^t,>o}wivk-i)Uf-ldai (203) 

{■v 1/2 1/2 
/ l{t,>o}'nf-ldai\ \ [ w\vk-i)dak-i] 

< 2Vi sup {e''°(^-*)||/i™(s)lt,>o||oo}, 

0<s<t<l 

We have used (187) for 6 < Oq ^ \. 

On the other hand, inserting J^^ dak = 1 into the main contribution in (200), 
and combining with the first term in (202) yields: 



/ y] l{t,>o,t,+i<o}|/i'"(0, a;; - *;?;;,«;)! 

"'n* ,v, , -, 



'i=i''' 1=1 

X {n^^^;+ 1 da J }{w{vi)e-''^'"'>*'dai }n^=\ {e''^"' ^''+'~*' ^ daj } 
- J 22hk>0fy+^<0}^^=iW^Vj)daj, 

since w{vj) > 1 by (21) and (187). Note Yl^^i l{<,>o,t,+i<o} = l{tfc+i<o}- By 
(201), we can further choose near 1/4 such that [^] < (1 - 2^)e 2 
in (187). We therefore get 

= / l{tfc+i<o}nf=i«)^(wO<^o-i 
< nf^i 1^ «;'(^;;)dai| < (1 - 2v^)e'^. 
for fc = ko{s). Hence, combining with (203), we have 

snp{e''''%m{t,x,v)^t^>o}\} <^V^ sup {e^°^||/i'"(s)lt,>o||oo}+(l-2Vi)e^' ||^| 

x.-y 0<s<l 'ii' 

Taking suPq<j<i and absorbing the first term on the right hand side, we obtain 

sup e''«*||/i„(i)l{t,>o}||oo < e^ll^lloo- 

0<t<l w 

Letting t = 1, and we deduce (195) uniform in m. We then deduce our lemma 
by letting m — > 00. ■ 
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4.4.2 Continuity of Diffuse G{t). 

Lemma 40 Let n be strictly convex as in (4) and (21) be valid. Let h and 
q satisfy {dt + v ■ Wx + v}h = q{t,x,v), with the diffuse boundary condition 
(184)- Assume h{0,x,v) = ho{x,v), continuous for {x,v) ^ Jq, and q{t,x,v) is 
continuous in the interior of [0, oo] x f2 x with supjo^oojxaxR^ I ^^^y^l < 
Assume 

ho{x,v)\j_ = —— / ho{x,v')w{v')da. (204) 

W{V) Jy 

Then for any t, {x,v) ^ 7q, h{t,x,v) is continuous. 

Proof. Wc fix (t, X, v) such that (x, v) ^ 7q, for any fixed k, wc recall (192) with 
s = for the expression of h{t, x, v). Now we take any point (i, x, v) near (i, a;, v) 
and evaluate h{i,x,v) by (192) with the same number of bounces (fc— bounces), 
and with corresponding t/, xi and Vi and dai. 

Step 1. Reduction to the approximate of phase spaces. Since i^{v) ^ 
\v\ for large v and 4 decays exponentially, by Lemma 39 and the Duhamel 
principle, 

\\h{t)\\^<\\G{t)ho\\oo+ /*l|G(i-s)g(s)|Urfs<C(i,||/io||oo, sup |^|), 

Jo [0,oo]xnxR3 

(205) 

For any e > 0, by Lemma 37 and (205) , we can fix k{e,t) sufficiently large, 
such that the last terms in (192) for both h{t,x,v) and h{i,x,v) are bounded 

by 

{\\h{tk)\\., + \\h{tkm I l{t,^^}Il1lldai + [ l{f,>o}nf-iida( < |. 

For the remaining sets 'i-{tk<0} and l{tj,<o}) for ei « s, define the non-grazing 
sets as 

Vf' = {vi : vi ■ n{xi) > e\ and \vi\ < — }, 

= {vi : vi ■ n{xi) > e\ and \vi\ < — }. 
We further split the integration region in (192) as 

/ l{t.<0} = / l{tfe<0} + / l{t,<0}; 

J J{there exists a Die VAVfi} ^{allDieV, ^} 

h-{h<0} = _ lfe<0}+/ _ l{tfc<0}- 

J J{there exists a w,eVi\V; i} ^{allwieV, !} 

Clearly, by (190), /y^^yei c?cr; + /yj\y=i c?*''; < C'ei, so that from the boundedness 
of ho and ^, the integrals in (192) over the almost grazing sets are small: 



/ l{t,.<o}- 

J {there exists viEVi\Vi } 



+ 



{there exists GVi\V^^"'^ } 



l{tfc<0}- 



<C(/io, J,fc)ei < |. 
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Therefore, it suffices to compare only the integrations over the non-grazing 
sets U^~^Vf^n{tk < 0} and nfj,^ Vf ' n{4 < 0}. For any 1 < I < k-1, recalhng 
the back-time diffuse cycle (188), if a{ti) = {vi ■ n(a;;)}^ > ^ > 0, then from 
convexity and the Velocity Lemma 5, we deduce 

a{ti+i) = {vi ■ n{xi+i)y > Ca{ti) > > 0. 

Hence a;(_|_i, fj+i are smooth functions of xi and vi from Lemma 6. A simple in- 
duction for I implies that xi, ti are smooth functions of (wi, ...vi--\) G Hj^^Vj^ : 

\ti - ti\ + \xi - xi\ ^ (206) 

as {t,x,v) {t,x,v), uniformly in U'fJ^V^ . Clearly nfr^/Vf' C U^'^V^ . 

Since xi ^ xi by the Velocity Lemma 5, Vf^ C V^^ . A simple induction leads 

to (206) and n^=\Vf C n^=\V/ for 1 < I < k - 1. Therefore, (206) is valid 

on both Uf'^Vf' and Tlf'^Vf', subsets of nfj^^V;^. 
Moreover, we have 

\ = {'"I ■ "(^0 >^i^vr n{xi) < ei, and \vi\ < — }. 

El 

Vf' \ Vf' = {vi ■n{xi) < £i,vi ■n{xi) > £i, and \vi\ < —}. 

By continuity (206), for {t, x, v) {t, x, v), xi xi, and both sets are contained 
in ^ 

{si - C\xi - xi\ <vi- n{xi) < ei+ C\xi - xi\, and \vi\ < — } 

£l 

which have measure £hj~j!lA_ ^qw define the approximate phase-spaces as: 

^1 

B = uf-l [Vf ^ n Vf ^] n {tk < 0, tk < 0}. (207) 
To estimate nJ^Tf^Vf ^ \ B, by an induction on k, we get 

\utlvt^\ufil[v!^nvn\ + \^tzlvt^\-ntilvn<c^^^^ sup \xi-xi\. 

l<l<k-l 

Notice that Uf'^Vf' C [Uf-^Vf' \Uf-^V^' ] UUf'^ Vf ^ , Uf'^V^' \B is contained 
in 

nt7 vf ^ \ nt7 [vr n vn u i^tl^r n {tk > o}] 

u[nt7vrn{tfc > 0}] u [nfr>- \ nti^v-] 

From Lemma 37, both /nfr/v,n{t.>o} nfr/da; and /nfj/v,n{t„>o} nfji^dCT; are 
bounded by Cs. By similar splitting for the set Ilfz^^^Vf ^ \ as|a; — xl-j-jt- 
fl -|- |t; — w| — > 0, we deduce 

/ ^i=id<^i+ nfjr/daz < 4C£-FC(£i, fc) sup \xi-xi\<bCe. 

Jnf-iVi^\B Ju'^-^vf^\B l<l<k-l 
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By bound for w in (187) and Cauchy-Schwarz's inequality: 

/ IL'[-ld^i{s)+ [ U'l~ldti{s)<Cit)V^. 

Thanks to the boundedness of ho and ^, to prove the continuity, it suffices to 
estimate the difference of h{t, x, v) and h{t, x, v) in (192), where the integrations 
are over the same set B. 

Step 2. Continuity of h{t, x, v) over B. 

Case 1: ti{t,x,v) < 0. In the case ti < 0, then < by continuity over 
the set B in (207). Then both h{t,x,v) and h{t,x,v) are given by the same 
formula (191). The continuity now follows from {t,x,v) {t,x,v) and the 
continuity of ho and q. Same argument also applies to the situation ti = and 
tl < 0. 

We only need to study the case ti = but fi > in which h{t, x, v) are given 
by the different expression (192). Over the set B, since \vi ■ n{xi)\ > ei > 0, 

from (40) that ti — t2 > But h ^ ti = 0, we therefore deduce that t2 < 0. 
This implies for k large 

B = {ti> 0, t2<o}r\B = nfji^ vf ^ n vf ^ . 

Applying (192) to h{t, x, v) over the set B with fa < 0, by Step 1, we obtain 
h{t,x,v) ~ / e''^^-'^'>q{T,x - v{t - T),v)dT + 

+ — rrr^ / l^i,>OM<o}ho{xi - hvuVi)w{vi)e-'^'''>*' dai (208) 

+ / / %,>o,t-2<o}9(a;i + (T-fi)t;i,vi)^I;(vi)e'^(''i)("-*i)daidr. 

w(v) iv^nv^ Jo 

Since h ^ ti — 0, it follows the last term above is small from the boundedness 
of ^. The first term on the right hand side of (208) tends to 



Jo 



T,X — V{t — T),v)dT, 

as second part of h{t, x, v) in (191), from the continuity of q. Since l{ti>o,t2<o} — 
1 over Ili~^Vi^ n Vf^ in this case, and by ti ^ 0, Xi ^ Xi G dfl, the second 
term on the right hand side of (208) tends to 



■w{v) 



ho{xi,vi)w{vi)dai - e-''('')*/io(x'i, w) e-''^''^*ho{x - tv,v), 



by the continuity of ho away from 7q and the compatibility condition (204). 
Therefore, we have shown h{i,x,v) h{t,x,v) by (191). 
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CASE 2: ti{t,x,v) > 0. From continuity, ii > and 



e''^^-*U{T,x-v{t-T),v)dT ^ / e''^^-''>q{T,x-v{t-T),v)dT. 

Jti 

It thus sufficient to only study integrals (192) over B for both h{t,x,v) and 
h{i,x,v). We further split 

B = J2Bn{U+i < 0,ti > 0;im+i < 0,im > 0} = Y^^im, 

i,m i.m 

and h{t. x, v) — h{t. x, v) ^ J2i m Ib ■ suffices to estimate the difference over 
each Bim, which can be rewritten from (192) as: 



gi/(t))(ti-t) 



■w(v) 



X { ho{xi - Uvi, Vi)di:,i{0) 



+ [ I q{T,Xi + {T-ti)Vi,Vi)d^i{T)dT} (209) 

gi/(i3)(ti-t) r 

X { ho{Xm - tmVm,Vm)dT,m{s) 



m-l 

+ 

3 

ri 



/ X] /- li'^^^j + i'^ -'tj)Vj,Vj)dT.j{T)dT 
JBim j — l Jtj + l 

I I q{r, Xm + {T- tm)Vm,Vm)dT^m{T)dT}. 

Jo JBi^ 



By (40) in Lemma 6, ti — U+i > ^ > 0. For £2 << Si, we further split 
{U > 0,ti+i < 0} = {U > £2,ii+i < -£2}U{0 <ti< e2,ti+i < 0}U{ii > £2, -£2 < k+i < 0}. 

CASE 2a: On the set Bi^, n {U > £2, U+i < -£2}. By continuity (206), for 

{t,x,v) {t,x,v), 

- £2 - £2 

^» > Y'^-^+i - — 2"' 

then h{t, x, v) has the same expression as h{t, x, v) with m = i in (209), and the 
difference in (209) is small over this set, from the continuity of ho and q. 

CASE 2b: On the set Bim (1 {0 < U < £2,U+i < 0}. Now t, ~ U ~ £2 
and ii+i < for £2 << £1. If > and ii+i < 0, we the again have the 
same expression so the difference in (209) is small on this set again, from the 
continuity of ho and q. Otherwise we have U < 0, and as t ^ i, x ^ x, v ^ v, 

ti-i - U-i > ^ > 0, 
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from (40). We define 

Btm = Bim n {0 < ii < e2,ti+i < 0, but ti < 0,ti-l > 0}. 
By the continuity of ho and q, the first term for h{t,x,v) in (209) is close to 

J^^ /io(a;i,i;0{n,t,^W'^i}{«i(^i)'^«^i}n}lUe'^^''^'^'*^+^-*^-^rf^,} 

+ / + E / ?(T,a;, + (r-t,)^;,-,t;,)dS,(T)rfr. (210) 



Since < 0, and > 0, tience m = i — is empty except for m = i + 

The second term for h{t,x,v) in (209) is given by ^— ^j^^^y— ■>< 

J^^ ho{xi-i - ii-m-i, Vi-i){n']zldaj}{w{vi-i)e-''^''*-'^''-'dai-i^^ 

+ 22 <iiT,xj + {T-tj)vj,vj)<n:j{T)dT 

-'bL j=l •'h+i 

+ [ [ q{T,Xi-i + {t -ti-i)vi--i,Vi-i)dT,i-i{T)dT. 
Jo Jb+ 



Since Xi-i - U-iVi-i ~ Xi-i - {U-i - ti)vi-i = Xj ~ Xi, U ~ 0, /J' ' ~ J^*' 
the above is simplified as 



j^^ ho{xi,Vi-^){n';zldaj}{w{vi.^)e^^-^--^^'^-'^^^^ 

+ 22 + - tj)-"j^'>'j)d^ji'^)dT (211) 



By the boundary condition (184), 



ho{xi,Vi-i)w{vi-i) = / ho{xi,Vi)w{vi)dai ~ / ho{xi,Vi)w{vi)dai. 



For ?;i e Vj^' n V-\ we have - U+i > But U < S2 « ei in B^:^, so that 
ti+i < 0. Hence 



n vf^ = {vf n } n {o < < £2, t»+i < 0, and u < 0, > 0}, 

i.e., no restriction on Vi G V^^ fl ^ in Sj^, because ti+i < and ti, ti only 
depends on vi, ...Vi-i, not on from (188). Moreover, since J^ei^^ei dai ~ 1, 
the first term in (211) 

—-— /zo(a;„z;,)^(^0{n,t,^Wa,}n}zUe''(''^-^(*^+^-*^'da,}, 
w{v) Jb+ 
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which is exactly the first term in (210). 

CASE 2c: n {U > £2, -£2 < U+i < 0}. Clearly U > by continuity 
over Bim- If fi+i < 0, then we again have i = m in (209) and the difference 
is small on this set from the continuity of ho and q. We therefore only need to 

consider the set 

BT^ = Bim n {ti > £2, -£2 < U+i < 0, but < ti+i,ti+2 < 0}. 

Since B~^ is empty except for m = i + 1, the contribution for h{t, x, v) is given 
by - — wTTT — X 

^ w{v) 

[ ho{xi+i - U+iVi+i,Vi+i){U';zl^^daj}{w{vi+i)e-''^'''+'^''+'d^ 

^ 111 1- % + ('^ - *i ' (^)'^^ 

+ / / q{T,Xi+i + {t -ii+\)vi+\,Vi+i)dY.i+-i{T)dT. 
Jo Jb7 



Since on -B~„, wc have t-i^i fi+i ~ £2, the last term is small from the 
continuity of q and the second term tends to the second term for h{t, x, v) in 
(209). Since fj+i ~ £2, xj+i ~ xj+i and ij+i ~ ij+i, the first term above takes 
the form 



ji/(i))(ti-t) 
'w{v) 



f _ ;io(x,+i,i;i+i){n^tiV2C^^^}{ii'(^^«+i)d'^On}^i{e''(''^)(*^+i-*^)da,-}. 

J B. 



(212) 

On the other hand, consider the first term for h{t,x,v) in (209). Since 
Xi+i -Xi = -Vi{ti - U+i), and U+i ~ £2, Xi - Uvi - Xi+i ~ £2, by the continuity 
of ho (away from 7g), it takes the form 



e 



v{v){ti-t) 



w{v) 



I ho{xi+r,Vi){n)zl+^dai}{w{vi)e-''(^'^''^dni^^^ 



(213) 

Since fj+i ~ £2 from continuity, for £2 << £1, tij^i < 0, and 

Vfli n V^i n {ti > £2, -£2 < U+I < 0, but < ti+uU+2 < 0} = v^i n v^i^, 

where U+i only depends on Vi, ...Vi, but not on Vi+i by (188). Prom the diffuse 
boundary condition (184) we have 

ho{xi+i,Vi)w{vi) ~ / ho{xi+-i,Vi+i)w{vi+i)dai+i. 
Jvf},nvfLnB- 



i+r "i+r "^im 



Since Lei or- 'i^i+i ~ 1> (213) reduces to (212) and we conclude our 
proof. ■ 
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4.4.3 Decay for Diffuse Reflection U{t) 

Theorem 41 Let Hq € L°° and assume (21). There exits a unique solution to 
both the (23) and (27) with the diffuse boundary condition (184)- Let U{t)ho 
be the solution the the weighted linear Boltzmann equation (27) with the diffuse 
boundary condition, then there exist A > and C > such that the exponential 
decay (14-2) is valid. 

Proof. Once again, with the L°° sohition h{t) = U{t)ha to the weighted Unear 
Boltzmann equation (27), from Lemma 39 and Duhamel principle (30), we de- 
duce from Ukai's trace theorem, h-y e L°° . So f = satisfies the original linear 
Boltzmann equation (23) with / € L^ and \\f{s)\\'^ds < oo. Hence unique- 
ness follows for / by using the standard L^ energy estimate. The well-posedness 
follows from the exact argument in the proof of Theorem 28. By Lemma 29 and 
the L^ decay for the diffuse boundary problem, it suffices to establish the finite 
time estimate (143). 

We apply the double-Duhamel's principle (31). By Lemma 39, the first two 
terms in (31) arc easily bounded by CK{t+ l)e~''°*||/io||cx)- 

We concentrate on the third term 

/ / G(t - si)K-u,G{si - s)K^h{s)dsdsi. (214) 
Jo Jo 

We now fix any point {t,x,v) so that {x,v) ^ 7o. From (192) with q = 0, the 
integrand above is given by 

e-(^)('^-t)l^^^^^^^{K^G{si - s)K^,h{s)}isi,x - tv,v) 



€)(v) — J f-^ '^{ti>si,ti+i<si}{KwG{si - s)K^h{s)}{si,xi + (si - ti)vi,vi)di:i{si) 

gi/(w)(ti-t) r 

H ^z^^ — J l{t^>s^}{G{t - si)K.^G{si - s)K^h{s)}{tk,Xk,Vk-i)dT,k-i{tk)- (215) 



where dEi(si) = {n*z/^^dCTj}e''(''')(«i-*')w(wOdCTin^{e''("^)(*^+i-*i)d<Tj-}, and 
the exponential factor in dS;(si) is bounded by e^°^^^~*^\ By Lemma 39, 

\\{Gitk - si)K^G{si - s)K^h{s)}itk)\\oo < Cxe-'?(*'=-^)||/i(s)|U. 

Letting k = A;o(e,Io) large in Lemma 37, as in (123), the integral of the last 
term is bounded by 

eCKe-"^' sup {e'?^||/i(s)||oo}. (216) 

0<t<To 

To estimate the first and second terms in (215), we first derive the formula 
for KyjG{si — s)Kyjh{s). Recall the back-time cycles of: Xc\{s) = xi — {ti — s)vi. 
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Wc denote fj, = t;/(si,Xci(si),t;'), and a;;, = Xci{ti' , {Xci{si),v'), and GVf, 
for 1 < /' < - 1. By (192) again, 

{KwG{si - s)K^h{s)}{si, Xci{si),vi) 
= J K^{vi,v'){G{si-s)K^h{s)}{si,X,i{si),v')dv' 

= J K^{vl,v')e''^'''^^'-''h^t,^<,y{K^h{s)}{s,Xcl{s,) - (si - s)v',v')dv' 

+ / K^{vi,v') Yl 1*;.^, / K^{vi',v")h{s,x'i, - {i/i, - s)vi,,v")dv"}d^udv' 

J w{v') ' +1- J 

+ / K^{vi,v') — — / lt>^>,{G{si-s)K^his)}{ti,x'k,Vk)dEk-i{ti)dv'. (217) 

J w{v'j J 

Once again, since ||G(si - s)Kyjh{s){t'^,x'^,Vk)\\oo < CKe^^^'''~^^\\h{s)\\ao, the 
last term is bounded by (216) due to Lemma 37, when k > fco(e) large. 

By inserting the main terms in (217) back to (215), we deduce that, up 
to eCif e~"^* supo<t<7'Q{eT*||/i(s)||j^}, the third term (214) in the double- 

duhamel representation (31) is Jg J^^ *dsdsi, where * is 



J K^(y,v')e-(-')i^-^Oi^,^^{K^h{s)}{s,X^i{si) - (si - s)v,v')dv' 

J w{v') ''+1- 

{ j K^{vi.,v")h{s,x'i, - [t'l, - s)vv,v")dv"}di:'i,{s)dv' 



X ^ 



K^{vuv')x (218) 



1=1 



X e"^^'^^'-''^lt.^<,K^{v',v")h{s,X^x{sx) - (si - s)v' ,v")dv" dv' dlli{s^) 

k-i 

j X! ^{ti>suti+^<s^}Kw{vuv')d'>:i X 

J 1 1 1/ 1 



w{v) 

gi/(t>')(ti-si) 

xl{t;,>s,t;,^j<s} r-r-^N Kyj{vi',v")h{s,x'i, - {t'l, - s)vv,v")dv"dY,'i,dsids. 



We now estimate them term by term. For the first term in (218), the back- 
time trajectories never touch the boundary. This term can be easily estimated 
as the proof of Theorem 20 for the in-flow case (e.g. (126)) by 

eCxe-^* sup {e^'\\h{s)\\^} + C,,T, r \\f{8)\\d8. 

0<t<To Jo 

For the second term in (218), we fix /' and consider vi' , v' and v" integration. 
Weiecalld^,{s) = {U']zl,+idaj}e''M{^-t[,)w{vi>)da[,U^^^ 
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and iti{vi')dai/ = • vii}dvii. The exponential factor in rfE;,(s) 

is bounded by e''"*^*!"*'). Notice that with w{vi>)ijl{vi') bounded and integrable 
for ^ < i. Hence from (45), 

/ j \K^{vvy)\dv''w{vv)ii{vv){n{x'i,)-vv}dvv (219) 

J\vii\>N J 

< ^[ w{vi,)iiivi>){n{x'i,) -vi'jdvi' = 

-'V J|^,,|>Ar -'V 

By similar arguments as in Case 1 (119), Case 2 (122) in the proof of Theorem 
20, we only need to consider the case \v\ < N, \v'\ < 2N, and \vi'\ < N and 
\v"\ < 2N, for some large N. As in Case 4 in the proof of Theorem 20, we can 
also use the same approximation (124). Hence, for each fixed we only need 
to control: 

/ If f'\-'°^'-*^\h{s,x[,-{t[,-s)vi,,v'')\dv''dv'dvAu'-I^da'j 

JtT Jn''~^Vi I J\v"\<2N,\v,,\<N Jt',, . . 



^ j Ht[,-s>jr} + Ht[,-s<jr} 
< f" 1{,;,_,>^} + -^ sup {e-^o(«-t;,)||/,(,)||^}. 

Jt'., ' ~ -'^ 0<s<To 



't[,^^ ' 0<s<To 

Here we have used the fact J {^'^Zi> +1^'^ j} = 1- Notice that x'^, and t^, depend 
only on t, x, v, Vi, not on vi' . By making a change of variable y = x'l, — 

{t\i — s)vii, for the first part, we use Fubini's theorem and the fact J {[^jZidcrj} = 
1 to majorize it as (/ = ;^) : 

To ( ~1 ^^'^ To 

Cn,To [ { [ \h{s,y,v")\^dydv" I = CN,To,a [ " \\f{s)\\ds. 

Jo \^Jyen,\v"\<2N J Jo 

(220) 

For the third term in (218), for fixed I, we consider the vi,v', v" integration. 
Recall dEiisi) = {n''rl^^da/}e''^''>'>'^''-^''>w{vi)daiU^jL\{e''^''^^^^^ and 
the exponential factor is bounded by e''°^*^~^'^\ Because w{vi)^i{vi) is boimdcd 
and integrable, as in (219), we can use the arguments in Case 1 (119), Case 2 
(122) in the proof of Theorem 20 to reduce to \vi\ < N, \v'\ < 2N and \v"\ < 3N. 
As in the second term, it then suffices to estimate 



e-'^^*-'^ [ r [' [ \h{s,X^i{si)-{si-s)v',v")\dsdsidv"dv'U^r\daj 

Jn'-\V!, Jti + i Js J\v"\<2N,\v'\<N 



which is bounded by (220) by the change of variable y — Xci(si) — (si — s)v' . 

In the last integral in (218), both back-time diffuse cycles first hit the bound- 
ary. For fixed I and I', we consider the vi,vi' and v',v" integration. We again 
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recall dS; and dSj, and their exponential factors arc bounded by eJ'"'^^^ and 
gi/o(ti-s) respectively. Notice that both w{vi)^i{vi) and 'w{vi/)iJ,{vi') are bounded 
and integrable. We use twice (219), and as in the reduction to (220), we can 
reduce to the case of \vi\ < N, \v'\ < 2N, and < N, \v"\ < 2N. Therefore, 
by using the approximation (124) and J {n.jZi' = 1) I {^jZi+id-f^ j} = 1) 
we only need to control: 

/ / f f' f e''°^'-*^\h{s,x\,-{t[,-s)vi>,v'')\dv''dvi>dsdsi\U'--lda'jU'j^^^^ 

which is again bounded by (220) by the change of variable y = x'l, — (i;, — s)vii. 
Summing over I, I' < k — 1, and letting A'' large, we summarize: 



||/i(t)||cx, <CK{l+t}e-^*||/io||cx>+£CKe-T* sup {e^-||/i(s)||^}+C;v,To,e r°||/( 

0<t<To Jo 

Choosing e small such that eCk = 5, and Tq large so that 2Cif{l+To}e~^-^° = 
e~^'^°, we deduce (143), and by Lemma 29, we deduce our theorem. ■ 

5 Nonlinear Exponential Decay 

Proof, (of Theorem 1): Step 1. Existence and Continuity: Let = 0, we 

use iteration 

irn 

{dt + v\/, + i^- = wTi — , — ) (221) 

w w 

with /i"+i|t=o = ho, = wg. We further split = h^+'^ + 

for TO > 1, where satisfies the homogeneous linear weighted Boltzmann 

equation (27) 

{dt + v-V^ + iy- K^}h^+^ = 0, h^+\- = wg, h'^+^^o = ho; 

while /i™^^ satisfies the inhomogeneous linear weighted Boltzmann equation 
(27) with zero-boundary condition: 

{dt + v-W^ + u- K^}h^+' = wr{ — , — ), h'P+'l^- = /i™+i|t=o = 0. 

WW 

Clearly, from Theorem 20, for some < A < Aq, 

sup e^*||/i^+i(t)||oo <C{||/io||oo+ sup e^«*|H(t)||oo}. 

0<t<CX3 0<t<oo 

On the other hand, denote Uo{t,s) and Go{t,s) to be solution operators for 
linear weighted Boltzmann equation (27) and (29) with zero boundary condition 



%s. 
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respectively, then Uo{t, s) = Uo{t—s, 0) and Go{t, s) = G{t—s, 0) are semigroups. 
Hence, by the Duhamel Principle, 

ft j^m um 

/i™+i= / Uo(t-s)wT{ — , — )(s)ds. (222) 

Jo w w ^ ' 

To avoid the extra weight function v(v) ~ {1 + \v\]'^ in Lemma 9, we further 
use the Duhamal Principle [/o(i— s) =Go{t — s)+ J* Go{t— si)KyjUo{si — s)dsi 
to bound (222) by 

II / Go{t-s)wr{ — ,—){s)ds\\^ + \\ / Goit-si)K^Uo{si-s)wri — ,—){s)dsids\ 

(223) 

For any {t,x,v), by Lemma 9: 

\{njT{ — ,—)}{s,x-{t-s)v,v)\<C,.{v)\m, 

By the explicit formula (109) with g = 0,we therefore can bound the first term 
in (223) as 

I / e-''('')(*-)l,_,,(,,,)<,{«;r( — ,— )}(s,x-(i-s)^;,z;)ds| 
Jo WW 

< C /*e-"('')(*-^)z/(t;)||/i™(s)||^ds 

Jo 

< C [ e-'^^'-'K{v)ds X sup {e-'^(*-"^||/i"(s)||^} 

J 0<s<t 

< Ce-^* X sup{e^«||/i™(s)||^}2. (224) 

s 

Here we have used the fact i'{v) > vo , and / e'~^^*~^''v{v)ds < oo. 

On the other hand, for the second term in (223), let the semigroup U{t)ho 
solve 

{dt+vV.+u- K^^^^^}{U{t)ho} = 0, (225) 



with {C/(t)/io}|7- = and U(0)ho = ho. By a direct computation, ^/l+p\v\^U (t) 
solves the original linear Boltzmann equation (27). By uniqueness in the L°° 
class, we deduce 

U{t)ho ^ ,/lT^Uit){^=^=}. (226) 

Therefore, we can rewrite Ky^U{si, s)wr{!^, ^)(s) to get 

/ / e-''''^'-^^^\\\ K^{v,v')Wl + pm\u{si-s){ ri — ,—){s)}\\^ds^dsdv'. 

Jo Js U ) + p\v'\^ W W 
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Since + \v\f € L^, {-^^=}-2(l + \v\) e so that Theorem 20 is 

vahd for U. Since ."^^ ^ < Cp, from the proof of Lemma 9, 

j Ky,{v,v')^/TTp\^dv' <C j Ky,{v,v'){\v-v'\ + \v\}dv' <+oo. (227) 
Hence the second term in (223) is bounded 

C / e-°(*-''^)^(5i-s)||{ r{—,—){s)}\\^ds,ds 

Jo Js ^1+ p\v'\^ W W 

< C f f e-^^*-''^\\K^{s)\\l^ds^ds 

Jo Js 

< C7e-^*x{sup e^^||/i"(s)||oo}^. (228) 

0<s<oo 

Collecting terms for both and h^~^^, we obtain for < A < Aq, 

sup sup {e^*||/i'"+i(t)|U}<C{||/io||oo+ sup e^°«||5(s)|U}. 

m 0<t<oo 0<s<oo 

for \ \ho\ \ sufficiently small. Moreover, subtracting — /i™ yields: 
{dt + vV. + i^-K^}{h"^+^ - /i™} = w{T{—, — ) - r(^^ , )} (229) 

WW WW 

with zero initial and boundary value. Therefore, splitting 

r(!^, - T{- — , - — ) = T(- — - — , — ) - r(^^ — , —) 

WW WW WW WW 

we can bound ||/i"+i(i) - ^"(Olloo as in (224) and (228): 

/•t um um—1 um 

C\\ Uo{t-s)wT{- ,^^)(s)ds|U (230) 

Jo WW 

ft um—l Lm— 1 um 

+C\\ Uo{t-s)wT{- ,- ^)(s)ds||oo 

Jo WW 

< Csup{e^^{||/i™(s)|U + l|/i'""'(s)l|}} X e-^* x sup{e^«{||/i™(s) - h"'-\s)\\oo}. 

s s 

Hence /i™ is a Cauchy sequence and the limit /i is a desired unique solution. 

Moreover, if fl is strictly convex, by Lemma 19, inductively, ft.™ is continu- 
ous in [0, oo) X {f2 X \ 7o}- It is straightforward and routine to verify that 
wr(^, ^) is continuous in the interior of [0,oo) x f2 x R^. Moreover, from 

Lemma 9, sup | 1 is also finite. We therefore deduce that and 

hence h is also continuous in [0, oo) x {O x R^\7o} from the uniform convergence. 
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Step 2. Uniqueness and Positivity. Assume that there is another 
solution h to the full Boltzmann equation with the same initial and boundary 
condition as h. Assume that supo<t<ru ||/i(t)||cx) is sufficiently small. Then 

{dt + v-V. + u- K^}{h -h} = w{T{-, -) - r(-, -)}, 

w w w w 

so that \\h{t) — h{t)\\oo is bounded by as in step 1: 

C\\ r Uo{t-s)wT{f^,-){s)ds\\^+C\\ f Uo{t - s)wT{-,f^){s)ds 
< Ct, sup {{\\h{s)\\^ + \\h{s)\\}} sup {\\h{s)-h{s)\\^. 

0<s<To 0<s<To 

This implies that supo<4<j,^j \ \h{t)-h{t)\\oo = if supo<j<T'^j ll/^WIU and supo<j<rj, 
sufficiently small. 

Finally, we address the positivity of the F = ^/Jif. We use a different 
approximation for the original Boltzmann equation (1) 

{dt + v- V + = Qgain(F'", F"), (231) 

with the inflow boundary condition = /(x + > and -F'""'"^|t=o = 

fj, + y/Jlfo, starting with = ^. Here 



u{F"')= j F"'{v)\v-u\^dudw. 



In terms of /™ = by (9) and (10), we have 

{dt+vy,+ = Kr + rgain(/™, D - Tiossir, r+'). (232) 

It is routine and standard to show that ft.™ = wf^ is convergent in L°° , locally in 
time [0,To], where Tq depends on ||/io||oo = ll'^^lloo and supo<t<r^ ||.9(i)||oo- 
By induction on m, we can show that if F"^ > 0, the Qgai„ {F"\F"') > 0. 
Denote the integration factor as 

I(t, S) = e- HFnir,Xir),Vir))dr^ (233) 

SO that £{I{t,s)F"'+^} = /(t,s)Qgai„ {F"^,F"') almost everywhere along the 
back-time trajectory [X{s),V{s)] of t,x,v, inside SI. As in (109), we express 
F™+i(i,a;,u) = 



'■t-tb>0 



lt_t^<o{I{t, 0)Fo{x - vt, v)+ I I{t, s)ggain (F" , F™) (X(s) , V {s))ds} 

Jo 

{I{t, t - ib){/x + ^g}{t - ib, ar - vt^, v) + [ I{t, s)Qgain (F™, F'")(X(s), V{s)))ds} 
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so that i^^+i > on [0,ro]. This impUcs that F > in the hmit with h = 
G L°°. By the uniqueness of our solutions with this class, F is the same 
solution we constructed earlier. We obtain F > for all time by repeating 
[0, To], [To, 2ro]...[A:To, (fc + l)To], from the miiform bound of supj ||/i(t)||oo- ■ 
Proof, (of Theorem 2 and Theorem 3): We use the same iteration (221) with 
either bounce-back or specular reflection for h™~^^. By the Duhamel Principle, 

h"'+^ = U{t)ho + / U{t-s)wT{ — ,—)is)ds. (234) 
Jo WW 

Applying either Theorems 28 or 35 respectively, by Lemma 9, we have 

ft i^m um 

||/i"+Hi)lloc <Ce-^*||/io||oc + || U{t-s)wT{^,^){s)ds\\^. 

To avoid the extra weight function i/^v) ~ {1 + \v\}^ in Lemma 9, we use 
U{t — s) — G{t — ,s) + G{t — si)Ku,U{si — s)dsito estimate the second term 
above. For any {t,x,v) ^ 7q, and its back-time cycle (Xci(s), 14i(s)), we use 
(131) and (159) respectively to get 

rt urn um rt um um 

I / G{t-s){wT{ — ,—){s))ds\ = \ / e-'^('')(*-«){^i;r( — ,— )(s,X(s),y(s))}ds| 
Jo WW Jq WW 

< C\ f e-'^^^^*-'^iy{v)\\h"\s)\\l^ds<Ce-'^* X { sup e'^*||/i'"(s)||oo}^ 

Jo 0<s<oo 

where we have used Lemma 9 and (224). 

On the other hand, for the second term, we use U as in (225) with cither the 
bounce-back or specular reflection. Hence (226) still holds. Since w^^(l-|-|i'|)^ G 
L^, Theorems 28 or 35 are valid for U, and we get from (228) 



Jo Js 



um um 

e-oit-^^)\\K^U{si - s)wr{ — , 'J-){s)\\^dsdsi 

w w 



< C 



+ p\v'\^ 



< Ce- = *x{ sup e^^||/i'"(s)||oo}', (235) 

0<s<cx) 

by (227). This implies that sup„ supo<t<oo{e^*||/i"+i(t)||oo} < C||/io||oo for 
1 1 /loll sufficiently small. Moreover, by subtracting h^^^ — h™, by (229) and (230), 
we deduce that /i™ is a Cauchy sequence and the limit h is the desired unique 
solution. If fl is strictly convex, inductively, by Lemma 25 and 32 respectively, 
/i™ is continuous in [0, oo) x {fi x \ 7o}, and «;r(^, ^) is continuous in 

[0,oo) X X R^. Moreover, from Lemma 9, sup | \ is also flnite. We 
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therefore deduce that h is also continuous in [0,oo) x {f2 x \ 7q} from the 
uniform convergence. 

As for the positivity of F, we fohow the argument in the proof of the inflow 

case by using a different iterative scheme (231) with cither bounce-back or specu- 
lar reflection boundary conditions, so that = ^ satisfies (232). Again, it 
follows from a routine procedure to show that /i™ = is a Cauchy sequence, 
local in time [0, To]. Assmne F™ > 0. For any (i, x, v) ^ 7q, we denote its back- 
time cycle (either bounceback or specular) as [Xc\{s), Vc\{s)\, and ti,t2, ■■■ti > 0, 
so that U+i < 0. Recafl (233), so that £{I{t, s)F™+i} = I{t, s)Qgai„ (F™, F"^_) 
almost everywhere along the back-time cycles [X(s),y(s)] of t,x,v, inside Ct. 
As in Lemmas 131 and 159, we can express F"^~^^{t,x,v) as 

I{t,h)F"'+\tuXuVi)+ fl{s,t)Q^^in (F™,F'")(Xei(s),Fci(s)ds > I{t,h)F"'+\tu X^Vi). 
Jtl 

For 1 < i < I, similarly, 

F^'+^tuXuVi) = I{tuU+r)F^+\u+uXi+uVi+i)+ f I{s,ti)Q^^^ {F^ ,F^){X^,{s),Vci{s)ds, 

Jtl 

F-^+^tuXuVi) = /(t;,0)Fo(Xci(0),yei(0))+ / ' /(s, 0)ggai„ (F'", F'")(Xei(s), V;i(s)ds > 0. 

Jo 

Hence by an induction over i, F™+^ > and we deduce F > by the uniqueness 
as in the proof of the in-flow case. ■ 

Proof, of Theorem 4: We use the same iteration (221) together with the 
diffusive boundary condition (184). By (234), we further bound its last term by 
the Duhamal principle as before: 

ft Lm um pt pt im um 

II / G{t-s)wr{ — ,—){s)ds\\^ + \\ / / G{t-si)K^Uisi-s)wT{—,—){s)dsdsi\\^. 
" Jo ^ ' ^ w w Jo J, WW 

We estimate the flrst term above in two parts: 

I / G{t-s){wr{ — ,—){s)}ds\<\ / l + l / I 

Jo WW Jf_i Jo 

For we use estimate (195) and Lemma 9 to get 

I /* ^ e--(*-)l^,,<.){!^^iSlBM}rf,| + I ^* ^ e-W(*-)l{,,<.){^.r(^, ^){s)}ds\ 

< Csupl /* e-'^('')(*-'*V(t;)ds| X ||/i"||^+Ce-'^*sup{e'^^||/i(s)||^}2 

< Ce-^*sup{e^1|Ms)|loo}'- 

For Jq ^, we use (196) to get extra decay in v as: 

/ e-^(*-){| r ^ - 1 - >^ ^ ||^ + ||e-W»;r(^,^)(.)|U}da<Ce-^*sup{e^-||M^)||oor- 

^0 W WW 
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Wc have used the fact i decays exponentially by (187). 

On the other hand, for the second term, we define U again in (225) with 
new weight Wi = r and the new diffuse boundary condition as 

{i}ho}{t,x,v)\^_ = / {Uho}{t,x,v')wi{v')da. 

^1 Jv'-n(x)>0 

Hence (226) still holds. By letting p further sufficienlty small in (21), we can ap- 
ply Theorem 41 for the diffusive C/ exactly as in (235). Hence, sup^ o<t<oo{6^*ll'*'"^^(^)lloo} — 
C'll'iolloo for 1 1 /loll sufficiently small. As before, we can construct the desired 
unique solution and establish continuity when O is convex. 

As for the positivity of F, we follow use the difl[erent iterative scheme (231) 
with diffuse boundary condition 

F™+i(t,a;,u)|^_ =c^ijl{v) j F"'+\t,x,v'){n ■ v'}dv' , 

and = ^ satisfies (232). Again, it follows from a routine procedure that 
ft,™ = w/™ is a Cauchy sequence, local in time [O.To]. Assume F"^ > and 
recall (233) so that ^{/(t, s)F"'+i} = I{t, s)Qg^in {F"\ F™) almost everywhere 
along the back-time trajectory [X{s), V{s)] oft, x, v, inside Ct. Recall (186). For 
any {t,x,v), consider its back-time diffusive cycle {Xci{s),Vci{s)). By a similar 
derivation as in (192), iiti{t,x,v) < 0, 

F"^+^{t,x,v) = I{t,0)Fo{x-vt,v)+ r/(s,0)Qgain {F"',F"'){X^i{s),V^i{s))ds > 0. 

Jo 

Mti{t,x,v) > 0, then for k>2, 



{t, X, v) = f I{s, h)Q^^^ (F", F'")(Xe,(s), V,,{s))ds+I{t, h)c^n{v) I 

Jti Jri';ziVj 



where H"^ is given by 

fc-1 



E l{t,>o,t,+,<o}i^o(^ci(0), Vci(0))ciSr(0) 

k—1 I 

+ E /' l{t, + i>0}Qgain(i^'",F'")(Xel(s),Kl(s))rfi:r(s)<is 

1=1 Jti + i 

+ J2 [' l{t,>0,t,+i<0}Qgain (F'",F™)(Xel(s),Fc,(s))rfSr(s)rfs 

7_1 S 



+l{t,>o}F'"+i(ife, Xk, Vk-i)di:]^_^{tk), 

^i(s) = {Il';zl,da,}{Iis,ti)[nixi) ■ vi], 
any e > 0, by Lemma 37, J^k-2y 'i-{tk-i>o}^jZid'^j < £ for A; large. Notice that 

j = l 3 



where rfSj(s) = {O^.J^^da j}{I{s,ti)[n{xi) ■vi]dvi}Tl]^^{I{tj,tj+i)daj}. For 

Tfc-2, 
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{tk > 0} C {tk-i > 0}, and by (233), I{s,U) < C{snp^^o<s<To I U, To). 
Since Fo>0 and Qgai„ (F™,F'") > 0, we conclude F"^+'^{t , x , v) > 

Ci sup 1 1^^(5)1 loo, To) / l^t^yo}F"'+\tk,Xk,Vk)dvh-iU']Z^daj 

m,0<s<To Jn'^zlVi 

> -Ci sup ||/i'"(s)||oo,To) / l{t,_i>o}|/ {n + ^r+^}{tk,Xk,Vk)dvk-i\ll';zldaj 

m,0<s<To Jn^I?Vi yVk-i J 

> -C( sup ||/l'"(s)||oo,To)£. 

m,0<s<To 

Since e is arbitrary, we deduce that F^+i > and this conchide the positivity 
of F over [0, To]. We then conclude -F > by the uniqueness. ■ 
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